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ABSTRACT 


This report addresses research into active (feedback) control of 
the Satellite Power System (SPS), a large mechanically flexible aerospace 
structure. The basic objective of this analysis is to develop an overall 
assessment of the interaction of the SPS structure and its active 
control systems. 
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Executive SiM«rv; NAS-9-16053 


In order to take maximum advantage of the apace environment, especially 
the available solar energy, future spacecraft and satellites will tend to 
be large (0.1 - 100 kilometers in their largest dimensions) and very mech¬ 
anically flexible structures. It is likely that such large space structures 
(LSS) will exhibit the following characteristics: 

(a) high total mass but low mass density (mass distributed throughout 
the structure over a large area); 

(b) a very large number of low frequency resonances (many below 
1 hz); 

(c) very low natural damping (hence, easily excited resonances) 

Because of their high level of flexibility, they will tend to react to 
forces, such as thermal and gravity gradients, solar pressure, and small 
impulse loadings, in ways which were heretofore igncrable in smaller struc¬ 
tures. In particular, a concept like the Satellite Power System (SPS) 
will be especially vulnerable to these distrubence forces end yet, simul¬ 
taneously, be expected tc meet stringent performance requirements for 
shape fidelity, orientation alignment, vibration suppression, and pointing 
accuracy; however, all LSS will face similar difficulties. 

Active control must be used to augment such structures in order to 
meet these specifications. The basic element of active control is feedback : 
sensors at various critical locations on the structure produce electronic 
signals proportional to local disturbance-induced errors in position, 
velocity, or acceleration; these signals are processed by control algor¬ 
ithms implemented in an on-line digital computer; the computer generates 
new electronic signals which are control commands applied to control 
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actuators such as thrust engines or torquing devices; the control 
actuators, placed at various locations on the structure, exert forces 
to bring the structure hack into line and reduce the effects of any dis¬ 
turbances. Although the placement of appropriate actuators and sensors 
is an important item, the control algorithm , implemented in real-time by 
the on-line computer, is the principle component in this feedback link 
from sensors to actuators; it must be carefully designed to carry out its 
complex tasks within the physical limitations of the actuators and sensors 
and the computational limitations of the computer. 

The essential element of the control algorithm is a model of the 
structural dynamics . This model is used to predict the behavior of the 
structure and compare it vith the actual sensor measurements. It uses 
any errors it detects to improve its perception of the instantaneous 
structural behavior and generate control commands to cause this behavior 
to meet the system specifications, e.g. to keep the SPS microwave system 
accurately pointed at the Earth. For the rather rigid structures which 
comprise current spacecraft ana satellites, such models are relatively 
simple; however, LSS require very complicated sets of partial differential 
equations and corresponding boundary conditions to predict the behavior 
at every location on the structure. A LSS is an example of what is called 
a distributed parameter system (DPS) which in theory must be described by 
an infinite number of modes of vibration, but even in practice a very 
large number of these modes is required to adequately approximate the 
dynamical behavior of the structure. 

The number of nodes included ir. the structure model is related to the 
computer memory capacity and its access time, i.e. how fast the computer 








cub calculate If only a few modes are included in the controller's 
model, the computer can do the calculations very fast and there is 
imperceptible delay in the feedback system. However, if a large number 
of modes is included in order to accurately model the structural behavior, 
the computer must respond more slowly and the actuator feedback commands 
are delayed and arrive too late to correctly reduce the distrubance effects. 
In some cases, they arrive so late that they increase the disturbance 
effects and cause the system to be unstable. This dilemma of how many 
nodes to keep (and which ones to throw out) and still be able to compute 
control co m mands rapidly is called model reduction ; it is the crux of the 
LSS active control problem. Even if all the low frequency modes of a 
structure like SPS were retained in the controller model, they would 
more than exhaust the available computer capacity. Consequently, a model 
reduction must be performed and a trade-off made: Seme modes are retained 
and others (thought to be less critical for performance) are left out of 
the controller's model. 

how we cone to the controller-structure interaction problem. By 
necessity, we are forced to design a controller which actively controls 
the structure using a reduced model. The controller computer is only 
aware of the modes included in its model of the structure. The residual 
(i.e. unmodeled) modes are an intrinsic part of the actual structure, but 
the controller has no internal knowledge of them. When a particular 
sensor detects a local displacement error, this error is made up of contri¬ 
butions from all modes - modeled and residual; the part of the sensor out¬ 
put caused by residual modes is called observation spillover . The control 
computer accepts the sensor signal and compares it with its reduced internal 
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model to produce control commands. Such a command causes a force to he 
applied at an actuator; this force excites all modes - modeled or residual. 
The excitation of the residual modes is called control spillover (control 
energy "spills over" into the residual modes). 

Consequently, through the structure controller we have feedback 
around the modeled modes of the structure (which we want) and feedback 
around the unmodeled modes through the control and observation spillover 
terms. This latter is the controller-structure interaction and it must 
not be ignored in the design and evaluation of active large space structure 
control. It may sound foolishly obvious but an actively controlled struc¬ 
ture is quite different from a passive structure; by electronic feedback, 
the active controller changes the structure's characteristics, and its 
response to disturbances. It is intended that this change will improve the 
performance of the structure; however, if the controller-structure inter¬ 
action assessment is not considered as a fundamental part of any active 
control design for LSS such as SPS, the sought after performance improve¬ 
ments will be lost. In fact, in some cases, the controller interaction 
with the umcdeled modes can cause them to be unstable which would mean an 
eventual less of control of the structure. 

This report presents a "top-down" approach to the design and analysis 
of active control systems for large aerospace structures, like the SPS. 

The introduction in Section I describes the LSS control issues and the 
objectives of this study in greater depth and gives some current references 
on the general problem. Sections II and III describe the main research 
performed: 




v 








(a) development of a mathematical framework for the LSS control 
problem; 

(b) development of methods to synthesize active structure con¬ 
trollers from distributed parameter or finite element models 
of the structure; 

(c) development of methods to analyze and assess the controller- 
structure interaction and its impact on the successful 
operation of actively controlled structures. 

Section IV describes the application of this approach to a generic model 
of the SPS which emphasizes the most stringent SPS control requirement: 
accurately pointing the microwave power transmission sub-system. Section V 
proposes and assesses methods to compensate and reduce the deleterious 
effects of the controller-structure interaction; this discussion includes 
structure, as well as controller, modifications. Section VI discusses 
an important issue for actively controlled structures: the control is only 
as good as the model. We have already seen the effects of the unmodeled 
modes; hovever, even the modeled modes may be ratner poorly described due 
to our lack of knowledge of the zero gravity behavior of large structures. 
Consequently, some of the parameters of the structure may be poorly modeled; 
hence, there is a need for cn-board estimation of structural parameters and 
the ability of the control algorithm to make use of these parameter esti¬ 
mates to adapt (i.e. tune) itself to the structure it is controlling. 
Sections VII and VIII present discussions and conclusions of the results 
of this study and recommendations for future research directions both in 
systems theory and experiment design. The papers published 'under this 
contract appear as Appendices A-S and G; Appendix F contains a paper which 
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I. INTRODUCTION TO THE 8PS ACTIVS CONTROL 


1.0 Introduction 

The Satellite Power System (SPS), as described in [l], is an extremely 
large aerospace structure; its operation will require integration of struc¬ 
tural design and active control systems to generate and transmit, via 
microwaves, large quantities of electrical power from geosynchronous orbit 
to receiving antennae on the earth's surface. The solar collector portion 
cf the satellite is the largest piece of the structure, see [l] Figure U; 
it has a lew structural mass (in proportion tc the total satelxite mass) and 
is mechanically flexible with many lew frequency resonances. In addition, 
the sicrcwave transmitting antenna, although much smaller in site, is a 
flexible substructure which must be fine-pcinted at the earth with minimal 
structural dynamic interaction with the solar erray. 

The following areas are most likely to inquire some active control to 
augment the structural design during in-space operation: 

(a) attitude control for solar orientation cf the array structure; 

(b) shape cr figure control cf the array structure to maintain the 
nominal form; 

(c) fine earth pointing control cf the microwave transmission 
antenna- 

Active (cr feedback) control would be implemented by on-line digital pro¬ 
cessing }f signals .from a variety of sensors located cr. the structure.tc 
produce control actuator commands,for thrust and momentum devices,to make 
the necessary corrections to achieve the above ccntrcl objectives. These 
control tasks might be achieved with separate, cr ''decentralized”, control¬ 
lers cr by a central control processor or seme combination cf these; the 
tc be accomplished simultaneously cr. the S?S. 


tasks would have 








2.0 Technical Discussion 


2.1 Problem Statement 

Sine* the SPS is a mechanically flexible structure, its dynamic be- 
haviour is modeled by a continuum description, i.s. a sat of coupled par¬ 
tial differential equations which model the flexible substructures. «.g. 
beams, plates, membranes, and a set of boundary conditions which model 
the interconnection of these substructures to each ether and to rigid masses 
in the total SPS structure. This means that SPS is in the category of 
distribute d parameter systems (DPS), i.e. those described by partial dif¬ 
ferential equations as distinguished from lumped parameter systems which 
are described by ordinasy differential equations. 

Such a DPS description cf SPS would have the fers: 

■St * ifc * V T (, - 1) 

vhere u(x,t) represents the (vector of) displacements (rotational and 
translational) from the equilibrium position of the structure H at the 
point x and time t due tc the external force distribution ?(x,t). The struc¬ 
tural internal restoring forces A u and damping Du* are generated by partial 

O O w 

differential operators I^with appropriate boundary conditions or. u, e.g. 

A jU + AgU* ■ 0, on boundary of H (2.2) 

The mass distribution of the structure is given by m(x). The force distri¬ 


bution is composed of disturbance forces F_ and control forces P^: 

F(x,t) ■ r D (x,t) ♦ F c (x,t) (2.3) 

Control may be introduced via !•'. force or momentum actuators with 
amplitudes f i (t): 


r b,(x) f 4 (t) 
i»i ■ 


where b 4 (x) are the influence functicns cf the ir.divic.ue_ actuators and 
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f(t) It the vector of actuator commands f<(t). Observation of the system 
it dona via P tensort: 

y(t) ■ CjU ■*» C g u t (2.5) 

where din y(t) ■ P and (2.5) allows for potitIon, velocity, acceleration, or 
combinations of these measurements. There vould alto be tentor and actuator 
(fynanict attociated vith (2.1*)-(2.5) but thete dynamics have not beer, thovn; 
they vould be part of the aaaeetnect. 

The performance desired of the SPS in terms of pointing accuracy, shape 
control, and tolar orientation say be stated in terss of closed-loop tran¬ 
sient ar.d steady-slate behavior or as a performance function, such as modal 
energy or pointing error, to be sinicized. Both types of performance speci¬ 
fication Will be taken Into account ir. this assessment. 

Pres DPS theory, the tynasical behavior of the uncontrolled SPS stric¬ 
ture is described by an infinite number of vibration mode* - the rigid 
body (zero frequency) sedes plus all the independent (resonant frequency) 
flex modes. In practice, finite elenent struct-oral analysis techniques 
are used to produce approximate modal data (mode frequencies and sode shapes) 
which yields a large-scale, lumped parameter model of the SPS (i.e. a l:uge 
but net infinite number cf codes'. This large scale model suffers from 
model error (i.e. error in the modal data) and unmodeled residual modes cf 
the SPS. 

In addition, although the large scale model can be used for analysis 
and computer simulation purposes, it is usually too large to use for an 
on-board, real-time controller. 7r.e limited capacity of such controllers 
will require that a reduced- o rder code- of the SPS be produced; this en¬ 
tails selection of a smaller subset of codes (or an aggregation of modes) 
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that gives a reasonably faithful representation of the total system be¬ 
havior. The active controller design would be based on this reduced-order 
model. Such a controller would take the form: 


f(t) ■ y(t) + H lg z(t) 
i(t) * H gl y(t) + H 22 z(t) 


( 2 . 6 ) 


where din z(t) ■ S < ». The controller dimension S and the parameters 
Hi 2 , Hg 2 would be determined directly from the reduced-order model 
of the SPS. The size of S and the operations involved in (2.6) are directly 
related to the capacity of the on-board control computer (or subcomputer). 


2.2 Fundamental Research Problems for Active Control of SPS 

The most fundamental control research problem for the SPS is the 
interaction of any on-board active controller (2.6), based on a reduced- 
order model, with the actual SPS structure (2.1)-(2.5). This structure- 
controller interaction occurs because of modeling error, control-device 
dynamics, and spillove-* due to modes not included in the reduced-order 
model (for "spillover", see [2]-[3]); this interaction is the most crucial 
research issue to be considered in this analysis. Other important issues to 
be assessed are the following: 

(1) selection of appropriate structure vibration modes to control to 
achieve desired system performance; 

(2) the effect of residual (uncontrolled) modes on system performance; 

(3) the effect of modeling error (i.e. medal approximation) on system 
performance; 

(!) the use of appropriate compensation, if needed, to counteract 
residual mode and model error effects 
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(5) the affect of actuator and sensor dynamics on controller design 
and system performance; 

(6) controller design with multiple distributed actuators and sensors; 

(T) control device (actuator or sensor) location, quantity, and char¬ 
acteristics necessary to achieve the basic control objectives; 

(8) examination of the trade-offs available in the basic control 
approaches, including structural design, passive dampers, direct output 
feedback, and modal control with compensation feedback and 

(9) separation control of the SPS array structure from the microwave 
antenna to reduce the effect of vibrations in the main structure on the 
beam-pointing accuracy. 

2.3 Accroach 

This analysis will use the flexible structure control theory of 
[2]-[6]to develop a general mathematical framework for the assessment of 
basic active control systems procedures for the SPS. The development of 
this framework will be the principle element of this assessment. It will 
be able to begin with a continuum description of SPS, or a finite element 
description, or even a simplified continuum model, such as a plate model, 
and fellow through, in a "tep-dov-" manner, to the basic controller logic 
and provide an assessment of the fundamental issues raised in Section 2.2. 
Particular emphasis will be given to making the framework broad enough to 
include new control trends, e.g. [?]-[9], and, yet, capable of yielding 
answers to the above issues and to the assessment of ether approaches to 
(reduced-order) modeling, e.g. traveling-wave vs. modal descriptions of 
SPS. It will not be limited to the current SPS reference system [l]- 
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This approach will emphasize: 

(a) generic problems of SPS control; 

(b) the basic design logic for controller synthesis; 

(c) the trade-offs inherent in the basic control approaches for the 
attitude, figure 4 and pointing control tasks* 

(d) the interrelationship and interaction of the overall structural 
aesign end its models with the active controller design and operation. 

(e) the distributed parameter nature of the modeling, control and 
interaction problems. 
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II. ACTIVE STRUCTURE CONTROL: CONTOQ TT.gR PTSTftw 
ATO • INTERACTION ANALYSIS 


Our basic objectives in this study are to develop fundamental con¬ 
troller synthesis approaches for SPS structure control and to assess the 
interaction of the SPS structure vlth its active controller(s). An imple¬ 
ment able active structure controller is one which can make use of a vide 
variety of possible sensors and actuators to synthesize on-line . via 
digital computer(s), the feedback control commands necessary to maintain 
the structural shape, attitude, and/or pointing requirements in the pres¬ 
ence of disturbances. Such a controller must be finite-dimensional. The 
controller dimension must not be too high because this dimension corre¬ 
sponds to the on-line computer capacity (i.e., the size of the memory and 
the speed of memory access). 

The fundamental issue is the finite-dimensional . implementable con¬ 
trol of an infinite-dimensional , distributed parameter structure. Any 
controller synthesis technique must produce finits-dime^sionaJ controllers 
for the distributed parameter structure; this will always involve some 
kind of model reduction of the infinite dimensional system. However, due 
tc this model reduction, there will always be a controller-structure in¬ 
teraction with residual (unmodeled) structure dynamics. In general, the 
sources of controller-structure interaction are the following: 

(a) structure modeling error (i.e., poorly known model parameters) 

(b) unmodeled control device dynamics 

(c) spillover due to structure dynamics unknown to the controller. 

Therefore, it is not enough to produce active controller synthesis 

procedures for structural control; the controller-structure interaction 
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mu*t be analyzed for any auch procedure. This analvals la an lntrinalc 
part of the problem and should he done with a view toward controller re¬ 
design for reduced interactions. Without a general mathematical frame¬ 
work which takes the above issues into account, it is impossible to 
handle the active structure control problem in any but a trial and error 
fashion. 

In Paper A: Feedback Control of Dissipative Hyperbolic Distributed 
Parameter Systems with Finite Dimensional Controllers (presented at the 
1979 IDEE Control and Decision Conference; revised and expanded version 
submitted to SIAM Journal of Control and Optimization) we have presented 
such a general mathematical framework for model reduction, controller 
synthesis, and interaction analysis. This framework is, used to assess 
modal control approaches to active structures. In modal control , it is 
assumed that the vibration mode frequencies and mode shapes of the struc¬ 
ture are known. A group of critical modes are chosen, e.g., by their 
excitability by the disturbance frequency spectrum cr by their contribu¬ 
tion to a performance index such as the microwave antenna pointing 
accuracy; the choice cf critical modes is not a trivial issue by any 
means and merits further investigation. These critical modes form the 
reduced-order model cf the structure. The nodal controller is synthe¬ 
sized with only the critical mode data. In Paper A, the most general 
linear controller is considered; this allows trade-offs in the controller 
design, the number and location of control devices (actuators and 
sensors), and the on-line computer capacity. Very simple modal criteria 
are developed to assist the designer in assessing the available design 
trade-offs. The controller-structure interaction is analyzed, and 
bounds on spillover and mod-: error are produced; these bounds indicate 
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the tolerable levels of these interaction tern, for stable operation 
of the active structure. When the hounds ere satisfied it makes sense 
to use the modal controller as designed; when they are not satisfied, 
the designed controller should be used vith caution and, more reasonably, 
it should be redesigned or spillover compensation should be added. We 
will say more later about spillover compensation. The mathematical 
framework presented in Paper A is very general and can handle many other 
types of active structure control; it is not limited to modal control . 

The most ccannon tool available to structural dynamicists is HASTEN 
which is a large-scale computer code based on a finite element analysis 
of a structure. Approximate mode shapes and frequencies are available 
from HASTEN end can be used in active controller synthesis. The fact 
that the modal data is approximate adds the effect of model error to that 
of spillover to further complicate the controller-structure interaction. 

In Paper E: Finite Element Models and Feedback Control of Flexible Aero¬ 
space Structures, the structure control synthesis and analysis problem is 
considered in detail for finite element models of the structure . Bound.- 
on tolerable model error and spillover are produced to indicate when an 
imp!ecen:able active controller based on a consistent finite element 
approximation of a structure like S?S makes sense and can be expected 
to operate stably . 

Among the practical constraints on structure control, we are inclined 
to assume that point actuators and sensors will be used; these may be 
placed at many locations on the structure. In some special situations 
a distributed actuator, e.g., use of electromagnetic field interaction 
of the SPS and the earth, may be available. Such distributed actuators 
and sensors are not excluded from the framework established in Paper A 
and may be used to advantage when they exist. 
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It should be mentioned that in all of our work on active structure 
control ve are talking shout feedback or closed-loon control of the 
structure. In "Control-Structure Interaction in a Free Beam” by S. 
Shrivastava and R. Ried (NASA Report JSC 16699, May 1980), an energy 
approach is used to study the controller-structure interaction for a 
free beam with on-off thruster control. The paper illustrates the 
excitability of structure modes to such nonlinear control forces. This 
is certainly an important consideration; however, one should be careful 
not to draw many detailed conclusions from such an analysis because the 
control is open-loop . Many of the controller-structure interaction 
effects appear because of the feedback nature cf active structure control. 
The analysis of Shrivastava and Ried aids in understanding the complex 
way in which a flexible structure responds to open-loop thruster loads; 
it is quite complementary to our analysis which considers linear forces 
generated by feedback control systems. 

The above mentioned papers A and E are Appendices A and B cf this 
report. 
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III. ACTIVE STRUCTURE CONTROL USING SINGULAR PERTURBATIONS 

In the previous section* we presented s genersl framework for active 
structure control synthesis and analysis. It was based on a regular 
■perturbation viewpoint, i.e., the spillover and model error terms were 
taken as perturbations of the desired system; this is how we arrived 
at our spillover and model error bounds for closed-loop stability. A 
different viewpoint which uses singular (as opposed to regular) perturba¬ 
tions is described and analyzed in Paper C: Beduced-Order Feedback 
Control of Distributed Parameter Systems Via Singular Perturbations 
(presented at the 1°80 Princeton Conference on Information Sciences 
and Systems; revised version submitted to Journal of Mathematical Analy¬ 
sis and Applications) and Paper D: Stability of Distributed Parameter 
Systems with Finite Dimensional Compensators Via Singular Perturbation 
Methods (presented at 1980 International Congress on Applied Systems 
Research and Cybernetics). Applications cf the general approaches to 
aerospace structures are discussed in Paper I: Closed-loop Stability 
cf large Space Structures Via Singular and Regular Perturbation Tech¬ 
niques (presented at 1?80 IEEE Control and Decision Conference). 

The principal difficulty in applying the singular perturbation 

approach is to arrive at a physically meaningful singular formulation. 

When this can be done the system separates into a critical subsystem 

that must be controlled and a residual subsystem whose overall effect 

can be approximated statically. This makes the controller synthesis 

and analysis very different from that of Section II and we feel this 

Justifies a separate section devoted to this approach- • «. . 

‘•v. 

So single approach to the active structure control problem can be 
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expected to yield all the answers . Each one he* ite special advantages 
and disadvantages in a given application. In this and the previous 
section, ve have presented two very general viewpoints which incorpor¬ 
ate the controller synthesis and the interaction analysis which we feel 
is crucial to the solution of this complex problem. 

Papers C, D, and E are Appendices C, D, and E of this report. 


IV. CONTROLLER-STRUCTURE INTERACTION TOR SPS: 
APPLICATION TO A GENERIC MODEL 


Although the assessment approaches described in the previous sec¬ 
tions may be applied to a vide variety of flexible structure models, 
including finite element models, ve have chosen to look carefully at 
their application to a simple model of SPS and the (potentially) most 
difficult control problem: accurately pointing the microwave antennae 
on the flexible solar array structure. Of course, more detailed models 
will exhibit more aspects of this complex problem but ve are operating 
on the trite-but-true principle: You have to start somewhere! 

The generic model ve have used for SPS is a flexible beam with 
lumped masses at both ends; this model vas also considered in 
"Projected Models of Solar Power System" by J. Juang (Jet Propulsion 
Laboratory Engineering Memo 3^7-58, Peb. I960). 

The dynamics are given by: 


l 


u tt + u* + A Q u ■ B q f 
y - C 0 u ♦ E 0 u t 


(!*.l) 


A C U 


a u 


xxxx 


- S u. 


XX 


(U.2) 


with appropriate boundary conditions and u(x, t) is the transverse dis 


placement off equilibrium of the beam of length 
atcr 2 r is chosen later; its form is not of g 
only adds a small amount of dissipation to the 


l . The damping oper- 
reat importance since it 
structure. We consider 


several actuator-sensor arrangements. 




The performance index it given 'by: 


J ■ f [£*(0, t) ♦ E 2 (f, t) ♦ f ? R f] dt (i*.3) 

J 0 

where the pointing error E(x, t) it given by: 

e(x, t) * e(x, t) - e r (x) (fc.M 

where S^fx) it the desired pointing direction at the location x and 


it the angle cf the normal for pointing) vector at location x at tine 
t . The control penalty tern hat a petitive definite weighting matrix 
E ; this penalty tern can be eelected to beep the actuator cc maan da 
within practical licit s. 

The systea (-.1) together with the performance index («.2) fora* 
a distributed Parameter crtinal control tro’cler ; this can be brought 
into the fern cf an optinal regulator by considering: 


e(x, t: ■ u(x, t) - u (x) 
» r 


( 1 . 5 ^ 


where u r (x) is chosen so that.' 


3x 


:a.6) 


It is easily shown that 


{ 


•tt + # t * k o * " s c 1 

y " C C e * £ 0 e t * C C U r 


(*.7) 
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when (m 1* the eu« by proper choica of u r ): 

A oV° 

we shall ehooss to measure 


y ■ y - C Q u r (1*.9) 

instead of y ; this can always l> don* because C c u r is a known tine- 
invariant tern. Also, note that the pointing error becoaes: 

E(x, t)-|| U.1C) 

Combining U.5), (*.7), U.9), (1.10), and (1**3). we obtain a 
distributed parameter optical regulator problem , as desired. This prob¬ 
lem may be analyzed using the methods of ?aper A. We have developed 
active modal controllers by twe straightforward methods: 

(a) produce the infinite dimensional controller (which cannot be 
implemented) and truncate it to obtain the finite dimensional controller; 

(b) truncate the distributed parameter system and performance 
index, and generate the finite dimensional controller from this reduced- 
order model. 

The latter method seems the nest natural from an engineering view¬ 
point. When the exnct modes are known (as they are in this model), the 
two methods yield essentially the same controller. 

This application and the study of its properties form the basis for 
the masters thesis of my research assistant, K. Maahayekhi; this thesis 
will be completed and available by June 1, 1981. 










function of configuration and the damping it a function of material! and 
construction, this may be attainable with atructurtl reconfiguration and 
use of composite materials. Also, we need to know the system parameters 
as well as possible; perhaps, an attempt should be made to structurally 
isolate subsystems or to yield a structure composed of well-approximated 
subsystems. One thing is very clear, the best compromise of structural 
stiffness vs. damping would be made in the preliminary design phase; such 


decisions require a basic understanding cf the effects cf 
figuration and materials tc meet the mission requirements. 


:ructural ccn- 


Structural add-ons include coatings to increase high frequency damp¬ 
ing and (passive or active; member dampers tc incres.se lev frequency 
damping. The latter are especially desirable because it is the low fre¬ 
quency damping that most needs augmentation; perhaps, some combination of 
dampers plus a more sophisticated active control can be used tc achieve 
the structural control requirements. 

Of course, much of the structural design is constrained by factors 
other than the 'ultimate function of the overall structure, e.g.. trans¬ 
portability and construction in orbit, nonetheless, in the acre distant 


ORBINAVfAGE IS 
OP POOR MALITY 













18 


future ve would expect to see control-configured structures . i.e. t 
structure and active control designed and optimized as one integrated 
package. Such an integrated optimization is not simply a matter of 
producing a "super computer program"; a much deeper understanding of 
structural behavior in the presence of active control will be required 
to generate control-configured aerospace structures. 


♦ 


V 





* 7 


2.0 Compensation by Controller Modification 

In general, we feel that structural modification alone is quite 
limited as a means to reduce the controller-structure interaction. The 
interaction problem occurs because of the combination of mechanical 
structure and control electronics; this suggests further use of elec¬ 
tronic compensation as an alternative. 

Actually, there are many ways to introduce this spillover compensa¬ 
tion . The basic requirements are the following: 

(a) compensation should be added-on to the existing controller 

(b) compensation should not substantially increase the number of 
control devices or the dimension of the controller (i.e., the on-line 
computer capacity). 

The methods we have investigated are prefiltering, cancellation, 
residual aggregation, and enhancement. Each of these has its advantages 
and disadvantages; none is a panacea. The most natural seems to be 
prefiltering. i.e., filtering the sensor outputs tc. reduce the spillover 
frequency components; however, phase distortion due to the filter can 
have as deleterious an effect on the stability as the spillover it 
removes. 

We feel this area has not received sufficient attention; there are 
many promising directions for electronic compensation methods which 
need further investigation. 

-•^ Special Purpose (Localized) Controllers 

Decentralized control in general sounds like a good bet for active 


structure control. It has been used successfully in the area, of electric 









power generation and distribution; this application area has some sur¬ 
prising theoretical similarities to flexible aerospace structures. 
Therefore, ve do not want to discount the possible benefits of a de¬ 
centralized approach to structure control; however, ve also do not want 
to forget the extremely difficult technical problems that* continue to 
impede the development of a general theory of decentralized control. 

Still, a localized approach has a great deal of merit if ve exploit 
properties of the particular application. Specifically, the SPS micro¬ 
wave pointing problem might be attacked in the following localized way: 
treat the flexible antenna as the precision system, i.e., the one in need 
of accurate control, and the solar array sls a disturbance system acting 
upon the antenna; try to design a disturbance counteracting control to 
point the antenna and simultaneously reduce disturbances from the solar 
array structure. Such an approach was first considered in "An Active 
Modal Control System Philosophy for a Class of Large Space Structures," 
by M. Balas and J. Canavin (presented at 1977 Blacksburg Symposium). 

Of course, a centralized controller with a complete dynamic model 
of the array structure with antennae attached car. do a much better Job 
than the above-described localized control. However, the centralized 
controller is much more complicated and may prove impractical tc imple¬ 
ment while the more simple localized controller could be implemented 
and may prove to be nearly as effective as the centralized scheme. The 
trade-offs in such localized structure control schemes remain tc be more 
carefully examined. 
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VI. PARAMETER ESTIM/SION ARP ADAPTIVE CONTROL OF STRUCTURES 


Active control of structures is based on a complete knowledge of 
certain system parameters; for exaaqple, in modal control, the mode fre¬ 
quencies and mode shapes must be known (at least, the mode shapes must 
be known at the actuator-sensor locations). In practice, these param¬ 
eters are poorly known and may vary with tine. The control designer has 
two choices: 

(a) make the design robust , i.e., capable of carrying out the 
control task even though the parameters have substantial error; 

(b) make the design adaptive . i.e., use the sensor information to 
estimate poorly known parameters on-line and have the controller tune 
itself up. 

It is difficult to produce robust controllers for highly accurate 
tasks without introducing unacceptably high gains; therefore, we have 
concentrated on the adaptive approach. However, one should certainly 
check the sensitivity of a siven structure controller to parameter inac¬ 
curacies before proceeding with parameter estimation and/cr adaptive 
control. 

A general approach to parameter estimation and adaptive control for 
aerospace structures is presented in Paper F: Toward Adaptive Control 
of Large Structures in Space by M. Balas and C. R. Johnson, Jr. (from 
Applications of Adaptive Control , K. S. fJarendra and R. Monopoli, eds.. 
Academic Press, 19?0). This paper was net written under this contract; 
however, we include it because it discusses the special problems of this 
field within the context of our general mathematical framework. In 
addition, seme preliminary investigation (accomplished with RASA support 
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this year) on adaptive control of reduced-order systems is presented in 
Paper G: Reduced-Order Adaptive Controller Studies by C. R. Johnson, 
Jr., and M. Balas (presented at I960 Joint Automatic Control Confer¬ 
ence ). 

Papers F and G are Appendices F and G of this report. 
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VII. DISCUSSIOH AND CONCLUSIONS 

No one who has worked in the field of aerospace structures would 
argue that active control of a structure as large and unruly as the 
proposed Satellite Power System is anything but a monumentally diffi¬ 
cult tasx. Furthermore, no one but a fool or a madman would undertake 
to solve the basic SPS control problems in a one-year study. Hoping 
not to be judged the former and unsure about the latter, ve have never¬ 
theless undertaken the study of some of these problems under this con¬ 
tract. 

Throughout this work, our emphasis has been on the generic problems 
and the distributed parameter nature of SPS active structure control 
with the goal of obtaining a deeper understanding cf structure control 
technology by focusing on an extremely complicated and difficult structure. 
We have gone in many directions during this study; often, in order to 
cover as much ground as possible, we have developed a given direction 
only to the point where it began to reveal its potential. Still, this 
work is not without substance. We have developed a mathematical frame¬ 
work for the assessment of the active control of aerospace structures 
like ( although net restricted to ) the SPS. Prom this ve have generated 
quite general techniques to: 

(a) synthesize structure controllers from distributed parameter 
or finite element models of the structure; 

(b) analyze the controller-structure interaction on the closed-loop 
stability; 

(c) compensate . by structure or controller modification, to reduce 
the interaction effects. 
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This has been done vithin the constraints of amicus use of exist¬ 
ing control and structure knovledge and tools, allovance for a variety 
of SPS models and configurations, and ability to include new develop¬ 
ments in control and/or structure technology. 

We have accomplished our principal objectives in this study: 

(1) to develop basic active structure control approaches for SPS 
and assess the design trade-offs 

(2) to develop mathematical methods tc assess the effects of con¬ 
troller-structure interaction on the closed-loop stability 

(3) to develop numerical methods to compute the stability bounds 
of (2). 

This work has produced seven published papers and a masters thesis. 
We have presented various aspects of the work in several national con¬ 
ferences and in seminars at L. B. Johnson Space Center and the Jet 
Propulsion Laboratory throughout the year. 

We hope our contributions and the perspective presented here vill 
further the understanding of Satellite Power System control and active 
control of aerospace structures in general. 
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VIII. RECOMMENDATIONS FOR FUTURE RESEARCH IK LARGS AEROSPACE STRUCTURES 

The research reported here has been on the development of systems 
theory. It is In this area that ve shall make our principal recommenda¬ 
tions. 

Two areas, out of many which we have touched on in this report, 
stand out as major issues in the successful development of active struc¬ 
ture control: 

(1) Structural modeling for active control , i.e., modeling of the 
structure explicitly for feedback control purposes as opposed to assess¬ 
ment of dynamical behavior under external loads; 

(2) On-line parameter estimation and adaptive control . 

The first of these really is an indication of the different view¬ 
points of the structural dynamicist and the structural control designer; 
emphasis on this problem would go a long way toward a much-needed 
rapprochement cf these two disciplines. One cannot realistically expect 
tc operate large aerospace structures without some resolution of the 
second issue; yet, on-line parameter estimation and adaptive control 
techniques lead to some of the most complex theoretical issues in dis¬ 
tributed parameter systems. 

There is a clear need to look beyond the systems theory aspects of 
large aerospace structures and to begin to consider the development of 
gr'-und-based, and eventually in-space, experiments in actively controlled 
structures. Of course, no successful ground-based experiment can begin 
to approach the dimensions cf something like the S?S; however, this is 
net necessary to gain the basic understanding that is needed tc cope with 
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this new area of technology. 

Instead, we say begin to investigate particular structural behavior. 
For example, ground-based experiments to determine the dynamical behavior 
of flexible membranes attached to semi-rigid structures would help us to 
understand a major component of the SPS: the Interaction of the solar¬ 
cell blankets with their supporting structure. Furthermore, active 
suspension control could be introduced into such an experiment to inves¬ 
tigate its capacity to suppress the blanket vibrations without seriously 
increasing vibrations in the supporting structure - a clear application 
of the systems concepts developed in this report to design and assess 
the controller-structure interactions. 

Other experiments in active control of beams and plates could be 
designed and carried out to assess the potential of active structures 
and to illuminate the problems. Although computer simulation is an 
excellent tool for analysis of structural behavior, it is only as good 
as the models of the phenomena to be studied. 7nese models require sub¬ 
stantial improvement in order to successfully predict the complex 
behavior cf iarge aerospace structures. Such improvements can only 
grow out cf a better understanding cf structures and their interaction 
with control - an understanding which car. be greatly aided by laboratory, 
and space-shuttle-based, experiments with generic structures. 
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A FEEDBACK control of dissipative hyperbolic distributed parameter systems 

WITH FINITE DIMENSIONAL CONTROLLERS 


Mark J, Balas 

Electrical and Systems Engineering Department 
Rensselaer Polytechnic Institute 
Troy, New York 12181 


ABSTRACT 

The class of dissipative hyperbolic distributee parameter systems 
Is defined here; It Includes many physical phenomena such as mechanically 
flexible structures and certain linear wave propagation problems. Feedback 
control theory for this class of systems is investigated. Although some re¬ 
sults on infinite dimensional controllers are presented, the emphasis is on 
practical (l.e., finite dimensional) controllers based on reduced-order 
models of the distributed parameter system and the closed-loop operation of 
such controllers in the actual system. Sufficient conditions, in terms of 
"spillover'' bounds for the unmodeled residuals, are oresentao for exponen¬ 
tially stable closed-loop operation. 

1.0 INTRODUCTION 

Many physical phenomena are best described by partial differential 
equations, e.o., [1], Chapt. 1, and hence, are distributed parameter systems 
(DPS). Several books have been written on DPS and control, e.g., ll'.-'J], as 
well as, a large number of papers, e.g., the surveys [5’-[3]; applications o* 
DPS control range from mechanically “exible structures to wildfire s*ocres- 
sion, e.g., the survey [10]. Our special area of interest has been applica¬ 
tion to large flexible aerospace structures, e.g., tne sur/eys [!l]-[!2]; 
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although the results obtained In this paper may b« applied to such struc¬ 
tures, they can be applied to a much larger class of problems, namely, 
dissipative hyperbolic DPS. 


2.0 DISTRIBUTED PARAMETER SYSTEM DESCRIPTION: DISSIPATIVE HYPERBOLIC SYSTTMS 


We consider DPS described by the following: 


■ A v(t) + B f(t), v(0) • Vg 

y(t) ■ C v(t) 


(2.D 


where the system state v(t) is in a Hilbert space H with inner product (•,•) 
and associated norm J| *j •. The differential operator A is linear, time- 
invariant, has domain 0(A) dense in H, and generates a C Q semigroup U(t) 
on H. The DPS is called dissipative hyperbolic here when the semigroup 
U(t) has the following "growth* property: 


- *0 # 


t > 3 


(? ? \ 


where z > 0 and Mg > 1; when e > 0, we say the DPS 1s strictly dissipative . 
I* •: » C, we assume M r ■ 1. For mere details on semigroups, see 
[133-CU] - 

The control Is Introduced via M inputs: 


M 

3 f(t) » Y. b i V- t5 
i*1 


ORIGINAL PAGE IF 
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wnere the actuator influence func 
via P sensors wncse outputs y,(t) 


tions are in H. ’he system is ooserved 
form tne vector y(t; • [y-^t).y p (t)]' 
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and art given by 

yj(t) • (c y v(t)) (2.4) 

where the stnsor Influence functions Cj are In H. Th« rank of the Input oper- 
ator B Is M and that of the output operator C is P. In most UPS applications 
the control and observation must be done with a finite number of devices; 
hence, these restrictions on S and C are natural. Usually, the control de¬ 
vices will be localized, l.e., the influence functions bj and Cj will be 
nearly Dirac delta functions. 

The class of dlssloative hyperbolic systems Includes mechanically 
flexible structures, where £ Is related to the natural or materials damping 
In the structure, and certain wave propagation problems which are described 
by symmetric hyperbolic systems of partial differential equations, e.g., [15]; 
In fact. It Includes most oscillatory DPS. It does not Include (so called) 
parabolic problems like the heat equation; for these systems,the semigroup 
U(t) is holomorphic [16] and the control oroolems are somewna: different. 

Control of parabolic DPS nas been considered in, for example, [5^, [loM’?'- 

Our control aporoacn to dissipative hyperbolic DPS emphasizes ^m- 
plementabU1 tv . l.e., from tne sensor cutouts y(t) i" (2.4), the controller 
must synthesize control commands f(t) In (2.3) to adequately staoilize ana 
improve the performance of the DPS, and these commands must be synthesized 
with an on-line computer. Consequently, the control synthesizer must be a 
finite dimensional dynamical system (usually of low dimension). This o’aces 
reasonable, practical constraints on tne DPS control problem [203- 

Infinite dimensional controller results are presentee ‘In Sec. 2, re¬ 
duced order modeling and implementable controllers are ciscussed in Secs. - 
and 5. Dur main results on closed-loco stability *-e presented in Sec. z wi t.o 










Ic 

j ■ 2a 

^ application of these concepts to control of large flexible structures in 
Sec. 7 and general conclusions in Sec. 8. 

ft 

t 

t 

» 

I 

I 

t 
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3.0 INFINITE DIMENSIONAL CONTROLLER 


Although the results presented In this section may not be Implem- 
table (and, in many applications, they are not), they do give some insight 
into the OPS control problem. 

Control laws for OPS, as In the case for lumped parameter systems, 
often take the form of gains multiplying full state feedback: 

f # (t) ■ G v(t) (3.1) 

where G Is a bounded control gain operator; this is especially true when a 
performance index is optimized, e.g., [6], Chapt. 14. However the states 
v(t) for a DPS are in an infinite dimensional space and, in most cases, can 
not be recovered with a finite number of sensors as in (2.4). Consequently, 
one is led to a state estimator (again, as in the lumped parameter case) of 
the form: 


f ;V 


* A v( t) + 8 *( t) + K y( t) • yi t)), v(G) * 0 

y(t) ■ C v(t) 


(3.2) 


where < is a bounded estimator gain ooerator constructed so that v(t) remains 
in D(A). The estimator error, e(t) 2 v(t) - v(t), satisfies 


■sSiil « (a - K C) e(t), e(0) * -v 0 (3.3) 

It is natural to take the control law 

f(t) * 3 v(t) * G v(t) * 3e(t) 

in olace of (3.1) and hooe for tne best. 

The following theorem gives sufficient conditions "or the exponen¬ 
tially stable operation of the infinite dimensional controller (3.2) and 
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(3.4) In closed-loop with the OPS (2.1), (2.3), (2.4): 

THEOREM 3.1 

If (A, B) and (A*. C*) are exponentially stabilizable, then there 
exist bounded operators G and K such that the infinite dimensional 
controller produces an exponentially stable closed-loop system 
and v(t) converges exponentially to the OPS state v(t). 

The proof of this theorem uses results from [21] and appears in Appendix I; 
a OPS (A, B) of the form (2.1) is exponentially stabilizable if a bounded 
gain operator G exists which causes the semigroup U^(t) generated by A + B G 
to satisfy 

~ £ l 1 

! i U 1 (t) |! < M 1 e , t ^ 0 (3.5) 

where M^ >_ 1 and > 0. 

It is known from [22] that, if s * 0 in (2.2), then neither (A, B) 

* * ★ 
nor (A , C ) can be exponentially stabilizable because 3 and C are finite 

rank (and, hence, compact). This means tnat, if the open-loop DPS is not 

strictly dissipative (no matter how small t is), then no amount of controller 

feedback will make the closed-loop system exponentially stable. 

The infinite dimensional controller may also be written in convolution 

form: 

t - 

f(t) * / GU (t - t) K y(-) dT (3.5) 

0 

where U(t) is the Cq semigroup generated by A - 3 G - < C. However, it may 
not be possible to implement either (3.6) or (3.2) and (3.-) with finite 
dimensional devices. Cf course, either version of the infinite cimensiona! 
controller may be approximately implemented and this is the subject of the 
next section. 
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4.0 REDUCED-0R0ER MOOELS 0 F OPS 


Even though the underlying structure of a OPS nay be described by 
(2.1)-(2.4) and, hence, is infinite dimensional, the controller may have to 
be based on a finite dimensional approximation of the original OPS, in order 
to be Implementable in the sense of Section 2.0. Such an approximation or 
reduced-order model (ROM) is a projection (not necessarily orthogonal) of 
the OPS on H onto an appropriate finite dimensional subspace of H. The 
ROM subspace has dimension N and its projection is denoted by P N ; the 
residual subspace H R associated with completes the decomposition 
H * an< * its projection is denoted by P R . The total OPS state v 

is the sum of the ROM state v^ = v and the residual state v R = P R v: 

v(t) » v N (t) -r v R (t) (4.1) 


The choice of subspace and the projection P^ completely specifies 
the residual subspace H R and its corresponding projection P R ; this choice 
is usually cictated Py the physical aoolication and knowledge (or lack of 
knowledge) about the scecific DPS involved. In many situations, the partial 
differential eouations or the corresponding boundary conditions are too com¬ 
plicated to permit simple closed-form solutions. Consequently, an acoroxi- 
mation of DPS, such as the Rayleigh-Ritz-Galerkin or finite element method, 
[23]-[24], must be used to deal with the problem numerically; this is one way 
to produce a ROM of the OPS but it is by no means the only way to generate 
such approximations. 

We will restrict our attention in this paper to consideration of 
ROM's obtained from orojection (not necessarily orthogonal) onto a pair of 
reducina subsoaces H„, H„ for the differential operator A, i.e., the following 

' ii K 
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conditions hold: 


H - h n 0h r 

(4.2a) 

P„ (0(A)) 5 0(A) 

(4.2b) 

p N A P R v . 0 

(4.2c) 

P R A P N v • ° 

(4.2d) 


From [25], p. 172, these conditions are equivalent to the condition that 


P^ commutes with A: 


v e D(A) ♦ Pjj v e 0(A) 


A P N v » P N A v 


(4.3a) 

(4.3b) 


Of course, P R also comnutes with A. 

Throughout the rest of this paper, any ROM considered for (2.1) 
will be based on some pair of reducing subsoaces and_ H R for A . The 
projection onto these subspaces may or may not be orthogonal; when it is 
orthogonal, we will have H R = and 


V 1 S * + 

V * 


v R :; 2 


(4.4) 


ana 


i .e. 


i i p ] 
i ‘ K N 1 


R ; 1 


; 4.55 


An example of reducing subspaces occurs when H„ is a modal subsoace , 


* span {3,, . . . , 3 n ; 

wnere A for eigenvalues with correspcnaing eigenfunctions 
these eigenvalues must be separated from the rest of the spectrum of 
(i.e., a rectifiable, simple closed curve can be drawn in the complex p 


V 

A 

lane 
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so as to enclose an open set containing . . . , \ N in its interior and 
the rest of the spectrum of A in its exterior). From [25], p. 178, and 
H r form a pair of reducing subspaces for A and s A is a finite 

dimensional operator whose spectrum is exactly {X^. X N > while the 

spectrum of A R s P R A P R is the rest of the spectrum of A. Note that the 
projections are orthogonal when A is normal [25], p. 277. 

We emphasize that other reducing subspaces may exist for A but 
the modal one is useful when a finite set of separable eigenvalues, and their 
corresponding eigenvectors, can be found; this is the case, for example, when 
A has compact resolvent ([25], p. 187). Also, the modal subspace is physi¬ 
cally reasonable for producing ROMs of many DPS such as mechanically flexible 
structures. It should be noted that, when the exact modes are not immediately 
available (even though they may be known to exist) and approximate modes are 
used for H^, the subspaces and H R are not necessarily reducing subspaces. 

We have the following result for dissipative hyoerbolic DPS: 

ORIGINAL PAGE IS 

THEOREM 1.1 OF POOR QUALITY 

If and H R are reducing subspaces for A in (2.1), then this 

OPS decomposes into 

v M (t) * A n v N (t) + 3 n f(t) , v N (0) * P N v Q (A.6) 

3vp(t) 

—~— 3 A r v r + B r f , v R (0) 3 P R Vq (4.7) 

y(t) * c N v N (t) + C R v-(t) (4.3) 

where v N 3 P N v, v R « ? R v, A^ * A ? N , 3 ? , 3 ? N 3, C N 3 C ? N , etc. 

In addition, the projections and P R commute with the semigroup 

U(t) generated by A and A R generates the Cq semigroup U R (t) = P R U(t) P R 






9 




i 

I c 

I 


I 



I 

i 


; C 

§ 


r 


i 


r 


r 


R 


f 


with the growth property: 

-et 

l!U R (t)|| < M R e , t > 0 (4.9) 

where ^ s ||P R || M Q with^, e)as in (2.2). 

The proof of this theorem uses results in [25] and appears in Appendix II. 

The ROM corresponding to the projection of (2.1) onto the reducing 
subspaces and H R is defined by (4.6) and (4.8) with C R assumed zero: 

v N (t) » A, v N (t) + B n f(t) , v N (0) » P N v Q 

y(t) * c N v N (t) (4.10) 

or (4.10) may be abbreviated as (A^, B^, C^). In particular, when is the 
modal subspaca for N(separated) eigenvalues of A, we call the corresponding 
ROM - the modal ROM . 

The interaction terms B R f and C R v R are called control and observa¬ 
tion spillover , respectively. It is through these terms that a feedback con¬ 
troller can interact with the residual subsystem (4.7). This would not be the on 
interaction if the ROM had been based on nonreducing subspacas for A; in that 
case model error terms - A v^ and P R A P^ v^ - would have appeared in 
(4.6) and (4.7), respectively. However, these terms are zero when reducing 
subspaces are used for the ROM (as pointed out in Theorem 4.1). 

Henceforth, the ROM (A^, B^, C^) will be assumed to be controllable 
and observable . It is not difficult to verify whether a given ROM satisfies 
this assumption via the usual finite dimensional space technicues; when the 
modal ROM is used, it is especially easy to verify, e.g., [12]. 

5.0 IMPLEMENTABLE FEEDBACK CONTROLLERS g 0R OPS 

The performance of the controlled DPS will usually be expressed in 












10 


terms that can be accomplished with ROM state feedback, e.g., pole placement 
or optimization of a performance index for the ROM. For example, when the 
modal ROM is used, the modes chosen would be the most critical (i.e., easily 
excited by whatever disturbance sources are present) modes ; Increased 
damping and/or stiffening of these inodes would be the main goal of feedback 
control. We will suppose that the desired control law for the OPS (2.1) 
is given (or adecuately approximated) by 

f*(t) « Gjj v (N (t) (5.1) 

where the gain operator G N can be selected via pole placement or optimal 
regulator methods, [26], applied to the ROM (4.10). 

If this control law (5.1) could be implemented, the closed-loop 
system (2.1) and (5.1) would have the form (using (4.6)-(4.7)): 


v N (t) * (A n + B n G n ) 

v„(t) 

(5.2) 

3v a (t) 



St * 3 R G N 

1 + A r v R (t) 

(5.3) 

where v(t) * v^(t) + v^(t). 

The closed-loop stability of (5.2)-(5.3) 

would 


follow from the result: 


fHEORSM 5.1 


Let the control law be given by (5.1). If (A^, 3^) is controllable, 


then Gjj can be chosen so that - e > 0 and 


! v( t)! | < m q m n m e 


'V 


's.i) 


where 

( . 


R : 


' P ' I : D 1 

1 k R' ■’ N 


if nonorthogcnal projections P^ and are used 

1/2 


(1 + 2 ^ + ) 

^ if orthogonal projections are used 


ORIGINAL PAGK i. 
0? POOR QUALITY 
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and ctjj 5 8 /Ajj, 3 s |jB R G^H (the control spillover coefficient ), 

Sut 

A n s A n + B n G n generates e with 


'N 


.^N*, 


<"n 




t > 0 


» w — 


(5.5) 


and (^, e) are as In (2.2). 

The proof of this theorem appears In Appendix III. Note that only control 
spillover 3 Is present in (5.4); this is because direct access to v„(t) has 

ll 

been assumed. 

Of course, one rarely has direct access to the ROM states v N (t); at- 
best, one could try to recover them from the outputs of P * N sensors but 
observation spillover would interfere with this approach. Therefore, the 
control law (5.1) can not be exactly implemented (except in the special case 
that the observation spillover in (4.8) is zero). Instead, estimates of 
G^j v N (t) must be generated from the sensor outputs y(t). This can be done 
via the following imolementable, finite dimensional controller : 

f(t) » y(t) * H 12 z(t) (5.6a) 

z(t) « F z(t) * H y(t) + £ f(t),.z(G) « 0 (5.5b) 

where dim z * S £ N * dim v^. 

Let e N (t) 5 z(t) - T n v N (t) and, from (4.6), (4.8) and (5.5b), we 

obtain: 


e N (t) « F e N (t) + H C R v R (t) 

, 4 N !0! • ? N v 0 

when the following assumptions a*** made about and £. 
F T N ’ \N A N * H C N * 0 


f ; ?\ 


(5.3) 


'N S N 


,0.3] 
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\k will assume that the matrix F can have anv desired stability margin 
a F , i.e., 

lie II < M F e F , t > 0 


(5.10) 


This means that P can be chosen so that It shares no coRmon eigenvalues with 
A n ; therefore, by [27], Chapt. 8, there will exist a unique solution T N of 
(5.8) and hence, E can be chosen to satisfy (5.9). 

From (5.6a) and the definition of e^t), we obtain: 


f(t) ■ Gj, v N (t) + H 12 e N (t) * H n C R v R (t) (5.11 

when we make the assumption that and can be chosen to satisfy: 


H 11 C N + H 12 t n “ g n 


(5.12) 


This last assumption, plus (5.10), can be satisfied when N - P<_S<,N, e.g., 
[28], but it can often be satisfied for S < N - P, e.g., [29]; the special 
case S * N is the usual full state estimator or Kalman filter based on the 
SOM, e.g., [30]. Therefore, the assumotions (5.10) and (5.12) can always be 
satisfied by some implementable controller (5.6) whose dimension S is less 
than N; the actual design of such a controller is a well-studied problem in 
finite dimensional space and will not concern us further; such a controller 
may be synthesized directly with a digital comcuter, [26], Chaot. 6. However, 
please note that the design of the controller (5.6),and all necessary assump¬ 
tions,can be done entirely in terms of the ROM. This is a reasonable aporoach 
since the residual subsystem parameters (A R , 3 R , C R ) may be much less accurately 
known than those for the ROM; in most cases only bounds on these parameters 
would be available in practice. 

The question remains: How well can this i mo''emen table controller 


be made to perform? This is the subject of the next section. 
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6.0 CLOSED-LOOP OPERATION OF THE IMPLEMENTABLE ROM-BASED CONTROLLER 


Even though the design of the controller (5.6) Is based entirely 
on the ROM(4.10), this controller must operate In closed-loop with the 
actual OPS (4.6)-{4.8). Consequently, the residual subsystem (4.7) Is In¬ 
volved In this operation via the control and observation spillover terns. 

It Is a rare situation when these spillover terms are zero; this occurs If 
the actuator and sensor Influence functions b^ and c^ lie entirely in the 
ROM subspace Hjj. Most practical DPS problems will Involve a certain amount 
of spillover even though It may be small; this can lead to degraded perform¬ 
ance and, In some Instances, Instability [30]. 

Let q(t) s [vj(t) e^(t)] T . From (4.6)-(4.8), (5.7) and (5.11), we 
obtain the closed-loop system : 


q(t) » q(t) + X| 2 v R (t) 




v„(t) 


3V R 


rz— ■ X,, q(t) + 

L i 


22 


wher 




12 


A N + S N 3 N H 12 


®N H 11 C R 
H C- 


( 6 . 1 ) 


^21 * ^ 3 R S B P. ^12^ ’ and ^22 * * 3 R w R 

with the initial conditions Cg ■ [ 3 4J Vg, -’ 4J Vg]‘ ar.a /-!:) « 3e* 

!!fl(t)n • C'v^t)!' 2 * ;;«„(*.) :: 2 ) V2 
have the following useful result: 
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THEOREM 6.1 

If no observation spillover Is present (C R ■ 0) In (6.1), then F 
can be chosen so that A c * c p - o N > 0 and 

i n(t)n*(iiv N (t)n 1 +ii. N (t)if)’‘< M e 

where 

!TJ 




q(0) M . t > 0 (6.3) 


M c - *N M F (1 + 


iw 


(6.4) 


and (Mjj, a^), (M p , o p ) are as In (5.5) and (5.10), respectively. 

The proof of this result appears In Appendix IV. Note that (6.4) is not the 
only bound for M £ ; in fact, the "sharper" we can make the bound on M c , the 
better our later results become. 

The next two results are the principal ones In this paper; they give 
sufficient conditions for exponential stability of the closed-loop system in 
terms of spillover bounds: 


THEOREM 6.2 


Assume that (A^, 3^, C N ) is controllable anc observable, that 
is chosen in (5.12) so that ^ - •: > 0, and that F is chosen 

in (5.6b) so that j p > If no observation spillover is present 
(i.e., ■ 0), then, the closed-loop system (6.1) satisfies: 


* 


>(c)’' < M. M m e 

“ w J 


t > 0 


15 . 5 / 


where 







IS 


«(t) i Cv N T (t) .„ T (t) v R T (t» T . 

v(t) • v N (t) ♦ v R (t) , 

IMOII » (||v(t)|| 2 + ||. N (t)|| 2 ) 1/2 md 
M* Is given by (5.4) In Theorem 5.1,and Is bounded as In (6.4) 

In Theorem 6.1. 

The proof of this theorem appears in Appendix V. 

THEOREM 6.3 

Let the open-loop system (2.1) be a strictly dissipative hyperbolic 
OPS (l.e., e > 0). Under the same assumptions as Theorem 6.2 except 
that observation spillover Is present , the closed-loop system (5.1) 
stability margin is at least e ^here 

r (6-6) 

when r 5 M3!! i !H^ C R !' + | |H C R * | is sufficiently small £ / > 0. 

In that case, 

✓ 

' >(t)H < M c M 0 M* e‘ £ : "w y ' : , t>0 £«./) 

The proof of this theorem appears in Appendix 71. There is no ooint in con¬ 
sidering such a result when (2.1) is not strictly dissipative hyperbolic,be¬ 
cause s' could not ever by positive in (6.6). 

7.3 APPLICATION: CONTROL OF r l£XI3Lc $T"UC~ ! JF.£S 

Peecbacx control of vibrations in mechanically flexible structures 
nas acoIication to many cur-ent engineering orcolems. Control of flexible 
structures, in general, was addressed in [201, [2t] ana [28]. A new aoplica- 
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tlon area - control of large aerospace structures - Is very active at present; 
the current control theory trends in this area have been surveyed In [12]. 

The results obtained In this paper have direct application to control 
of flexible structures, In general, and largo aerospace structures, In parti* 
cular; this application will be discussed In this section In order to Illustrate 
the use of the results of Sec. 6.0. This section will follow the format 
established in [38]; for clarity's sake, some of the format will be repeated 
here. 


We consider the class of flexible systems that can be described by a 
generalized wave equation; this represents the Idealization of many mechanically 
flexible structures. The generalized wave euqation Is given by: 


3^u 


+ + a u 

at 0 


(5.1 


which relates the vector of displacements u(x,t) of a body a bounded open 
set with smooth boundary in n*dimensional Euclidean space 3 n , to the applied 
control forces F(x,t). The operator A is a time-invariant, symmetric differ¬ 
ential operator with compact resolvent and lower semi bounded spectrum, "he 
domain 0(A ) of A Q is dense in the Hilbert space L 2 C.) with (•,•) denoting the 
usual inner product and '' • 1 denoting tne associated norm, "he natural camping 
in the structure is modeled bv tne term br-r where 0 is a non-negative constant. 

* 9L 

The control forces 


M 

F(x,tl ■ 3f(t) - : b^ (x)f t) 
1«1 

are provided by M actuators with influence 
are measured by ? averaging sensors 
y(t) ■ C fl u(x,t) 

•where y 4 (t) ■ c,(x)u(x,t)dx with j ■ 

j •• J 

sensor functions b,.(x), c^x; are in L’v.T/ 


‘/RIGDSAL PAG. 1 .' 

OP POOR QUARTS (5.2 

functions b^(x). "he disolacements 

,2.“. The actuator and 

and normalizes to nave unit integral. 
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When the support of b^x) is in a small neighborhood of a point x.., we say it 
is a point actuator and, similarly, we define a point sensor . The point 
actuator and point sensor situation is of special interest here. This class 
of distributed parameter systems includes interior and boundary control of 
vibrating strings, membranes, thin beams and thin plates. Although, only 
displacement sensors are considered in (6.3), this is not a restriction - more 
general sensors may be modeled. 

It is well known ([25], p. 277) that the spectrum of A contains only 
isolated eigenvalues X k with corresponding eigenfunctions such that 

x 1 < x 2 < . . . 

and Ao (< = x, We will assume that X^ is positive. The resonant mode fre- 
quencies of the structure are given by = (X^) and the corresponding 
eigenfunctions ^ are the mode shapes . Thus A satisfies 

1 a||u|i 2 , a > 0 

2 

Every vector usL (P.) has a unique representation 


where u. = .Vuc.dx and we define the orthogonal projections ? ' , p ' by 

.< ’..X u u 

N 

a £ u.p. 


A new operator A is 
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The energy inner product (•,•)£ is defined on by 





+ (w^ § w 2 } 


( 6 . 8 ) 


for u 2 e V and w.,, w 2 e W , and the Hilbert space H is defined as the 
closure of in this energy inner product. The associated energy norm is 
denoted by Jj*j|g and is a measure of the total potential and kinetic energy 
in (u,|y) where u is a solution of (6.1). Let v * [u 1 *, |^-] be in H and 
write (6.1), (6.2) and (6.3) as 


f " Av ♦ Bf , v Q cH 1 
y * Cv 

where B J S and C * [ Z n 0]; this is in the form (2.1) 

l°oJ 0 

It is easy to see that 


(6.9) 


v * 


k*l 



(5.10) 


n ? . o 

and ;iv !i t ■ Z [\.u. “ + u, ]. Also, A generates a C semigroup 'J(t) 

" k»l * K * 3 

by [1 A] p. 11-59. The following result, due to Goldstein and Sosencrans 
[39], gives conditions under which this semigrouo U(t) satisfies (2.2) and 
hence (6.9) is a (strictly) dissipative hyperbolic system: 
if b < 2a then 


( 5 . 11 ) 


v(t) 


< Y e‘ ot ;'v(O) 


t > 0 


where y * [l-ib/Ca)][l-(b/2a)]“ 1 . Consequently, 'J(t) satisfies (2.2) w’th 


M « ( Y ) 1/2 and c * b/2. Nore that the hypothesis for this result means 
that e < a, i.e. the damping must be smaller than the souare c* the lowest 
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mode frequency . Another bound of the type (6.11) is available in [36] without 
the restriction e < a; however, the damping coefficient b is not quite so 
directly related to s in (2.2). The bound used above will be very convenient 
to illustrate our results. 

Recalling (6.6), we define the projections P N> P R on H by 


_ 


r R N l 

u 

_ 

p o u 



N 

.w. 


o 

a. 


' 


r r i 

“u 

- 

p o u 



R 

_W. 


Lp 0 w J 




( 6 . 12 ) 


and note that they are orthogonal in the energy inner product and that P^H and 
H R *P R H form reducing subspaces of H, as described in Sec. 4.0. In fact, is a modal 
subspace, and v N »P^v and v R »P R v satisfy (4.1), (4.4)-(4.5). Consequently, the 
system (6.9) partitions into the form (4.6)-(4.8) ; and the corresponding modal 
ROM is given by (4.10) where 3^, may be identified with the matrices 




‘N 

-bl 


[C N ° 0] , respectively, with A,, * diag [\ 1 




fnus, the controllability and 

r 


v N J> 8 n 3 . and C N = U^,o k )] 

observability of the ROM can be determined from (A^, C^ u ) wnen the damping 
b is small ([38], Theorem 2.1). Of course, it is quite easy to assess the con¬ 
trollability and observability of this latter system and the conditions deoend 
on the interaction of the actuator-sensor influence functions with the mode 
shapes o,, . . . ,3^. 

If the ROM is controllacle and observable, the usual ~lnite d~mens 1 era’ 


state space controller can be designed based on the medal ROM (e.g. [26]). 
has the form: 








(6.13) 


f(t) * G,, v N (t) 

v M (t) * a n + b n + ^(y(t) - y(0) 

(6.14) 


|y(t) » c N v M (t) , v N (o) - o 

This is a particular case of (5.6) with dimension S«2N (twice the number of 
modes in the ROM) and parameters: H^*0, H 12** * * A N ” H *S an< * 

£«B n . Also, note that e N (t) * v N (t) - v N (t), i.e. I N ■ I 2N - Since ROM is 
controllable and observable, and can be chosen to achieve (5.5) and 
(5.10) with any desired stability margins 0 ^ and dp. 

Thus, (6.13)-(6.14) represents an imolementable (finite dimensional) 
modal controller for the flexible structure (5.1)-(5.3), or equivalently (6.9). 
It is based on a modal RGH, i.e. it is designed to control the first N modes 
of the structure (of course, different sets of N modes could have been used - 
it is not necessary that they be the first N modes). The basic assumption 
made in synthesizing such a modal controller is that the N modes chosen for 
control will adeouatelv , *eoresent the important vibrational behavior of the 
whole flexible structure ; this assumption can often be satisfied in oractice, 
and it is the basis for many structure controllers. 

Having synthesized an implementable modal controller to adequately 
tailor the vibrational response of N critical modes in the flexible structure, 
we must not think we are done. The closed-loco benavior of our controller is 
affected bv the resicual modes in H„ for which we have not designed, as well 

" ft 

as the controlled mcces in for which we have designed. These interactions 
occur through the control and observation soil'cver terms defined in Sec. -.3 
and can drastically altar the ces^red performance of the controlled structure; 
in fact, an example is presented in [30] where the closed-loop system become' 


unstable. 
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Consequently, the analysis of closed-loop behavior discussed in 
Sec. 6.0 is especially pertinent. The principal results Theorems 6.2-6.3 
give an indication of how much spillover can be tolerated in the closed- 
loop system. In the context of flexible structures. Theorem 6.2 may be 
compared with the results of [38] where observation spillover was assumed 
negligible (however, note that no damping was present in [38]). Combining 
the results of Theorems 6.2 and 6.3 for this application, we find that the 
closed-looped system,consisting of the flexible structure (6.1)-(6.3) and 
the finite dimensional modal controller (6.13)-(6.14) J will remain ex¬ 
ponentially stable with stability margin s * given by (6.6),if the observa- 

/ 

tion spillover coefficient r is sufficiently small to make e >0, or 
equivalently: 

r < £/m c M o M * (6.15) 

In this application, e * j where b<2a in (6.1) and (6.4), M Q is given by 


M s 

'o 


2a^b 

2a-b 


1/2 


Also, from (5.4) since orthogonal projections are used: 


M * (1 + - - N “) 


1/2 


where ^ * 3/4^ and 4^ * z^ -£, and, from (6.4), 


M c i -Vf (1 * -4-^ 


,0 . 10 ) 


. C . i / 


( 6 . 18 ) 


where 1 * c- - and i^, Mp come from (5.5) and (5.1C) respectively. 
The '’single bar" norm denotes the matrix norm comcatible witn the euclidean 
norm. Finally, the control and observation spillover coe _ "icient3 5 arc T 
are given by 


/ M ,\ 

3 . \'*^\\ <( Z I (t•) w ) 

* * c \i*l k-N*1 * / 


? \l/2 


N 


'£ 13' ; 


ORWENAL PAOK 
OP POOR QUALTs 











V 


22 


r 


t 


r 


f 


r 


r 


» 


» 




and 

. P * 2 \ 1/2 

r » ! U N C R | I < 1 K n |/ t Z (c,,40* ) (6.20) 

N K E * \j-l k«N+l J * / 

It may seem like quite a complicated business to obtain all the data 
to check (6.15); however, most of the information can be obtained from the 
finite dimensional ROM and its corresponding controller, and the rest comes 
from the structure properties a,b. The spillover coefficients (6.19)-(6.20) 
depend on the controller gains G^, and the actuator-sensor influence on 
the residual modes. Consequently, (6.15) gives an indication of how much 
spillover can be tolerated in the closed-loop system. 

Additionally, (6.15) gives the control system designer an aid to 
redesign the system for better overall stability; Quite simply,either the 
controller gains must be reduced or the actuator-sensor influence functions 
b., Cj must be chosen to minimize their influence on the residual modes. 

Of course, this latter could be accomplished if one had completely free 
choice of b., c„.; one would choose them to be in H N which is orthogonal 
to and both 7 and 3 would be zero. Please note, for a flexible structure 
this would mean the influence ^unctions were linear ccmoinations cf tne con¬ 
trolled modes only. Yet, in practice, most actuators and sensors nave a 
very localized influence (they are nearly oplnt devices), and this practical 
constraint would make it very unlikely that zero soil lover could be ac.oieved. 
Nonetheless, adjustments in the location o* actuators anc sensors could be 
made in order to reduce the spillover. Further, one could reduce spillover by 
controlling .mere modes at the cost of increasing t.oe controller ccmciex-ty 
and correspono'ing computational burden on one :n-‘ - re computa'*. jene” 
approaches to spillover reduction ana compensation (e.g. prefiltering) are 
surveyed in [12]. The apriori bound (5.151 can be used to assess the trade¬ 
offs in these approacnes for “lexible structure control. 
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3.0 CONCLUSIONS 

In order to better understand the behavior of highly oscillatory dis¬ 
tributed parameter phenomena such as mechanically flexible structures, the 
class of dissipative hyperbolic distributed parameter systems (OPS) has been 
Introduced. The fundamental property of such OPS Is that the open-loop 
system (2.1) must be exponentially stable, i.e., satisfy (2.2) for e > 0. 
Conditions have been presented in Section 3 for exponentially stable opera¬ 
tion of the infinite dimensional controller, (3.2) and (3.4) or (3.6), in 
closed-loop with the OPS; however, the principal concern of this paper has 
been with finite dimensional controllers for the DPS. Such controllers are 
practical in the sense that they can be imolemented with a small number of 
actuators and sensors and one,or more,small on-line computers. 

The finite dimensional controllers (5.6) developed here in Section 
5 are based on reduced-order models (ROMs), i.e., finite dimensional approxi¬ 
mations, of the actual CPS. The class of ROMs used here has been derived 
in Section 4 from reducing sufcspaces; a large subclass of reducing subspaces 
contains the modal subspaces which have been used extensively in flexible 
structure modeling and control. The controller synthesis is dene entirely 
in terms of the ROM. 

Since the finite dimensional controllers developed in Section 5 
must operate in closed-loop with the actual DPS and not just the ROM, it is 
crucial to analyze the spillover , i.e., controller interaction with the 
residual (unmcdelled) part of the DPS. One cannot expect to be able to :o 
this in utmost detail since the resicuals of the DPS are the part that is 
least <ncwn; however, spillover bounds can be prescribed for stable operation 
of the closed-loop system and this has been done in Section 5. The principal 
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results are Theorems 6.2 - 6.3 ; the value of these results Is that the 
stability conditions may be checked using only knowledge of the ROM and 
bounds, 3 and r, on the spillover present in the open-loop DPS. The bounds 
3 and r may be calculated directly, or estimated where necessary, from the 
ROM projections of the actuator and sensor influence functions b^ and c y 
When the spillover bounds are small enough to satisfy the stability conditions, 
the finite dimensional controller will function to enhance the behavior of 
the ROM subsystem of the DPS without causing any instability in the residual 
subsystem. When the stability conditions are not satisfied, the system may 
be stable or it may not; however, this should be a warning that some spill¬ 
over comoensation be added and the conditions rechecked. 

Several closely related DPS topics seem to require further in¬ 
vestigation: 

(1) sharper spillover bounds 

(2) bounds obtained from singular perturbations 

(3) nonreducing subspaces and model error 
(4.) non!inear control 

(5) adaptive control 

Since the spillover bounds have been obtained from regular perturba¬ 
tions, they may be sharpened by making tighter estimates, e.g., for M. in 
(6.4); connections need to be made between our results and other stability 
results for large scale and DPS obtained by Siljak [31] and Michel and 
Miller [32]. In addition, bounds obtained *or DPS via singular perturbations 
look premising for many applications, e.g., [33]. 

Many ROMs are projections onto nonreducing sucspacss; for example, 
when approximate modes are substituted for exact modes, the approximate modai 
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subspaces are not reducing subspaces. In the case of nonreducing subspaces, 
the controller interaction with the residual subsystem can occur through non¬ 
zero model error terms A^ ■ P N A P^ and A^ * P R A P N , as well as the spill¬ 
over terms. The bounds obtained here for spillover should be extended to 
model error bounds, as well. 

Finally, preliminary results In nonlinear and bilinear control of 
OPS have been obtained in [34]-[36] and the surface has barely been scratched 
in adaptive (or self-tuning) control of OPS, e.g., the survey [37]. Much 
more remains to be done in both of these areas. 
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Since (A, 3) and (A , C ) are exponentially stablllzable, there 
exist bounded G and H such that A * B G and A* ♦ C* H generate U^(t) and 
V(t), respectively: U^t) satisfies (3.5) and V(t) satisfies 

-c,t 

l|v(t)|| < Mj e 2 , t > 0 (1.1) 

where Mg >. 1 and > 0. Choose K ■ -H and from [21], p. 251, A - ,< C gene¬ 
rates Ujtt) '"here Uj(t) ■ V*(t); consequently, J * Ug^t) i | ■ MV (t)j| ■ ||V(t)j 
Therefore, 

il«(t)ll < e** 2 ' ||v Q || (1.2) 

Also, from (3.3)-(3.4) and (2.1), 


• (A ♦ 3G) v(t) ♦ BG e(t) 


or, from (3.5) and (1.2), 


I lv(t)! ( 


-C,t Z At 

< M, e 1 (la / e AT dT ; IBGI !)!>-; 


1 - €““ 


wnere A » s, - . Let = M, and : 5 min( , £,). 

I - c i c * 

if i > o, ; *v(?) ! < x ff * mit i :v (i * mbg;; ■ 1 * ■ -- ?. 

X *> 

v(:)J < m, » mH 11 v, ; ! (i * 1 -bg: 5 ^ V -—) . 

In either case, for t > 0, 

1!v(t)H < e"" 1 ! |v Q n (1 * MBSJI/iai) 

From (1.2), for t >_ 0, 

■:«(:)" < m, e“ :: "v 

*■’ m W 

"here’ore, the closed-loco system *s exocrerval'y stacle. 


f A < G, 


r 
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APPENOIX II: PROOF OF THEOREM 4,1 


By the definition of reducing subspaces (a.2c and d), A NR s P N A P R ■ 0 
and A^ i P R A P N ■ 0 and (2.1) becomes (4.6)-(4.8). Since the definition 
of reducing subspaces Is equivalent to P N and P R conruitlng with A (see 4.3b), 
then, by [25], p. 173, they also comnute with the resolvent operator R(X, A) 
of A for any X In the resolvent set; because A generates a semigroup, the 
resolvent Is non-empty. But R(X, A) is the (unique) Laplace transform of the 
transform of the semigroup U(t), e.g. [25], p. 482; therefore, the projections 
must commute with U(t). 

Let U R (t) 2 P R U(t) P R then U R (t) Is a C Q semigroup on H R because 
U R (t) U R (C) - p R u(t) P R U(t') P R - P R U(t) U(f) P R ' Pp ij(t * f) P R » 

U R (t + t‘) for t, t' >0; this uses the fact that P R commutes with U(t). Of 
course, U R (0) ■ P R which is the Identify on H R and the strong continuity follows 
from that of U(t). It is also clear that A 9 generates U 9 (t) because, for v in 

It f tw - * * 


W'" * / * • » ’ w ■ - - — - 9 - ^ ' ' 

U R (t) v ■ v _ 11m p 'J(t) v -.V « p H* vjJi . 

* * « t t—0 r R t R t-0 t 


*R AV 


> R v ■ A R v by the continuity of D R . The growth orooerty (1.3) *cl'cw$ 


from (2.2) and ;;U {t)) 
U(t). 


? R i ! ,p q;' 


because P R commutes wi t.n 
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APPENDIX III: PROOF OF THEOREM 5.1 


Since (A^, 8^) is controllable and finite dimensional, can be 
chosen so that a N > e. From (5.2)-(5.3) and (5.3), 

11v N (t)l1 < M n , h 1 |v n ( 0 )|| (III.1) 

ll* R (t)|l < M r e’ e£ (||V 0 i| * / e eT ||B r G n || ||v n (t) || dr) (III.2) 

because U R (t) v R (0) = P R U(t) v Q from Theorem 4.1. Combining III.1 and 
III.2, obtain 

i |v R (t)| | < M R r et (1 + .<yi - e'^) a N j |P N | I) (|v 0 !! (III.3) 
Therefore, using (111.1) and (III.3), 


11v(t)j| < !!v N (t)|| + 1)v R (t)|I < 


-St 


[m r + M N (e ^ + M R CI - ») i!P N I! 1 !!v 0 !! (in.4) 

Note that (III.4) is the same as Triggiani's bound in [17], Theorem 6.1. 

-A t 

From (II1.4), obtain (5.4) for nonorthogonal projections because e J - 


M ( 1 „ a 

R "N ' e 


< e 




+ Mr, < 1 + M- x. and * Ml IP- 1 
RN— K N R O R 


> 1 from 


Theorem 4 . 1 . However, when the projections are orthogonal, then ) ! v R (t)^ < 

e“ £ '' (;: v R (0) |! ♦ a N | |v N (0)i;) and :!v(-);' 2 * 1 ! v Rj (t) ' ; 2 * "v R (t): 


- (M 0 m n e " £t)2 C( I iv R (o) |! + a N ! Iv tN ( 0 );! ) 2 * i!v, N (o) ] i 2 ] 


2 . 


, .2 


< (i M N M N 2 ) (,m o m n e“ st ) 2 [ i ! v N (o)11 2 v |!v R (o) j ! 2 ] 


. 12-i 


iv( 


< (1 * ^ + x^ 2 ) ‘^ 2 M q M,j e’’* i ;V 0 


2 o o o o 

because (a * a b) * a" <_ (1 + a + z~) (a * b") which isLemma 3.1 in [19]. 










* 


32 


APPENDIX IV: PROOF OF THEOREM 6.1 


With C R - 0 in (6.1) 

*111 
q(t) - e 11 


From (5.5) and (5.10), it is easy to obtain 

lle N (t)|| ||T n v n {0)|| 


and 


“Cut t AuT 

*»(t)|| " [| l» M (0) 1! ♦ / e H 

0 


N U,T N' 

II®n ^] 2 1' Il e N^*)M ^ 

but, using (IV.2), (IV.3) becomes 

-V |!B„ h, 2 ! j ! |T n ! ! 

iiv M (t)|| < Mj| M P e N [1 * — N —^ -—31 |v N (0)| | 


because F can be chosen so that A c > 0. From (IV. 2) and (IV.4), the 
result follows because M R and are not less than unity and A. > 0. 


(IV. 1) 

(IV.2) 

(IV.3) 

(IV. 4) 
desired 
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APPENOIX V: PROOF QF THEOREM 6,2 

Let C R * 0 In (6.1) and note that * 0 and ^22 * A R' frrom 

(6.1), obtain 

r 

q(t) » q(t) 

< (v.i) 

av R (t) 

+a r v n (t) 

where !iq(t)jj satisfies (6.3) in Theorem 6.1 by appropriate choices of F 
and Gjj. From Theorem 4.1, A R generates the semigroup U R (t) with growth 
property (4.9). Using this in (V.I) produces 

j! v R (t)!' < M r e " £t (i'v 0 !! * / e £T MA^n !!q(t)M dt ) (V.2) 

where j |X 21 || » i|B R & N || 

Since (6.3) is similar to (111.1) and (V.2) is similar to (III.2), 
the proof of this theorem - Theorem 6.2 - can be carried out in the same 
manner as the orcof of Theorem 5.1 in Appendix III (where q(t) ta<es the 
place of v^(t) and replaces M^.) 
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APPENDIX VI: PROOF OF THEOREM 6.3 

Let w(t) - Cq T (t) v R T (t)] T and, from (6.1), obtain 


w(t) * J w(t) + *(t) 

— 0 ^ 


(VI.1) 


where A Q ■ 


*11 0 
*21 A R 


and AA q * 


u m 12 
0 B R H 11 C R 


From Theorem 6.2, A. generates a semigroup S Q (t) with growth property, from 


(6.5): 


IS (t)il < M e e 


-€t 


where = M^M*. Also, from (6.1), 




1/2 


(VI.2) 


(VI.3 ) 


Since $ Q is bounded, J Q + ±J Q generates a semigroup S(t) (e.g., [13] p. 30 
or [6] p. 210); therefore 

w(t) * S(t) w(o) 

However, *rcm (VI.:), 


(VI.A) 


w(t) * S a (t) w(o) + / S Q (t-t) d ” 

"Taking bounds in (VI.5) and using (VI.3), obtain 

t .. 

"w(t) M < M e e’ st [| !w(o)!! - :.* e"” 1 : *(") ! ! d ” J 
” E 0 

Let n(t) = e £t ''w(t) ana, -rcm (71.5), obtain 
nit) i n(o) * 7.' h(-) dt 


; 71.5; 


(VI.6) 




From the Sronwall Ineouality (e.g., [6] p. 124 or [14] p. 111-45) apolied to 
(VI.7), 
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h(t) <. Mg e‘* h(o) (VI.8) 

or, equivalently, 

I|w(t)|! < e‘ £/t ||w(o)|| (VI.9) 

where e' is given by (6.6) and this is the desired result (6.7). Clearly, 
the semigroup S(t) is exponentially stable when t is positive. 
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ABSTRACT 


Large flexible aerospace structures, such as the solar power satel¬ 
lite, are distributed parameter systems with very complex continuum des¬ 
criptions. This paper investigates the use of finite element methods’ 
to produce reduced-order models and finite dimensional feedback con¬ 
trollers for these structures. Our main results (Section 5.0} give con¬ 
ditions under which stable control of the finite element model will 
produce stable control of the actual structure. 
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1.0 Introduction 

Large flexible aerospace structures, such as the solar power satel¬ 
lite^, are distributed parameter systems with very complex continuum 
descriptions. Often, their dynamics are modeled by large-scale 
systems of ordinary differential equations which are obtained via the 
finite element method for structures. To obtain implementable, i.e. 
finite dimensional, feedback controllers for these structures, we shall 
use this model as the basis for the control synthesis. 

The question of closed-loop operation of such controllers within 
the actual structure is certainly a fundamental one. In'this paper, 
we shall point out some connections between the finite element models 
and the closed-loop stability. These results will be used to give an 
answer to the following version of the above fundamental question: 

When does stable control of the finite element model ensure stable con¬ 
trol of the actual structure? 

"or distributed parameter systems, in general, the finite element 
method (FEM) for time-dependent problems is very well explained in [2] 
ChaDter 7, [3] Chapter 9, and the excellent survey [4]. The related 
approaches called "spectral methods" are described in [5]. The use of 
FEM for analysis of distributed parameter control problems involving 
interior and boundary controls has been considered, for example, in [6] 
[10]; all of these references are concerned with parabolic (i.e. diffu¬ 
sion-like) distributed systems and, with the exception of [10], they 
are concerned with open-loop (i.e. non-feedback) control. 

In contrast, our concern is feedback control of hyperbolic (i.e. 
highly oscillatory) distributed parameter systems and the use of FEM 
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to develop low-order finite dimensional controllers for such systems. 
Control approaches for structures have been surveyed in [11]-[13]. The 
most popular approach for structural control is the resolution of the 
structural vibration into normal modes, (semi) discretization of the 
distributed system by truncation of this modal expansion, and controller 
design based on the truncated modal model; this is called modal control . 
Modal analysis of structures has been described in many references but 
an enjoyable survey is [14]. The use of modal analysis in feedback 
controller design depends on a complete knowledge of the modal data 
which is rarely available for something as complex as a Targe aerospace 
structure. Instead, this modal data is approximated by finite element 
structural analysis computer codes (e.a. NASTRAN), and it is this approx' 
imate mooal model that is used for the controller design, e.g. [15]. 

The stability of closed-loop operation of the finite dimensional 
modal controller in the infinite dimensional distributed parameter 
structure is questionable even when the actual modes are known, as shown 
by a simple example in [16]. The mechanism for instability is the 
combination of control and observation spillover in the residual (un¬ 
modeled) modes of the structure [13]. However, closed-loop stability 
bounds have been obtained for highly oscillatory distributed parameter 
systems with low-order finite dimensional controllers [17]. These 
bounds indicate the tolerable levels of spillover for stable operation, 
but they are based on "reducing" subspaces. For large aerospace struc¬ 
tures, this corresponds to exact knowledge of the modal data used in 
the controller, although the residual modal data need not be so pre¬ 
cisely known. This paper addresses the closed-loop stability question 
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for structure controllers based on finite element approximation of the 
modal data. The distributed parameter control problem for flexible 
structures Is described in Section 2 and the finite element approximation 
in Section 3. The corresponding finite element modal controllers are 
presented In Section 4, and our main results on the closed-loop stability 
are Theorems 1 and 2 in Section 5. 


2.0 The Distributed Parameter Control Problem for Flexible Structures 

We consider the class of flexible structures that can be described 
by a generalized wave equation; this represents an idealization of many 
mechanically flexible structures. The generalized wave equation is 
given by: 


S 2 u 

17 


♦ bljr + A u ■ F 
9t 0 


( 2 . 1 ) 


which relates the(vector of)displacements u(x,t) of a body n. a bounded 
open set with smooth boundary 9fi in n-dimensional Euclidean space R n , 
to the applied control forces F(x,t). The operator A Q is a time-invariant, 
self-adjoint differential operator with compact resolvent and lower semi- 
bounded spectrum. The domain D(A ) of A Q is dense in the Hilbert space 
L 2 (;l) with (.,.) denoting the usual inner product and li-!l 0 denoting 
the associated norm. The natural damping in the structure is modeled by 
the term b It where b could be any bounded self-adjoint operator commuting 

0 L 

with A_ but is taken here to be a non-negative constant, for convenience, 
o 

The control forces 


F(x.t) » B f(t) * I b.(x)f.(t) 

J i«l 

are provided by M actuators with influence functions b^(x). The dis¬ 
placements are measured by P averaging sensors 


( 2 . 2 ) 
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y(t) - C 0 u(x,t) (2.3) 

where y^(t) ■ Cj(x)u(x,t)dx with j ■ 1, 2, . . , , P. The actuator 
and sensor functions b^(x), Cj(x) are in L (f<) and normalized to have unit 
integral. When the support of b^(x) is in a small neighborhood of a point 
x^, we say it is a point actuator and, similarly, we define a point 
sensor . The point actuator and point sensor situation is our principal 
interest here. This class of distributed parameter systems includes 
interior and boundary control of vibrating strings, membranes, thin beams 
and thin plates. Although, only displacement sensors are considered in 
(2.3), this is not a restriciton - more general sensors may be modeled. 

It is well known ([18], p. 277) that the spectrum of A Q contains 
only isolated eigenvalues X k with corresponding eigenfunctions o k such 
that 


Xg ^ • 

and A Q c k * X^. We will assume that X^ is positive. The resonant 
mode frequencies m. of the structure are given by * (X^) IL and the 
corresponding eigenfunctions are the mode shapes . Thus A q satisfies 

(A o u,u) 0 > a!!u!, a > 0 (2.4) 

1/2 2 

and has a square root A . Every vector ucL [?.) has a unique represen¬ 
tation 


u(x) « : 'J.a(.x) 

k-1 * k 


(2.5) 


N * R 


where u^ ■ (u.c^Jq and we define the orthogonal projections P 0 , P Q by 


P ' u - : u. 


k-l 


( 2 . 6 ) 


P o Ru 


Vk 


k-N+1 
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Let V be the domain of A Q and W be the domain of A Q 
A is defined in H by 


1/2 


A new operator 


a) 0(A) - 

V x W S H 1 


f w 1 

b) 

1-bw - A q uJ for 

enerov inner product (...) is 

( h] • 

r t 

u 2 

w ) 5 <A fl 1/2 u. 

J 

1*2' o 

u-|. u 2 e V 

and w.j, Wg c W, < 


2'o 


V 2'o 


(2.7) 


( 2 . 8 ) 


space , is defined as the closure or in this energy inner product. The 
associated energy no rm is denoted by |M| and is a measure of the total 
potential and kinetic energy in (u,|^) where u is a solution of (2.1). 


T T 


Let v ■ [u^ , I*: ] be in H and write (2.1), (2.2) and (2.3) as 


ff * Av * Bf , v 0 eH 1 


(2.9) 


Cv 


where B ■ 


* 

0 

w 


and C » [C 0 0]. 


It is easy to see that 



2 2*2 

and |jv|| » l [X k u k + u k ]. Also, A generates a C Q semigroup U(t) 

hy [19], p. 11-59. The following result, due to Goldstein and Rosencrans 
[20], gives conditions under which this semigroup U(t) is (slightly) 
exponentially stable and hence (2.9) is a (strictly) dissipative hyperbolic 
system: 
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if b < 2a then 


v(t)n < m y et i |v(o) 


( 2 . 11 ) 


t > 0 


where M Q • (y) 1 ^, e ■ b/2, and y ■ [1*(b/2a)][1-(b/2a)]‘^ 


Note that 

the hypothesis for this result means that e < a, l.e. the damping must 
be smaller than the square of the lowest mode frequency . Another bound 
of the type (2.11) is available In [21] without the restriction t < a; 
however, the damping coefficient b is not quite so directly related to 
c in (2.11). The bound used above will be very convenient to illustrate 
our results. 

Recalling (2.6), we define the projections P N , P R on H by 


- - 


N 

u 


P u 


. 

0 

. W , 


P 0 N w 


r 2) 


r _ 


r B r 1 

U 


V 

w 

" J 


P P w 
m 0 J 


and note that they are orthogonal in the energy inner product and that 
H n - P n H and H R »P R H form reducing subspaces c* H, as described in [17]. 
In fact, H n is a modal subspace, and v N e P N v and v R » P R v satisfy: 


3v N 

1C 9 V N + B N f 
3v„ 


3t 


A P V R + B R f 


(2.13) 


y • C N v N + C R v R 
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where ■ P N A P N , B N ■ P N B, C N ■ C P N , and, correspondingly, for 
(A r , S r , C r ). Note that A^, B N , C H may be identified with the matrices 

. [C N ° 0], respectively, with A N • diag [X^ 

. . . , X N ], B n ° ■ [vbj,d k )] T , and C N ° * [(Cj,e k )], and, if point sen¬ 
sors and actuators are used, then B,^° ar.d C N ° involve only mode shapes 
$ k evaluated at the actuator and sensor locations. The terms S R f and 
C R v R are called, respectively, control and observation spin over , and 
it is through these terms that the residual modes v R are affected 
by any feedback controller which accepts the sensor outputs y and 
generates control commands f. 

Modal controller? can be based on the reduced-order model (ROM) 
whose states are v N and whose parameters are (A^, B N , C N ); such a model 
is obtained by ignoring (for the purpose of controller synthesis) the 
residual modes, i.e. setting 3 R * 0 end C R • 0. The controllability 
and observability of the ROM can be determined from U N , c N °, C N °) when 
the damping b is small ([22], Theorem 2.1). Of course, it is quite 
easy tc assess the controllability and observability of this latter 
system,and the conditions depend on the interaction of the actuator¬ 
sensor influence functions with the mode shapes e-j, . . . , $,j. 

If the ROM is controllable and observable, the usual f inite dime n¬ 
sional , state space controller can be designed based on the moaal ROM 
(e.g. [23] or [24]). It has the form: 


[■ 


0 


-bl 


N 


( 


f'(t) * v N (t) 

v N (t) * A n v N (t) * B n f(t) ♦ * N (y(t) - y(t)) 


(. y(t) ■ C N v N (t) , v N (0) * 0 


( 2 . 12 ) 


(2.14) 
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This controller has dimension 2N (twice the number of modes in the ROM). 
Since the ROM is controllable and observable, 6.. and K„ can be chosen to 

iN N 

achieve any desired stability margins o N and c p for 6^ and 

A N - Kn C N respectively; optimal control, as well as pole placement 
techniques, may be used successfully for the controller design. 

Thus, (2.13)-(2.14) represents an implementable (finite dimensional) 
modal controller for the flexible structure (2.1)-(2.3), or equivalently 
(2.9). It is based on a modal ROM, i.e. it is designed to control the 
first N modes of the structure (of course, different sets of N modes 
could have been used - it is not necessary that they be the first N 
modes). The basic assumption made in synthesizing such a modal con¬ 
troller is that the N modes chosen for control will adequately represent 
the important vibrational behavior of the whole flexible structure ; this 
assumption can often be satisfied in practice, and it is the basis for 
most structure control!ars. 

Having synthesized an implementable modal controller to adequately 
tailor the vibrational response of N critical modes in the flexible 
structure, we must not think w? are done. The closed-loop behavior of 
our controller is affected by the residual modes in Hp for which we have 
not designed, as well as the controlled modes in for which we have 
(. ’signed. These interactions occur throuch the control and observation 
spillover terms defined in (.2.13) and can drastically alter the desired 
performance of the controlled structure; in fact, an example is presented 
in [16] where the closed-loop system becomes unstable. Consequently, 
the closed-loop stability analysis presented in [17] is extremely impor¬ 
tant; these results give an indication of the amount of spillover that 
can be tolerated in the closed-loop system without loss of stability. 


- “ill) II ' III! III! ill ■ 






3.0 The Finite Element Approximation and Open-Loop Convergence 


Unfortunately, all of the controller synthesis and closed-loop 
analysis presented so far is based on reducing subspaces and and, 
in particular, on exact knowledge of the system modes ^, . . . , <|> N 
used in the ROM; this modal data must be approximated in most cases 
due to the complexity of the actual partial differential equations and 
associated boundary conditions which describe the structure (2.1). Such 
approximations do not lead to reducing subspaces, in general; therefore, 
further discussion is necessary. 

The most natural approach (and the one we shall use)' seems to be 
to approximate the modal data required for the modal ROM (A^ t , B^, C N ) 
via the finite element method (FEM) for time-dependent distributed para¬ 
meter systems (2.9); the data required will be the ROM mode frequencies 
and the corresponding mode shapes 4> k (x) evaluated in neighborhoods 
of the actuator and sensor locations when point devices are used. Once 
this FEM approximated modal data is available, the controllability and 
observability of the ROM may be tested,and the FEM based modal controller 
may be synthesized using the approximate modal data in (2.13)-(2.14) as 
though it were exact. Although the FEM approach requires a great deal of 
"number-crunching" to obtain the approximate modal data, this data is 
needed for the open-loop structural analysis and is calculated anyway, 
regardless of control considerations. Also, the large scale computer 
codes necessary for such FEM calculations are available and in general 
use (e.g. NASTRAN and others). 

This makes the FEM approach to structure controller synthesis 
appealing; however, the closed-loop stability analysis is much more diffi 
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cult. The fundamental question to be answered is - when does stable 
control of the FEM model produce stable control of the actual structure? 

To begin to answer this question, we must give a more precise descrip¬ 
tion of the FEM approach. 

A precise description of the FEM approximation of time-dependent 
distributed parameter systems like (2.9) is given, for example, in [2] 
Chapter 7, [3] Chapter 9, or the survey article [4]. We point out 
that, historically, the FEM was developed specifically to numerically 
approximate the behavior of a given distributed parameter system or set 
of partial differentia! equations; in control terms, it was intended 
to give close approximation of the dynamic behavior of open-loop sys- 
tems like (2.9). However, we will use it in a very different way, 
namely to develop reduced-order models of the system (2.9) for synthesis 
of finite dimensional, feedback controllers to be used in closed-loop 
with the system (2.9). In the usual aoplications, the stability of the 
ROM plays a big role, and hyperbolic systems like (2.9) are notorious for 
their lack, or very low margin, of stability when the FEM is used. In 
our application, the stability of the ROM is less important because the 
feedback controller will be designed to improve its stability; more 
important to our application is the stability of the residual (unmodeled} 
pa-t of the system over which we have no direct control. 

We describe the FEM approach to finite-dimensional controller syn¬ 
thesis for the distributed parameter system (2.9). Let be an increasing 
sequence of finite dimensional subspaces of the state space H for (2.9). 
Each subspace has dimension N. To make life easier, we assume that 
each is a subspace of D(A) so that its elements satisfy 


ORIGINAL PACK L 
OF POOR QUALIT* 
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the boundary conditions for A; however, so called non-conforming elements 
may be used in the more general case. In the FEM, each subspace H N 
consists of splines (i.e. piecewise-polynomial functions) of fixed degree 
defined over a mesh (usually, of triangles) laid out to approximately 
cover the structure il (see [3] Chapter 6). No matter how irregular the 
shape of the boundary of f2 such meshes can be fitted very closely; this 
is one of the principal assets of the FEM. To each mesh, a normalized 
mesh parameter h (where 0<h<_l) is assigned so that the mesh is refined 
as h -*• 0 and the dimension N of the subspaces increases indefinitely. 

3.1 The Galerkin Approximation 

Let P^ be the orthogonal projection from H into H^; this is called 
the Galerkin projection . The corresponding orthogonal projection onto 
hJ is called P R (i.e. P R ■ I - P N ) . The "rate of convergence" of H N 
to H is said to be of order q when 

j|P R v|| < K h q (3.1) 

for v in H (or D(A)); this rate is related to the ability of splines in 
H n to interpolate functions in H. We shall not be concerned with the 
rate of convergence ; consequently we write (3.1) as 

lim j|P D v!| * 0 for v in H (3.2) 

N-k*> f 

and suppress the dependence on h for our discussion. 

Let -^(x), . . . , ij> N (x) form a basis in H N (i.e. they are linearly 
independent). These functions are called patch functions or assumed mode 
shapes . An approximation of the solution v(x,t) of (2.9) can be formed 
in H N by 

N 

v N (x,t) « 1 v k (t) i|» k (x) 


(3.3) 
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i.e. assume separation of time and space variables with all spatial 
variation lumped into the patch functions ^ k (x). The choice of the 
coefficients v^(t) remains; these are obtained by substitution of (3.3) 
into (2.9): 

3v m 

^ * A v N + Bf + E (3.4) 

where E is the equation error , and the v^(t)'s are chosen so that 


p n U) - o 


(3.5) 


This is called the Galerkin approximation ; when it is carried out with 
the spline subspaces described above, it produces the FEM approxima ¬ 
tion (3.3) where the coefficients v^(t) are given by the entries of the 
solution vector v^(t) = [v^(t), . . . , v^(t)]' for the system of 
ordinary differential equations: 


m n ^ = a n ^ + b n f 


(3.6) 


where M N = [(i|» £ , ^ k )j , 

An = A » 

and § N f * [U £ , Bf)]. 

The matrix is symmetric and positive definite because H' k ( x ))k=l are 
linearly independent: 


- ~ N N N N 

z‘ M,, 2 * I I Z 5 (V^) Z I, = ( 1 1 and 

“ N ~ £*1 k=l 


■Jr £’ v k' k 
N 


£=1 11 k=l K K 


0 = 2 / M n z if and only if I z k ^ k = 0 which leads to the z k « 0 
due to the linear independence of the u> k « Therefore, (3.6) can be solved 
uniquely forj^it) whenever jv fJ (0) is specified, and hence the FEM approxi¬ 
mation (3.3) is obtained. It is assumed that v N (0) is given by the 


vector of coefficients of 
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v n(°) 


P N v 0 


(3.7) 


expanded In the basis x)}. Note that 


V N * P N V N 


(3.8) 


The approximation (3.3) is called a semi discretization of (2.9) because 
time t remains continuous in (3.6). 

It should be noted that to obtain the most analytical benefit from 
the FEM, the approximation (3.3) should be obtained from the “weak" form 
of (2.9); however, we omit discussion of this technicality and refer the 
interested reader to [3]. 

3.2 Convergence of the Approximation 

Let e(x,t) = v(x,t)-v N (x,t), the FEM solution error , and consider 
from (2.9), (3.4), (3.5) and (3.8): 

3v, 


—- s A v + B f 
3t N V N d N t 


ae 


v N (0) 


P N v 0 


at s a rn v n + Ae + S R f ’ * P R v 0 


C N V N + Ce 


(3.9a) 

(3.9b) 

(3.9c) 


where A N * F N A P N , * P N B, C N « CP N , B R = P R B and A RN = P R A P fr 
Since P N is orthogonal projection onto the finite dimensional subspace 
H n , (A n , B n , C n ) may be identified with appropriate matrices easily re¬ 
lated to the ones in (3.6) by choosing the orthonormalized basis for 
N 

obtained from ; i n addition, A R ^ is a bounded , linear operator. 

In general, A RN is not zero, and hence and H R are not reducing sub¬ 
spaces for A; furthermore, v^ * v, and so and e are not orthogonal 
to each other. Consequently, except for th>» special case where the exact 


modal data is known and H R is chosen to be a modal subspace, the FEM 
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does not partition (2.9) into the previous form (2.13) analyzed in [17]. 
The FEM partition (3.9) must be analyzed separately. 


From (3.9b) and the fact that A generates the Cq semigroup U(t) 


obtain: 


e(t) « U(t) P R v Q + U(t-i)[A RN v n (t) + B R f(t)] dr (3.10) 

Take norms in (3.10) and use (2.11) to obtain: 

||e(t)|| < M 0 e* ct (.||P R v 0 || ♦ e £T [|!A RN v(t)j| 

+ l|B R || ||f<r)|| ♦ 1!A rn || !!e( T )|j] dr! (3.11) 
Now use Gronwall’s Inequality ([25], p. 124) on (3.11) to obtain: 


!|e(t)]i < Mq e N (||P R v 0 M + £ e £T [|iA RN v(t)|! 

+ !!B R i! !!f(T)j i] dt) (3.12) 

where e N = e - Mq!IA RN i!- This can be used to prove the following 
convergence result: 

Theorem 3.1 : The following two conditions guarantee convergence, over any 
finite interval of time, of the FEM approximation (3.3) to the actual 
open-loop solution of (2.9): 

(a) 1 jA rn |I < a where a is independent of N 

(b) The approximation is consistent , i.e. lim !!A RN v|j « 0 for 
all v in D(A). 

Furthermore, conditions (a) and (b) are equivalent to 

11m M A rn !! s 0 (3.1: 

N-*-® 

Proof: From [18], p. 151 lemma 3.5, we obtain that (a) and (b) imply 
(3.13) because D(A) is dense in H. Conversely, it is easy to 
see that (3.13) imolies (a) and (b); so 


(3.13) 
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these conditions are equivalent. Use (3.12) with e N replaced by e-MqO 

enc let 0 < t < T. From (3.2) and (3.13), obtain 1im |je(t)l| * 0. # 

N*® 


This is only a slight modification of the usual open-loop convergence 
results (e.g. [5] Sec. 4), but it will be convenient for our needs. 

A Priori estimates of the convergence rate of e(t) as a function of the 
resh parameter h can often be made. 


3.3 An Alternate Form of The FEM Approximation for Flexible Structures 
We want to carry the development of the FEM approximation (3.3) and 
(3.6) a bit further. We can rewrite (3.3) as 


v..(x,t.) « [u N (x,t), w N (x,t)] T 


wnere 



N 




r u,(x,t) * z v. (t) 

1 ' k*'l k 

(3.14a) 

<j 

V*-*) = 2 v k (t) 

k a i 

t k 2 (x) 

(3.14b) 

wnere v k = O k 1 ^ 2 ] T is 

in H n . Now the entries 

in the matrices M fJ , 


A,j, in (3.6) are given by: 


■ (Kj 1 . *0 *k’) 0 * ( *H 2 ' (3 - 15i ) 

it..A* k ) . Aq t k ? ) - I** 2 , A 0 - b(*/, » k 2 ) (3.15b) 

U 0 0 

M ? 

(c.,Bf) * l (0b.) f. (3.15c) 

* i*l £ 1 0 1 


From "his data (3.6) may be solved for the coefficients v^(t) of the FEM 
approximation (3.3). 

An alternate form of the FEM approximation is obtained by choosing 

0 0 
spline subspaces $ in H 0 with 6,(x), . . . , e N (x) as a basis for S . 
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Instead of making the approximation in the state space formulation (2.9), 
we do it directly on (2.1): 


u N (x,t) * I u.(t) 6.(x) (3.16) 

n k«l K K 

where u k (t) satisfy the system of ordinary differential equations: 

^ u N + b u N + A n ° u N • B n f (3.17) 

where u N = [u 1 , . . . , u N ] T , M N * C(S £ , 6 k ) ], A N ° * [(© £ , Aq 6 k ) ]. 
and 3 n ° * [(e £ , b^)]. It is easy to see that M N and A^° are positive defi¬ 
nite matrices because Aq is a self-adjoint, positive differential operator. 
From the extended eigendata problem: 


A N ik * X k ik 


(3.18) 


we obtain the approximate mode frequencies 


“k * ( V 


(3.19a) 


and approximate mode shapes: 


: k (x) ■ j, e k (x) 


(3.19b) 


where z. » [z, . . . , z kN ] 1 . When we form the nonsingular matrix 


V N ■ til.*1(3 


(3.20) 


we find that 


V N *N V N “ 1 


(3.21a) 


V N A N V N * A N “ dia9 [X T • • • ’ X N ] 


(3.21b) 


The modal coordinates q f , are defined by 


-N ' V N *N 


and (from (3.17) we obtain 

£n + b S.N + A N° S-n V N B N f 


(3.23) 
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This system of ordinary differential equations is decoupled into N second 
order systems and, hence, is quite easy to solve for . From (3.22), 
we can obtain the coefficients for the alternate approximation (3.16). 

The alternate formulation (3.16) of the FEM for flexible structures 
is the standard one used in structural analysis; it is related to the 
one obtained in (3.3) or, equivalently, (3.14). Appropriate versions 
of (3.2) and (3.9) and a convergence result similar to Theo. 3.1 can be 
produced; fcr example, see [2] Sec. 7.3. For convenience we shall con¬ 
tinue to use (3.3); however, alternate versions of our later results 
hold for (3.16). 

4.0 Feedback Controllers Associated with the FEW Approximation 


The FEM reduced-order model associated with (2.9) is defined on H N 
and given by 
3v. 


N 


3t 


a n v n + V 


V N<°> * P N v 0 


(».l) 


* m C N V N 

where (A N , B N , C N ) are defined in (3.9) and are determined by approximate 
modal data. Since H N is a finite dimensional subspace, (A N , B N , C N ) may 
be identified with their matrices in an appropriate basis of H^, and 
(4.1) is equivalent to a lumped parameter, state variable system for 
which a well developed feedback control theory exists, e.g. [73] or [2d]. 

The controllability and observability of (A N , B N , C^) ere -idsily 
checked via simplified conditions obtained in [16]. These conditions 
are the same even though approximate, rather than exact, modal data are 
used; the conditions apply because of the low level of natural structural 
damping present in these structures (i.e. b is small in (2.1)). 
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by 


The FEM Feedback Controller is based on the ROM (4.1) and defined 


f " G N v N 


^ °N + B N f + K N 


(4.2) 


U N V N 


v N (0) 


where, due to the controllability and observability of the ROM, we can 
adjust the controller gains G N and sc tn: ^ + B N G N and A^ - K N C N 

have any desired stability margin. The controller (4.2) has finite dimen¬ 
sion N (where N * dim H N ); so we can abuse notation by using v N instead 


3v 


of . This controller consists of a linear feedback control law and 
full order state estimator (full order in the sense that it is matched 
to the full order ROM). Much lower order controllers than (4.2) may 
be developed, but we leave that issue for another time. 

We define the estimator error e^ = v,, - v^ and, from (3.9c) and 
(4.2), obtain: 


®N " (A N ’ K N V e N + K N Ce 
Also, from (3.9a) and (4.2), we have 

V N * (A N + B N V V N + B N G N e N ( 4 - 4 ^ 

If there were no FEM solution error (i.e. e ■ 0), then (*.3) and (4.4) 
would be designed with a desired stability margin. Consequently, the 
controller (4.2) would stabilize the model (4.1) by design ; however, 
our principal interest is the closed-loop stability of the actual struc¬ 
ture (2.9) with the controller (4.2). 
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5.0 Closed-Loop Stability; Main Results 


The closeo-loop system consists of the actual structure (2.9) with 
the FEM controller (4.2). The closed-loop system state Is [v e N ] T where 
v • v N ♦ e and e N * v N - v N ; equivalently, this state may be written 

“ * [V N «N 


where u satisfies (from (3.9), (4.3) and (4.4)): 


« ■ A p u: 

ot C 


where A £ ■ 


‘11 


A il A 12 

L A 21 A 22 


and 


a n * b k g n b n g n 1 

A. - * 

0 

0 A N ■ *N C N 

’ 12 

V. 


(5.1) 


'21 


£ a rn + b r g n b r S3 anfl A 


22 


Sines represents a finite dimensioral system, it is bounded and gene¬ 
rates a Cq semigroup U N (t) with the following growth property: 

-e..t 

1 :U N (t) * 1 < M n e r ' , t > 0 (5.2) 

where ^ ^ 1 and > 0; c f , is the designed stability margin obtained 
in Sec. 4.0 by the choice of controller gains 6^, and K^. We have the 
following closed-loop stability result : 


Theorem 5.1 : If the FEM model (A^, B N , C N ) is controllable and observable 
for eve-y N and if conditions (a) and (b) of Theo. 3.1 are satisfied, then 
controller gains G N and can be chosen so that A^ generates a C Q semi¬ 
group U c (t) such that 
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u(t) ■ U c (t) u(0) 


(5.3) 


and 


■‘C* 


I|U c (t)|t < M c « . t > 0 

mere ||u|| i (,||v N || 2 ♦ ||« N I | 2 ♦ ll «!| 2 )' 72 . 


(5.4) 


M C ’ % M 0 0 * 


IIS C I! !IS C H 


2 1/2 


e», - e 


e,, - e 


) 


£ C • E - M c S N 


“ * 11**1.11 * 2 li B Rl 


’N " H"RN" ' hi-rii i,G N i! , and Mg, e are obtained from (2.11). 

The closed-loop system is stable, if 6^ is bounded for all N and,for N 
sufficiently large, 

Mj. B n < c (5.5) 


Proof: Choose G N and < N so that e N > c since (A N , 6^, C N ) controllable 
and observable. Let A r * Ag + iA where 


A n A i2 ”| 



1 

O 

O 

_ J 



AA « 

1 

O A ' 

_ 22 - 

| and 

1 

[ A 21 0 J 


Note that Aq generates a Cq semigroup Uq(*) such that 
! |U Q (t) 11 < M c e 0 , t > 0 

where Eg = min (e N , c) ■ e; this follows from arguments similar 
to those used to prove Theo. 5.1 in [17]. Since !I&A|| • 

11A^i)! < 6 n . AA is a bounded perturbation of the semigroup gene' 
rator Aq; hence, by [25], Theo. 10.9, p. 210, A c generates a C c 
semigroup Ug(t) which satisfies (5.A). The conditions (a) and 
(b) of Theo. 3.1 imply that 
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lit MArhII " o 

and this, together with (3.2) and the boundedness of G^, Implies 
11m * 0. If N can be chosen large enough to satisfy (5.5), then 

c c > 0 and the closed-loop system (5.3) Is exponentially stable. 

Note that &£ may Increase as N ■* * so It Is not always certain 
that (5.5) can be satisfied even though S N - 0. * 

Theorem 5.2 : In Theo. 5.1, if there exist positive constants a, and 
&2 such that 

^Vn V ' (5 - 6 > 

then there is a finite integer N* such that, for all N _> N*, (5.3) is 
exponentially stable. 

Proof: This is a corollary to Thee. 5.1. From (5.6), we outain 
M. £„ < a, “ •* 0 as N - Therefore, choose N* such 

W N I It 

that o 1 < e and, for all N > N*, (5.5) is satisfied: 

the result follows from Theo. 5.1. * 

The above results provide conditions under which ths clcsec-loop system 

(5.3) is stable and the FEM solutior error e{t' satisfies 

-c r t 

!l«(t)l! < IMt;l! < 2 m c e L |! v Q 11 , t>0 (5.7) 

because 

"w(3) ! ' 2 - 2 m P n v c '! 2 * !! P r v q M 2 < 2Mv Q 1 ! 2 . 

This is an a priori convergence estimate for the closed-loop system. 






tr 

l 

I t 
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f ? » 

IS, 


Since e ■ v - v N ■ (v - P N v) + (P N v - v N ) and * - P N v -_P R v * 0 
as N « by (3.2), we may be interested in the behavior of the term 
P N v - v N as N -*» «. However, from the properties of the Galerkin pro¬ 
jection P N , 


P N * ’ P N P R * + P N (P N v ' V * P N V ‘ W K 
P D e » P B P D v + P D (P W v - vj - P„ v 


(5.8a) 

(5.8b) 


. R c ' r R r R Y ' r R' r N ’ ' ’N' ' R 
i. e . > p^v - v N * P N e and v - P N v * P R e form the orthogonal components 
of e in H n and , resDectively. From (3.9b) and (4.2), obtain: 


3P n e 
at 


3v r 


= A N P N e + A NR V R 


( TT = A RN P N e + A R V R + (A RN + B ° G ^ V N + B » e » 


R V N T R W N C N 


(5.9a) 

(5.9b) 


where v R = P R v R , A NR = P N A P R , etc. Also, (4.3)becomes 

*N ' (A N ' ^ C n) e N + K N C M P N e + K N C R V R (5 - 10) 

Let u » [v N ,e P e , v R ] T be the closed-loop system state, and, from 
(4.4), (5.9) and (5.10), obtain 


3l0 n 

f ">C“ 


where * 


Ml m2 


L Q 21 Q 22 J 


and 


Ml 


71 


a n + b n g n b n g n 

°1 


i ° 1 

0 A n - K f , C N 

k n c J 

1 * Q 12 * | 

k n c r 

o 

o 

A n . 

1 

a nr 

• a rn + b r g n b r g n a rn^’ 

and 

Q 22 * a r 



(5.11) 


C -Ji 
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The following alternate closed-loop stability result may prove useful in 
those cases where conditions (a) and (b) of Theo. 5.1 are not satisfied 
by the FEM approximation: 

Theorem 5.3 : Assume, for every N, 

(1) the ROM (A^, B n> C n ) is controllable and observable 

(2) A^ is (exponentially), stable with margin 

(3) Ap generates an exponentially stable Cq semigroup U R (t) with 
I |U R (t) 11 <. Mp e R , t > 0 

(4) Hal I |A nr ] I - 0 

f'K ® 

Then controller gains G N and can be chosen so that Qq in (5.11) gene¬ 
rates a Cq semigroup Uq( t) with 
-c c t 

!|U c (t)I| <K c e , t > 0 (5.12) 

where Mq = max (M^, M R ) , 

C c * c 0 ' r N ’ £ 0 s min ^ c N’ C R^ • 

r N s I 1 ' C r' I + ! I a nrI ! 

and Q^-j generates the Cq semigroup i(t) with 
! |u-,-|(t) | i < e ° wt , t > o. 

If is bounded for all N and, for N sufficiently large, 

M c r N < (5.13) 

the closed-loop system (5.11) is exponentially stable. 

The proof of this result is sufficiently similar to that of Theo. 5.1 
that we shall omit it. In some cases, the conditions (l)-(4) of Theo. 
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5.3 may be easier to verify than those of Theo. 5.1. Again, the stability 
of (5.12) leads to an a priori convergence bound on the FEM solution 
error: 

" e C l 

I|e(t)|| < 2M C e c ||v Q || , t > 0 (5.14) 

6.0 Conclusions 

The use of the finite element method (FEM) for numerical simulation 
of complex distributed parameter systems is well known. In this paper 
the use of the FEM to produce reduced-order models and finite dimensional 
feedback controllers for distributed parameter systems ha-s been investi¬ 
gated with emphasis on application to control of large aerospace struc¬ 
tures. Our principal concern has been to determine the stability of 
tne controller in closed-loop with the actual distributed parameter 
system even though the controller is synthesized using a FEM approxi¬ 
mation. Our main results. Theorems 5.1-5.3, give conditions under which 
stable closed-loop operation can be expected. Such conditions are 
especially desireable for large aerospace structure control because the 
FEM model is the tool most widely accepted for structural analysis and 
large scale computer codes already exist to produce such models for very 
complex structures. 
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ABSTRACT 


Implementable feedback control of distributed parameter systems must 
often be based on reduced-order models due to the infinite dimensional 
nature of the actual open-loop system. In this paper, we analyze the be¬ 
havior of controllers designed via reduced-order models obtained with singu¬ 
lar perturbations techniques. When such controllers are used in the actual 
distributed parameter system, the closed-loop stibility is in question. The 
results presented here provide bounds on the smallness of the singular per¬ 
turbations parameter to ensure stable operation; such a priori bounds may be 
used to evaluate potential reduced-order controllers for distributed para¬ 
meter systems. 
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1. INTRODUCTION 

Distributed Parameter Systems (DPS) are described by partial differen¬ 
tial equations. Examples of DPS include heat diffusion and chemical pro¬ 
cesses, wave propagation, and mechanically flexible structures, such as 
large flexible spacecraft and satellites, high-speed aircraft and surface 
vehicles, and large civil engineering structures, e.g., bridges and tall 
buildings. 

The state spaces for these systems have infinite dimension. Thus, it 
Is impractical or impossible to Implement feedback controllers based on 
complete models of these systems; hence, reduced-order models must be used 
for the controller design. A great variety of reduced-order modeling tech¬ 
niques exist for general systems, e.g., the bibliography of [1], and new 
techniques are presently being developed in the specialized areas 11k* large 
aerospace structures, e.g., [2], [3]. Of particular interest are the model 
reduction techniques based on asymptotic methods, such as multiple time 
scales and singular perturbations, e.g., [4]-[7] and the excellent survey [8]. 
These papers address the large scale, but finite dimensional or lumped para¬ 
meter, systems; to our knowledge, very little has been done with asymptotic 
methods for DPS with the exceptions: [9] Chapt. 5, [10], and the survey [11]. 

Our Interest here is less in the area of derivation of reduced-order 
models for DPS (although that is an Important topic in its own right); 
rather, we are concerned with the successful operation of the controller, 
designed on a reduced-order model, in closed-loop with the full system. Other 
work on this basic topic may be found, for example, in [2], [12]-[15], and, 
for singular perturbations methods, [16]-[24]; none of these singular pertur¬ 
bations papers deals with DPS, with the exception of [24]. This paper 
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extends the preliminary results on singular perturbations and stability for 
DPS obtained In [24] and obtains an upper bound on the singular parameter 
for stable, closed-loop operation. Our results for DPS are In the same spirit 
as those of [21]-[23] for lumped parameter systems; the method of proof 
Is different due to the Infinite dimensional nature of the DPS problem. 

We consider linear , time-invariant distributed parameter systems with 
the following form: 

v* a Av ♦ Bf , v(0) ■ v n 

, * 0 o.i) 

y « Cv ♦ Df 

where the state v Is In a Hilbert space H with Inner product (...) and 
associated norm ||*|| and the control vector f and observation vector y 
have dimensions M and P, respectively, which denote the number of (inde¬ 
pendent) actuators and sensors. When the system (1.1) Is a distributed para¬ 
meter system (DPS), the dim H ■ » and the operator A is (usually)an un¬ 
bounded, differential operator with domain D(A), containing all states which 
satisfy the boundary conditions of the problemSrdense in H, while the in¬ 
put and output operators B and C have finite rank M and P, respectively, and 
D is a PxM matrix. 

The operator A generates a C Q semigroup U(t) on H. This semigroup 
U(t) is usually exponentially stable, i.e., it has the growth property: 

IIKOII < «o * ot • ‘i 0 (i-z) 

where K Q > 1 and > 0. In physical systems of "hyperbolic"-type > energy dissi¬ 
pative mechanisms make (1.2) true, even for DPS; however, may be quite 
small , as it often is in the case of large aerospace structures [2]. 

In the following sections, we develop the general idea of reduced-order 
models for the systems of type (1.1) and obtain results on the operation 
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of controllers, based on these reduced-order models, In closed-loop with 
the actual system ( 1 . 1 ). 

2. REDUCED-ORDER MODELING UNO CONTROLLER DESIGN 

Let H n and H R be subspaces of the total state space H with dim H N * 

N < • and H ■ H R 0H R . Define the projection operators P N and P R (not 
necessarily orthogonal) and let v N » P R v and v R ■ P R v. This decomposes 
v Into v ■ v N + v R and the system ( 1 . 1 ) Into 


V N " *N V N + *NR V R + B N f 

■ V N< 0) * P N v 0 

(2.1) 

V R * A RN v N * A R v R + B R f * 

• V R<°> ‘ P R v 0 

(2.2) 

y ’ C N V N + C R V R + D f 


(2.3) 

Note that all parameters, , A^ R , 

etc., with the exception of A| 

R* * re 


bounded operators because they Involve projection onto the finite dimen¬ 
sional subspace H^. We abuse notation slightly by writing v R for 3v R /3t. 

Henceforth, we assume that A R generates a semigroup U R (t) on H R which 
satisfies the growth condition: 

• <S R t 

||U R (t)|| < K R e , t > 0 (2.4) 

where K R > 1 and 6 R > 0. In the special case of reducing subsoaces (i.e., 
both H r and H R are A-Invariant), 6 R ■ 6 Q and K R ■ ||P R || Kg where <5g, Kg 
are as given in (1.2). In general, K R and 6 R depend on the choice of H N> 

H r subspaces. The growth condition (2.4) will hold for "parabolic", as well 
as "hyperbolic" systems even though the more stringent ( 1 . 2 ) does not hold. 

The reduced-order model (ROM) for the system is (2.1) and (2.3) with 
C R and Aj^ r assumed to be zero: 

{ ’n ■ A H 5 N * B N f ■ V°) ■ P N v 0 


(2.5) 









Thus, the ROM depends on what choice of subspace H R Is made and what type 
of projection P N Is used (o»* alternatively, what H R Is). The subspace H N 
Is called the ROM subsoace and the subspace H R . the residuals subsoace . 

The terms A^ R v R and v R are called model error and B R f and C R v R are 
called control and observation spillover , respectively. 

3. REDUCED-ORDER MODELING AND CONTROLLER DESIGN VIA SINGULAR PERTURBATIONS 


For certain choices of subspaces H^, H R In Sec. 

produce ROM's such that (2.1)-(2.3) become: 

2, It Is possible to 

V N " *N V N + ^R V R + B N f 

(3.1) 

e V R " A RN + A R V R + B R f 

(3.2) 

y - ^ V N + C R V R + D f 

(3.3) 


where the singular perturbation parameter e > 0 represents some small para¬ 
meter dependence, such as electrical networks with parasltlcs or dynamical 
systems with small masses and time constants (see [25]}; alternatively, e 
may represent a ratio of time-scales In the system (see [7]). This approach 
Is especially valuable when the ROM Is based on vibration modes for a mech¬ 
anically flexible structure and the frequency separation of slow and fast 
modes can be used to provide the decomposition of (3.1)-(3.3). A survey of 
singular perturbation techniques for model reduction is presented In [8]. 

It should be noted that one of the most difficult tasks may be to formulate 
a particular DPS problem Into the singular perturbations format (3.1)-(3.3). 

The parameters such as A^ may Involve regular perturbations In c, as 
well, l.e., 

A„ ■ A® ♦ e Aj (3.4) 

where A|j Is bounded for all e and A^ ■ 0. Also, It Is possible to have multi 
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parameter singular perturbations where several e parameters are present; see 
e.g., [26]-[27], To simplify the presentation, we shall not consider either 
of the above situations here. 

The singular perturbations ROM Is ol*»1ned when e ■ 0 and Is given by 


[ *N " *N *N * *N f » V°> " P N v 0 


(3.5) 


l y “ Z H *N + * f 

where * A N " ^ir A r 1 A^, & N ■ B R - A NR ^ 1 B R . ? N • C R - C R A R 1 A^, 
and E » 0 - C R A R ’^ B R . This Is somewhat different from the usual ROM (2.5) 
due . the parameter correction terms Involving A R "^. In the special case of 
reducing subspaces H R , H R (as In [24]), the model error terms A NR and A RN are 
zero, and A^ • A^, B R * B R , and E R ■ C R ; however, 0 f D. Therefore, even In 
this special case the ROM (3.5) Is not quite the same as the usual ROM (2.5). 

Note that the ROM (3.5) Is finite dimensional (dim v R ■ N) and the para- 
meters [1^, E R , E R , E) may be Identified with their corresponding matrices In 
an appropriate basis for H^. Henceforth, we assume the ROM (3.5) Is con¬ 
trollable and observable ; for any ROM, this Is easy to verify with the the 
usual ran!* tests for finite dimensional systems ([28] Chapt. 11). 

Also, note that, although A R Is usually not bounded for a DPS, A R ”^ 

1$ bounded due to (2.4). In fact, from [28], Theo. 8.9, 


-1 



(3.6) 


The reduced-order controller based on the ROM (3.5) is given by 

'' -V 

z » z ♦ E r f ♦ l^(y-y) , z(0) ■ 0 (3.7) 

, y ■ Ejj z ♦ E f 
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where dim z * N. This finite dimensional controller can be Implemented for 
the DPS, and the gains G N , K N can be designed for closed-loop stability when 
e 8 0; however, successful closed-loop operation is in question when e > 0. 


4. MAIN RESULTS: CLOSED-LOOP STABILITY 


The fundamental question is whether a stable closed-loop reduced-order 
system (e = 0) will remain stable when the same controller (3.9) is used 
with the actual system (3.1)-(3.3) when e > 0 (although small). The answer, 
as given below, is that it will when e is small enough ; bounds on the small¬ 
ness of e are obtained. 

The closed-loop behavior of the total system (3.1)-(3.3) with the 
controller (3.7) is governed by 

w. 


_3E 

_f 


' H n 

H 12 


3 



EWg „ 


- H 21 

H 22- 


1 


2 J 


where w, = [vl eI] T , e N 8 z - v N , w, = v D , and 


'll 


'N N- 
A N + B N G N 


b n g n 




-1 


% s ^ K N C R ’ a nr^ a r ^ A RN + B R V 


(4.1) 


'12 


'NR 


L K N C R " A NR J 


'21 


- A RN + B R g n b r V 


and H 22 8 A^. Note that , H^> Hg-j are all bounded operators, and ^ 
is bounded by (3.6); therefore let 


-1 
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l|H 11 l| <M 1 . ||h 12 ||<m 2 , I|h 21 H < M 3 , 

(4.2) 

11^22 11 — ^4 = * 11^22 ^ 2111—^5 

When e = 0, (3.1) reduces to 

ft 1 * H 1 w-j (4.3) 

which represents the reduced-order design consisting of the ROM (3.5) and 
the reduced-order controller (3.7) in closed-loop. The reduced-order closed- 


loop operator is given by: 


H 1 * H 11 “ H 12 H 22 H 21 


a n + b n g n b n g n 


*N ’ K N C N 


which can be made stable by the choice of gains G^, K^, since the ROM (3.5) 
is assumed controllable and observable in Sec. 3. Therefore, let the transi¬ 
tion matrix U^t) associated with H-j satisfy: 


||U 1 (t)|| < e 


, t > 0 


Furthermore, is bounded: 


11^11 im 6 

’ M 6 = M 1 +M 2 M 5^ K 1 /6 1 
When e > 0, we can write (4.1) as 


w = H w = 


H^/e H 22 /e 


.T . ,T-iT 


where w = [w^ w 2 ] . The primary question to be answered in this section is 
whether the semigroup U c (t) generated by H in (4.6) remains stable for small 
positive e, i.e.. 


|U c (t)lI < K c e 


, t > 0 
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The answer is that it does remain stable and the following theorem — our 
I main result — gives a bound on the smallness required of c: 

THEOREM 1 : Let the controller gains (a^, be chosen so that H-j in (4.3) 
has stability margin 6-j as in (4.4). There exists §n > 0 such 

* that, for all 0 < e < e Q , the closed-loop system (4.6) — consisting 
of the full-order DPS (1.1) and the reduced-order controller (3.7) 
based on the ROM (3.5) — is (exponentially) stable and the controller 

* state z converges (exponentially) to the reduced state v^; this means 
(4.7) is satisfied. An upper bound for Cq is given by: 

6 R . 6 1 

» e 0 < (1 + M 5 ) K r M 2 (M 2 + M 6 ) ' min ^ * Y V 

where y is defined later in Lemma 1, and < c , 5 C , in (4.7), are given 
by: 

1 5 c = min (^ , 6 2 /e Q ) (4.9) 

2 1/2 2 

K c = K 1 k r 0 + ct + or) (3 + 3 M 5 + M^) (4.10) 

* where 

^ = 6 1 - e Q K 1 M 2 M 4 (1 + M 5 ) (M 2 + Mg) 

I 6 2 = 5 r - e Q k r M 2 0 + M 5 > 

o = M 2 /|A| , A = 6-j — 6 2 /£q ' 

> and the constants are obtained from (2.4), (4.2), (4.4), and (4.5). 

The proof of Theo. 1 requires the following two lemmae: 

LEMMA 1: There exists > 0 such that, for all 0 < e £ , the nonlinear 

1 mapping h(L) = H 22 _1 + e H^" 1 L (H^ - H-j 2 L) defined on 

n * {L | 11L - H 21 || < 1} has a unique fixed-point L* = L*(e) 
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in ft (hence, ||L*|| < 1 + Mg). Furthermore, z 2 * Wg + L* w^ trans¬ 
forms (4.1) into 

H-- Ifw. n 


*1 


ez. 


H 1 “ e W 1 


12 


H 22 + e L* H 12 


(4.11) 


-1 


where = H-|g H 22 L* (H^ - H-jg L*). An upper bound for is 

[M 4 (1 + Mg) ( Mg + Mg)]" 1 and we write e T * y[M 4 (1 + Mg) (Mg + M g )]‘ 
where 0 < y < 1. 

LEM1A 2 : There exists e Q such that, for all 0 < e < e Q , 


2 J/2 . -6 c t 

[|!* 1 (t)ll 2 + IU 2 (t)|| 2 ] it c ! c 


[||w,(o)|| z + ||z 2 (o)|| z 3 


1/2 


(4.12) 


1/2 


where K c = K R (1 + a + a ) and a, 6 C the same as in Theo. 1. 

An upper bound forces given by (4.8). 

The oroof of Lemma 1 involves a contraction mapping argument ([29], 
Theo. II. 1.1, p. 24); however, the results of [30] might be used to obtain 
a different (possibly,larger) bound for e-j. It is well known ([29], p. 24) 


that 


lim L. 


(4.13) 


where L fc is obtained via the algorithm: 

L k+1 = h ^ L k^ * k - 0, 1, 2, ... , 

and Lq any member of ft. This algorithm may be implemented to calculate the 
desired L*. The following gives an indication of the convergence rate: 


I Kjt " L *l I - tv I l L o ’ L *11 


(4.14) 


where y is defined in Lemma 1. 


ORIGINAL PAGEU 
OF POOR QUALITY 






-- — -- - To 

The proof of Lenina 2, and hence Theo. 1, makes use of Lemma 1 and 
results from [28]. The proofs of Lenina 1 and 2 and Theo. 1 appear in 
Appendices I-III. 

5. CONCLUSIONS 

Our main result (Theo. 1) provides an upper bound (4.8) on the small¬ 
ness of the singular perturbations parameter (e) which ensures stable closed 
loop operation of a finite dimensional controller with the full-order distri 
buted parameter system. Since the term and Hg-j in (4.6) appear because 
of spillover and model error, the bounds (and hence, the stability) improve 
when these terms can be reduced (.1.e., when M 2 and M 3 can be made smaller). 
Lenina 1$2 taken together form an infinite dimensional version of the 
Klimushchev-Krasovskii result [8]. 

When the distributed parameter system can be put into the singular 
perturbations format, (3.1)-(,3.3), the stability condition presented here 
can be checked with only a limited knowledge of the unmodeled residuals pre¬ 
sent in the full-order system. This makes it possible to synthesize low- 
order controllers for distributed parameter systems via general reduced- 
order modeling techniques and analyze their operation in closed-loop with 
the actual system. Such a result appears to be particularly applicable to 
distributed parameter systems with multiple time-scales or high-low fre¬ 
quency separation, as is often the case in large aerospace structures [2], 
[31]. The general modeling issue, i.e. obtaining a singular perturbations 
format, for large-scale or distributed parameter systems is quite complex; 
see [32] for further discussion. 
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Appendix I: Proof of Lemma 1 

We want to apply the contraction mapping theorem to h to obtain L*. 
Consider any L in n and note that 

IIDI < I|L - h 22 -' h 21 || + ||h 22 ’’ h 21 || < 1 + m 5 (1.1) 

Also, ||h(L) - H 22 _1 H 21 || - e||H 22 _1 L (H„ - H, 2 L)|| 

< e M, ||L|| (M 1 + M 2 ||L||) 

< e «4 (1 + M s ) (Mj + M 6 ) (1.2) 

because Mg * + Mg Mg in (4.5). Therefore h(L) is in fl when 0 < e <. 

where 

0 A ■ [H 4 0 ♦ M 5 ) (M 2 + Mg)]* 1 (1.3) 

because ||h(L) - H^' 1 H 21 || <1. 

Consider L and L 1 in n and, from (1.1), 

11h(L) - h(L')|l ■ c ||H 22 _1 [L(H n - H 12 L) - L'tH,, - H, 2 L')]|| 

1 £ M 4 (I |L(H n - H^g L) - L(H^ - H^g L')j| + 

11L(H ]1 - H 12 L*) - L 1 (H ]1 - H 12 L’)||) 

< e m 4 (11l h 12 || + !|H nl - h 12 L'll) ||L - L'M 

< e M 4 (2(1 + Mg) Mg + M ] ) | |L - L*|| < ||L - L'M (1.4) 


when 0 < e < Eg where 
e B " 


-1 


[M 4 (M 1 + 2(1 + Mg)Mg] 


(1.5) 
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But, from (1.3), 

e A ■ C(1 + Mg) M 4 (M, + M 2 (1 + Mg)))* 1 

< [0 + H 5 ) M 4 (M 1 + 2 M 2 (1 + Mg )))" 1 

■ +M sV 

Therefore, e A < (1 + Mg) e A * Eg-, hence, we choose and, from (1.4), 

h is a contraction on n. By the contraction mapping theorem ([29] Theo. 

II.1.1, p. 24), h has a unique fixed point L* in ft for any 0 < e < e-j, i.e. 

L* - h(L*) - H 22 -1 H 21 + e H 22 ' ‘ L*(H n - H, 2 L*) (1.6) 

and, from (1.1), 

||L*||<1+M 5 (1.7) 

Note that - H 12 L* s - e ^ and, by substitution of z 2 into 
(4.1), obtain w-j + H 12 (z 2 - L*w-j) s (H-j - e w 1 + H 12 z 2 and 

e z 2 * e(w 2 + L* w-j) 

* H 21 w 1 + H 22 w 2 + e L* ((H ] - e W.,) w-, + H 12 z 2 ) 

- H 22 [h(l*) - L*] + (H 22 + e L* H 12 ) z 2 

• (H 22 + e L* H 12 ) z 2 by (1.6). 

This is the result desired. # 








15 






Appendix II: Proof of Lemma 2 

Assume 0 < e < e Q < e-j and, from Leimia 1, obtain 
^ - e w 1 + H 12 Zg (II.1) 

ez 2 » (H 22 + e L* H 12 ) Zg (II.2) 

Choose t * t/e and Zg(t) ■ z 2 (t). Therefore, from (II.2), 

A 

3T * H £ ■ e *2 ■ < H 22 + eL * H 12> h ■ 

Hence, from H 22 - A R and (2.4) and [27], p. 215, 

A 

IlijWM < K r t ? 1 |z 2 (0)| | (II.3) 

where 6 g * 6 R - e Q (1 + Mg) M 2 K R because e ||L*H 12 || < e Q (1 + Mg) Mg. 
Choose 

e 0 < 6 R C;i + M 5 } M 2 K R ]_1 (II,4) 

so that &2 > 0* From (H*3) and ||z 2 (t)|| * ||z 2 (t)||, obtain 

4 * 

l|z 2 (t)|| < K R e e IU 2 (o)|| 

t 

<K R e E 0 | |z 2 (0) | ( (II.5) 

because e < Eg . 

Choose Cq < 5 1 [(1 + Mg) Mg M 4 (Mg + M g ) K ,]’ 1 (II. 6 ) 

where $-| and K 1 are given In (4.4); assume e Q satisfies (II.4), as well. Let 
5^ * - Eg K.j M^ M^ (1 + Mg) (Mg + Mg) which is positive by (II. 6 ). Since 

II W 1 11 < Mg M 4 (1 + Mg) (Mg + Mg), we have, from (4.4), (11.1), (11.5), and 
[28], p. 215, 









VI 
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|w,(t)|| <K, e’* 1 * [11^(0)11 *n 2 f t e* 1 * I|z 2 (s)l|ds3 


-6,t 


< ^ e 1 C||w 1 (0)|| + H 2 K r /* e As ds ||z 2 (0)||] 


■6,t 


< K, e'° n Ul«,(0)|| +Hj K„ £f-l ] 


(11.7) 


A 2 

where A ■ 6, - — . 

I E*» 


-? t 
E 0 


If A > 0, then (II.7) Implies Hw^tJII < ^ e u (||w 1 (0)|j + 


”2 


- 5 lt 


F* I |z 2 (0) ||). If A < 0, then (II.7) Implies ||w 1 (t)|| < ^ e 

A 

H K 5 

(i Iw^O) I! + 11z 2 (0)I!)* In either case, if 6 Q • min (Sj, ~), then 

||w 1 (t )11 < ^ e C ( 11 »» 1 ( 0)11 + l|z 2 ( 0 )||) (II» 8 ) 

and, from (11.5) and (11.8) and and K R >_ 1, 

-6 c t2 


Mw^tJir ♦ ||z 2 (t)ir < (K, Kr e c ) (||z 2 ( 0 )ir + (l|w 1 ( 0 )|| 

M 2 
w 


M, , -6 t 2 2 

+ W ll*2C°)ll) ) 1 < K i k r « ) 0 + a + a ) 


(I|w 1 ( 0)|| 2 ♦ ||z 2 ( 0 )|| 2 ) 

Mo 

where o ■ jjj- and the easily verified inequality: 

b 2 + (a + b ) 2 < (1 + a + a 2 ) (a 2 + b 2 ) 
is used with a * ||w^(0)|| and b * ||z 2 (0)||. Therefore, 


1/2 


-6 c t 


w^tjir + n* 2 (t)ir) iK c e c (iiw 1 (o>ir + 1iz 2 (°>ir> 


1/2 


(II*9) 






and this will meet all the requirements. However, we can refine (11.10) 
further because 0 < y < 1 and MgtMg + Mg]" 1 < 1; consequently, (II.10) 
becomes: 



which Is the desired result. # 
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Appendix III: Proof of Theo. 1 

Assume 0 < c < c Q , where e Q given by (4.8). Consider thet, from 


■ W 
w * 

♦ L* w 1 , 


m. . m 


- m 

' 

- Q 

"l 

CM 

M 

_ % 


- w 2. 


where 


'i o' 


Q -i . 

‘ I o' 


and 


.L* I, 



.-l* I. 


From Lemma 2,(11.9) and (III.1), obtain 

INt)|| < l!Q* 1 || 11Q11 K c e" 6ct | |w(0) 11 (III.2) 

But ||Qw|| - 4|w 1 || z ♦ ||w 2 + L* w^l* 

? 1/2 7 

< (1 ♦ ||L*|| + |!L*|! 2 ) ||w|| Z (III.3) 

by use of the easily verified Inequality: a 2 + (b ♦ aa) Z < (1 + a + a 2 ) • 

(a 2 + b 2 ) with a ■ ||w^ j|, b ■ ||w 2 ||, and a « ||L*||. Therefore ||Q|j < 

(1 ♦ | jl*| | + | |l*| | 2 ) 1/2 and, similarly, | IQ* 1 11 £ (1 ♦ !|L*|| ♦ !|L*M 2 ) 1/2 . 
Note that, from Lemma 1 (1.7), ||L*|| < 1 ♦ M g . When these Inequalities are 
substituted Into (III.2), It Is clear that the desired result is obtained. # 
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ABSTRACT 

Using s singular perturbation formulation of the linear time-invariant distributed 
parameter systea, we develop a aethod to design finite-dimensional feedback compen- 
saters cf any fixed order which will stabilize the infinite-dimensional distributed 
parameter systea. The synthesis conditions are given entirely in tens of a finite- 
diaensimal reduced-order aodel; the stability result depends on an infinite-diaen- 
slocal version of the Klimuahchev-Kraaovskii leasaa developed in 1101 * 


1.0 INTRODUCTION 

Many engineering systems exhibit a distributed paraacter nature and must b- des¬ 
cribee by partial differential equations. Exaaples of such distributed paraacter 
systecs (DPS) include heat diffusion 4 chemical processes, wave propagation, & 
mechanically flexible structures. Various aspects of the control of DPS have been 
considered in, for example, [1] - 15]; our experience in DPS has been shaped by 
applications in large aerospace structures [6]. 

The state spaces for DPS have infinite dimension; sc, at best, reduced-order models 
must be used in controller synthesis. However, the closed-loop stability of the 
infinite dimensional DPS with a finite dlaensicnal feedback controller becomes a 
fundamental issue. The synthesis of finite dimensional controllers for DPS snd the 
analysis cf their closed-loop stability by singular and regular perturbation techni¬ 
ques bare been our main areas of emphasis [7]; this theory has been developed with 
flexible structures and other highly oscillatory DPS applications in mind. 

Even in large-scale, lumped parameter systems, such as electric power distribution 
networks, it is necessary to perform model reduction and reduced-order controller 
synthesis and to analyze closed-loop stability. The use of asymptotic methods, 
especially singular perturbations, has been very successful in this regard (e.g. 

[S]). Va have extended certain of these singular perturbations methods for DPS 
to provide estimates of stability in an infinite dimensional setting [9] - [10] and 
applied them to mechanicall** flexible structures [11]. In this paper we will use 
these singular perturbation results to analyze the general finite dimensional com¬ 
pensator for linear DPS. 
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We consider linear, time-invarlent DPS with the fora: 

( v t » Av + Bf, v(0) • Vq 

y - Cv 


( 1 . 1 ) 


where Che state v(c) is in a Hilbert space H (called the state space) for all 
t*0; the state space has inner product (*,•) and associated norm |]-1|. The 
control vector f and the observation vector y have dimensions M & P, respectively, 
which denote the number of (independent) actuators and sensors in use; thus, rank 
B*M and rank C»P. The operator A is an unbounded differential operator with 
domain D(A) dense in H and A generates a C. semigroup U(t) ou H. Many practical 
DPS can be stated in this form [7]. 


2.0 REDUCED-ORDER MODELING OF DPS: A SINGULAR PERTURBATIONS FORMULATION 


Since the state space B of the DPS in (1.1) is infinite dimensional, we must obtain 
a reduced-order model (ROM) upon which to base the finite dimensional controller 
design. In general, this is done by selecting a finite dimensional subspace. 

(with dim. % » N < «) contained in D(A). This subsp-ce Hu is the ROM subspace ; 
its complement Hr is the Residual subspace, and together they decompose the state 
space: 

( 2 . 1 ) 

We define the (not necessarily orthogonal) projection operators Py and ?r and let 
vy » Pyv and vr * Prv. These decompose the state v ■ vy + vr and the DPS (1.1). 

In this paper, we will make two basic assumptions on model reduction processes 


(1) the subspaces Hy, Hr are reducing subsnaces (i.e. they are both A- 
invariant); 

(2) the subspaces Hr, Hr can be chosen to give a singular perturbations formula¬ 
tion with small parameter eaO and an exponentiallv-stable residual subsystem. 


These assumptions yield the decomposed DPS U.l) in the following form: 


av 


N 


3t 


S V N + V 


(2.2a) 



*R V R + V 


C N V N * C R V R 


(2.2b) 

(2.2c) 


where Ay * PyAPy, By * PyB, Cy ■ CPy, etc. and Ar generates the Co semigroup 
Ur( t) * PrU( t)P r with exponential growth condition : 


|| U R (t)|| < e“V, tiO (2.3) 

where Kr a 1 and Or > 0. All operators except Ar are bounded. The terms Bftf & 
Cr vr are called, respectively, control & observation spillover ; they represent 
the interconnections through which the controller can affect the residual sub¬ 
system. 


The ROM of the DPS, in this formation, is obtained by setting £-0: 














, 9«‘ - ** *K + V (2 . 4) 

■> *S v « + V 

where D[j ■ - CrAr Hr & Ar ^ is a bounded operator due to (2.3). The ROM is a 
finite dimensional system and the parameters (A^,B^,C^,%) may be identified with 
matrices by choosing a basis in %; this will be done wherever necessary. We 
shall assume , henceforth, that the ROM is controllable and observable (easily 
verifiable conditions exist for finite dimensional systems). The ROM is completely 
determined by the choice of the subspaces Hu & Hr. 

3.0 FINITE DIMENSIONAL COMPENSATORS FOR DPS 

The form of the finite dimensional feedback compensators used here will be the 
following: 

f f " L n- V + L i2 2 

i (3.1) 

l i- Hl y + L 22 2 

where the compensator state z has dim. z * a < N < ® and L]j,Li2»L2i»L22 are 
matrices of appropriate sizes. We say the compensator is output feedback when 
ct«0. The order a of the compensator is assumed to be fixed at some acceptable 
value which reflects the available capacity of the on-line computer being used to 
implement (3.1). In [10], the compensator order was a » N. 

The compensator design is synthesized as though the ROM (2.4) were the full DPS 
(1.1), i.e. as though £*0. Let 

r -V¥S < 3 - 2 > 


r*» ° 


Hi “ (I M " L u V L 11 
Hi * (I M " Hi V L 12 
Hi ’ Hi ( H + H Hi } 

H2 “ L 22 + Hi °N H2; • 


Let the composite state qjj - 
(for e-0) : L 2 


in H« x R ; from (2.4) & (3.1), this satisfies 


The following theorem gives the conditions under which a stable design can be 
synthesized: 


Theorem 3.1 : If the ROM (A^,3i5,Cjj,Dj5) is controllable and observable and 
M+P+aiN+1 (3.4) 

then L may be chosen so that F given by (3.2) has any desired pole locations in 


I 












» 


I 


I 


t 


» 


I 


I 


» 




» 


chc complex plane. 


The proof follows easily from results in [12] or [13] and depends on_the finite 
dimensionality of Che closed-loop system when e»0. Since Che gains L can be 
synchesized, Che compensacor gains Lll,Li2>L2i»I-22 can be found; however, we muse 
assume that jj || is sufficiently small to make boch (I>j- L^ Dji) -1 & (Ip+DN^il)” 1 
exisc. 


4.0 CLOSED-LOOP STABILITY FOR THE DPS WITH A FINITE-DIMENSIONAL 
COMPENSATOR 

Let e be such that 0 < e < eg, aQ d consider the closed-loop system consisting of 
the infinite-dimensional DPS (2.2) and the finite-dimensional compensator (3.1) . 
Let the closed-loop state w in H ■ Hm xR Q x Hp be given by 


where wjj 


■j V s|is in H>j x R a 4 wj ■ vj in H R . The norm on fi is defined by 

1/2 


( v. 


N 


i*ir + !iv R ir) 


From (2.2) & (3.1), the closed-loop state satisfies: 


St 

9w_ 


H 11 W N + H 12 W R 


e St E 21 W N + H 22 w R 


where 


E 11 


12 


h + ViiS 


L 2lS 


B1 C 
N 11 R 


L L 21 C R J 


B N L 12 


(4.1) 

(4.2a) 

(4.2b) 


E 21 tV'llS B R L 12 5 


H 


22 


*R + B R L U C R 


But, (4.2) is in the form of the infinite dimensional KlLmushcnev-Krasovskii len™« 
whose proof appears in the revised version of [10]; estimates on the upper bound 
eo are presented there, as well. From the infinite-dimensional K-K lemma, we 
obtain the following closed-loop stability result : 

Theorem 4.1 : If E 22 generates an exponentially stable Cg semigroup and the finite 
rank operator 

8 ■ H U - H 12 a 2 2 ' 1 *21 (4 - 3) 

has all its poles in the open left half-plane, then there exists positive to (whose 


I 







value can be estimated from the model data) such that, for all 0 < E < cq , (4.2) 
is stable. 


This theorem says that under certain conditions the finite dimensional controller 
based on the ROM (£-0) will continue to stabilize the infinite dimensional DPS 
when C is "small" but non-zero. The "smallness" required of £ can be estimated 
from bounds on the control and observation spillover terms, the constants Kr,Or 
in (2.3), and the BOM data. 


The following lemma greatly simplifies Theo. 4.1: 

Lemma 4.2 : H • F ■ A^ + BjjLCj, (4.4) 

The proof appears in the Appendix. Conseo.uently, Theo. 4.1 may be restated as 
our main result: 


Theorem 4.3: If the stable compensator synthesis conditions of Theo. 3.1 are 
satisfied for the BOM (2.4) and if the control and observation spillover coeffi¬ 
cients £ * ijBpj! & F » I) C R j| are sufficiently small that 


~> g r 




'li 1 


(4.5) 


where Kr & Or appear in (2.3), then there exists positive £o such that for all 
0 < e ^ £q the same finite-dimensional compensator (3.1) stabilizes the infinite- 
dimensional DPS (1.1). 


This follows from Lemma 4.2 and use of the semigroup-generator perturbation result 
(I14], Theo. 10.9) applied to 

E 22 * \ + B R L 11 C R 

because BqL i-| is a bounded operator. 


5.0 CONCLUSIONS 

Our main result (Theo. 4.3) gives conditions under which a general finite-dimen ¬ 
sional compensator , based on a singular perturbations reduced-order model, will 
stabilize an infinite dimensional (linear) distributed parameter system . This is 
an extremely useful result since all practical feedback compensators must be 
finite dimensional in order to be implementable with on-line digital computers. 
The result is valid for large-scale lumped-parameter systems, as well. 

The most difficult assumption to satisfy is the choice of reduced-order subspace 
to achieve a singularly perturbed formulation (2.2) of the distributed para¬ 
meter system (1.1); the other assumptions are reasonably easy to satisfy. Modal 
methods have worked for mechanically flexible structures, e.g. [11]. A general 
discussion of this modeling difficulty is given in [15]; it is a fundamental 
problem in large-scale or distributed parameter systems and should not be over¬ 
looked. 
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APPENDIX: PROOF OF LEMMA 4.2 


H “ H n - h 12 h 2 


-1 


21 


A^ + Z Lj^Cj, 

L 21 (I ? " C R E 22 -1b r L 11 )C N 
11 C R H 22 ~R 


where Z • 1^ - L„,C„H„„ J 7B 


H 22 " \ + B R L U C R‘ 


*n + b n l c n 
S + b nHi c n 
ihi c n 


B 2 L , 

N 12 

L 22 " L 21 C R H 22 ^rHz 


Vl2 


12 J 


I 








I 


« 


f 


f 


t 


r 


i 


t 


» 


> 


where ■ 

(I M ’ L llV L 11 

C 12 " 

(I h - hlV' 1 L 12 

“21 " 

L 21 (I P + D /ll ) 

T m 

l 22 

L 22 * L 21 D N^12 

& d n - ■ 

• W 1 h 


Therefore, H ■ F if end only if 
(.) I - (Ij, - L U D S )- 1 and 

<>» Vu'\ ' - V 1 * - HlV 


However, (a) (b): 

’ °N CI M - L llV ' 1 ’ - V - D N + °2hlV2l'\ 

* C th' 1( b. - b * l u c r h 22 'H 

■ - h i u c *y\ 

" C R A R (H 22 ' 3 R L 11 C R )H 22 E R 

■ C R-V 1(A R )E 22' 1 *S ’ C R H 22' 1 5 R 

Thus it remains to prove (a) which we recast as: 
id) (Ijj + QV)" 1 - ^ - Q(WQ * I r ) _L vJ 

where Q « L^C R and V ■ B R are both bounded (finite rank) operators. 

Since (WQ + 1^) (I R + 1 B^L^C^) 1 - ^ and H 22 ^ is assuoed t0 exist 

for snail enough spillover bounds, therefore 

= - Si * Q(VQ + V 1 * 


where the required inverse exists. Thus, to prove (a) we need only prove (a'); 
however, this follows directly: 

(1^, + OK) (1^ - Q(WQ + I R ) _1 V) 

- y + QK - Q(WQ + I R f 1 V 

- Q[WQ(WO + I R ) -1 ]W 

- Ijj + QW - Q(WQ + I R ) _1 W 

- Q[I R - (WQ + 

- ly + QV - Q(VQ + I R ) _1 V 

- QV * Q(VQ * I R ) _1 W - 

•• (Ijj + QV) -1 - In - Q(WQ + I^r 1 V 


and (a') is proved. 
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Abstract 

Large Space Structures are treated as a special 
class of highly oscillatory distributed parameter sys¬ 
tems. Practical controllers will need to be finite di¬ 
mensional; however, this calls the closed-loop stabili¬ 
ty Into question due to spillover interactions. Recent 
results, obtained using regular and singular perturba¬ 
tions theory, give stability bounds for closed-loop 
control of this type of distributed parameter system; 
this paper surveys these results and their implications 
for large space structures. 


Ignore these external disturbances to simplify this 
discussion and consider the control forces as generated 
by M Independent actuators: 

M 

F (x,t) ■ B.f ■ I bi(x) f.(t) (1.3) 

o o 1 1 

Observation is done with P Independent sensors whose 
output vector Is given by 

y(t) • C 0 u(x,t) * E Q u t (x,t) (1.4) 


1.0 Introduction 

Large eerosoace structures, such as satellites and 
spacecraft, are mechanically flexible structures which 
will usually recuire active (feedback) control systems 
to suopress'vibrations and satisfy stringent require¬ 
ments for pointing accuracy, orientation, snape, and 
station-xeepirg. The problems and current directions 
ir. large aerospace structure control nave been surveyed, 
for examole, in I"0-[3]. 


1.2 State Space Description 

The OPS (1.1), and (1.3)-{1.4) can be rewritten in 
many state space formulations (not all of tnem are equiv¬ 
alent); however, the most natural one is tne state 
v(x,t) ■ [u(x,t), u t (x,t)] whith satisfies: 


Sucn structures are hlgnly oscillatory, distributed 
paramete* - s'Stems (DPS). The viewpoint we have been 
pushing for"DPS *n general is a practical control the¬ 
ory, i.e. one that includes such practical constraints 
as a finite number of (localized) control devices and 
finite-dimensional (implementable) controllers; this 
has been discussed in [4], Other helpful references on 
aspects of DPS t.oeory are r 5]-[8j. In this paper, 
large aerospace structure control will be viewed in the 
frameworx estaolisned in [a], and recent stability re¬ 
sults us:nc -ecu'ar ano singular perturbations tech- 
ricues [9j-IlCj will oe examined in the large space 
structure context. 

,aroe Ae-cscace Structure Description 

Large aerospace structures may often be modeled by 
the following DPS description: 

u..!x,t) - :_u.(x,t) * A u(x,t) « F (x.t) (1.1) 

ww -J w V 0 

wnere F„;x,o) is the displacement of the structure 
< R n off its e-uilibrium position in space. The 
damping and stiffness differential operators D$ and A$ 
are oefinec on a domain 2(A«} which is a subspace of 
the Hilbert Space H e • L*(fl) with Inner product (•,■)«; 
D(A.) contains all smooth functions in H 0 which satisfy 
tne structural boundary conditions. Furthermore, the 
operators satisfy: 


with the state space being the Hilbert space 
H ■ DfAg^ 2 ) x H 0 with the energy norm 

Iiv!i 2 • Mu t ;; 0 2 ♦ ii* 0 1/2 ui; 0 2 n.«) 

and the operators (A, B, C) determined from (1.1) and 
(1.3)-(1.4). The input and output operators 3 and C 
have finite rank, and the system operator A has domain 

0(A) ■ D(A 0 ' /2 ) x D(A g ) and generates tne C 5 semigroup 
U(t). 

The camping present ir. the structure gives tne sys¬ 
tem a small amount of overall Dissipation. It is tne 
job of the finite-dimensional feedback controller to re¬ 
distribute this dissipation to critical structual mooes 
to improve their performance without creating instabil¬ 
ities in any residual modes. 

2.0 Mocel-Reduction and Controller Svrtresls 

The synthesis of a finite-dimensional controller is 
based on a reduced-order (finite-dimensional) model 
(ROM) of the DPS. In general this Is obtained by de- 
composing the original state space H Into the ROM sub¬ 
space and the residual subspace h r : 

h-h n £h r (2.1) 


0 < < k 

”0 *< " 2 '< J k’k 


(1.2a) where dim h N < «. Then projections v N • P,.v and 

v R » P.v are defined for these suospaces, and the DPS 
(1.5) is decomposed: 


where :. (x) are tne mode shapes . tne mode *reeuen- 
cies . anc ; the modal danb'nq (which is generally 
thought tb oe sr.all for large aerospace structures). 



Although the force distribution r 0 (x,t) in (1.1) is 
tne sum of :ist.rbance and control forces, we will 


V N V N “R V R 


(2.2a) 
(2.2b) 
(2.2c) 
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wh« r « A* • P n AP n , A.. r . P n AP r , 8 n • P^. etc. Th. 
t«rms A nr v R and A RN V N are called model error , and 
B R f, C R v R art called control and observation spillover . 

respectively. Th* parameter c > 0 appears In some 
cases for appropriate choices o? the subspaces H<< and 
Hr; when It appears and is small, w* are sealing with a 
singular perturbations formulation of th* OPS. In many 
eases c •l for the choice ot subspaces and this Is a 
regular perturbations formulation . Both of these form¬ 
ulations wii) be dealt with Here. A basic, and reason¬ 
able, assumption made In both formulations Is that A R 

generates a stable C« semigroup . 

2.1 Regular Perturbations ROW and Controller Synthesis 

When t • 1, th* regular perturbations RO M is ob¬ 
tained from (2.2) by ignoHng tn* residua' state v R : 

( v N ) t " A N V N * B N f ( 2 - 2 *) 

y • Cj, v N (2.3b) 


with A® - dlag [«*.wj], ojj • 2 dlag .5^], 

and B r , cjj have terms Involving the Inner products in 

Ho of actuator or sensor influence functions with mode 

shapes .: N - Corresponding parameters appear for 

th* residual subsystem (A R , 8 R . C R ) In (2.2). 

The modal controller associated with this modal ROM 
can be obtained d1r«ct)y using th* parameters (2.7) In 
(2.4). Modal controllability and observability condi¬ 
tions are available and easily verified, e.g. [3]; con¬ 
sequently, th* controller synthesis may be carried out 
by standard finite-dimensional techniques. Not* that 
any subset of critical modes may be used In this ap¬ 
proach. 

2.2 Singular Perturbations ROM and Controller 
Synthesis 

7h* singular perturbations formulation of (2.2) 
occurs wnene > o anc small. The singular perturottIons 
ROM Is obtained from (2.2) by setting c • o and is 
given by 


This is a finite-dimensional subsystem (since 
dim Hjj < ■); consequently, the following finite-dimen¬ 
sional controller may be synthesized: 


f ■ S ; N 

(2.4a) 

( V't * A N * V * 

(2.4b) 

5 ■ C N 5 N 

(2.4c) 


y * V N + ^ f ( 2 -8b) 

where « N • A^ - ^ R A R A RN> ^ - A NR A R 8 R , 

* C N * C R A s' A 9N’ and Z • - C R A r B r . This Is 
clearly different *rot: the ROM oDtaineo in (2.2). How¬ 
ever, it is a finite-cimenslonal ROM and the controller 
associated with it Is given by: 


The controller gains ^ and are designed to 
staoilize and improve tne performance of th* ROM (2.3) 
in closed-loop with th* controller (2.4). Of course 
the closed loop stability of the actual OPS (1.5) with 
the controller (2.4) will be a principal question; reg¬ 
ular perturoation theory is used to answer this ques¬ 
tion in [9], 


For large soace structures, tne most sobu'ar model 
reduction technique nas been tne modal one, i.e. 


P.,v • 


« u 


. H U t 
- ' 


(2.5a) 


(2.5b) 


f * y R (2.9a) 

(v N ) t * \ y N * # N f * t N (y-y) (2.9b) 

y -T N 0 N -?f (2.9c) 

This controller differs from that in (2.4); however, 
tne gains anc 2^ are oesignec by standard finite- 
dimensional techniques to stabilize and improve tne per¬ 
formance cf the RCM ;'2.8) in closed-loop with the con¬ 
troller (2.9). Again, tn* closed-loop stability of tne 
reouceo-oroer controller (2.9) with tn* actual DPS wn*n 
c is nonzero Is in question. To attempt to answer this 
fundamental question a singular perturbations stability 
theory was developed in [10] for DPS of th* form (2.2) 
with finite-dimensional controllers (2.9); this theory 
gives estimates of tne toleraele size of e. The corre¬ 
sponding singula- oe'turoations theory for finite dimen¬ 
sional systems is surveyed in [H]-[12]. 


where P“ u • 2 u,o,, u * : u^Sl. and u,, 

N k»l * * R k»N*l * * 
satisfies (by using (1.2) in (1.1)): 

„ 2 u • F ■ (f .a 
k k « O k 


♦ 25 fc u 


k k 




In this formulation, the projections P^ and Pr 
are orthogonal n H and tn#y •’reduce" A, i.e. tne 
model error te-'is A NR • 0 and A^ ■ 0; of course, this 
depends on kr>,«1ng th* true structure modes rather 
than aporcximations o' them. Sc the 'egular pertur¬ 
bations ROM 2.3) o' the structure nas parameters: 



For large space structures, as in most applications, 
tner* is no preferred singular perturbations formula¬ 
tion. The general approach has been applied to spact 
structures in [12]. However, as pointed out in [14], 
the hardest task in the application of th* singular 
perturbations approacn is casting the problem in an 
appropriate singular formulation (2.2) wner* tn* small 
parameter c will neve a connection with physical param¬ 
eters end it will be possible to verify whttner t is 
sufficiently smell or not; this is true for applica¬ 
tions of th* aporoacn to DPS, as well as lumped param¬ 
eter systems. 

Here «* will mention two directions - time scaling 
anc 'recuency scai'.-g - which lead to singularly s*r- 
turoed formulations of the mocalmodei eiscusseo 'or 
structures at th* ere of Sec. 2.1. 

In the time-scaling aooroecn (see [12] for a general 
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discussion), we assume that the structure modes of 
(2.6) seperite Into vlorstlons on two time scales • a 
slow one and a fast one - whose ratio Is c. I.a. (2.6) 
becomes: 

“k * 2 'kVk * J k u k ‘ %' 1 i k i N l2 - 1Cl) 

c 2 ^ * 025^1^ ♦ u£u k • F k , k ^ Hd (2.10b) 


Some relabeling of modes might be required to obtain 
(2.10). Our associated singular perturbations formula* 
tlon Is obtained by taking v N as In (2.6a) and, In 
place of v R , using: 


where v. is as In (2.5b). This yields a singular form¬ 
ulation*^.2} of tne modal model of the structure with 
c being the ratio of the time-scales; the mooal param¬ 
eters are the.same as those given In (2.7) with the ex¬ 
ception that C R « [C$ cEg] replaces C R in the residual 

subsystem. The ROM (2.8) is obtained by setting c • 0, 
and the corresponding controller (2.9) can be synthe¬ 
sized; this controller will differ from the original 
modal controller (Z.a). 


In the frequency scaling approach (see [15] for a 
general discussion), we assume that c ■ 1/w^ for N 

sufficiently large, l.e. the residual structural modes 
nave sufficient frequency separation from the ones con¬ 
tained in tne ROM. In this case tne associated singu¬ 
lar perturbations formulation (2.2) Is obtained by 
taking v R as in (2.5a) and. In plade of v R , using: 


? 



o! 


l/c 


( 2 . 12 ) 


where v. is as in (2.5b); the modal parameters are tne 
same as*those in (2.7) with the exceptions: 


h * > !o *o ' * n<i ' : r * [« 2c r c£ RJ r,9l4Ce 

: .* H ; "‘'R - 

A R anc C., wnere • diag [s^, ...], 

5r * 2 -1*1 a N*l ’ and - VlW for 

k > N*l. Something interesting happens here when we 
seT c • C to obtain the SOM ’2.8) ano corresponding 
controller (2.9): since ig pecomes zero wnen c ■ 0 
(and * 0 , Ap f . • 0 because of tne use of structure 
modes). tne singular perturbations ROM (2.8) is exactly 
tne same as the regular perturbationse one (2.3), i.e. 

3^ ■ Aj,, T N • Sg, • C„, and 7 • 0. Consequently, 

for the freouenev scaling approach the controller (2.}) 
i s exactly tne same as tne original mooal contrclTpr ] 
when sufficiently many mooes are controlled, however, 
the singular perturbations stability results of [1C] 
would seem to yield different conditions on the size of 
c • 1/u^ than the regular perturbation results of [9]. 


controller synthesis techniques for OPS have appeared 
In [9]-[10] with corresponding stability results ob¬ 
tained from both singular and regular perturbations 
formulations of tha problem. These have been surveyed 
here by considering their Implications for modal con¬ 
trol of large aerospace structures. Both singular and 
regular perturbations formulations yield useful results, 
but the singular formulation seems extremely versatile 
due to the variety of singular formulations available 
In large space structure applications; two of these 
singular formulations - time and frequency scaling • 
ware discussed here. 
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On-line adaptive control is essential for Large Space Structures 
(LSS).where the modal parameters are poorly known, due to modeling 
error, or changing, due to variable configurations. It Is especially 
Important that such adaptive controllers produce stabilizing controls 
during adaptation due to the small damping present In LSS. However, 
any such controller must be based on a reduced-order model of the LSS. 
The spillover from the unmodelled residuals, as well as the modeling 
error, can deteriorate the performance of the adaptive controller and, 
If uncompensated, this spillover can defeat the whole purpose of the 
adaptive control. 

This paper Investigates adaptive control for LSS using direct and 
indirect schemes and points out the mechanisms whereby observation 
spillover can upset the stability of the controller. The framework for 
nonadaptlve control of LSS Is reviewed and many of the generic problems 
of adaptive LSS control are pointed out within this framework. These 
generic problems must be overcome for successful operation of adaptive 
LSS control. 


1.0 INTRODUCTION 


This paper deals with the basic problems Inherent In adaptive con¬ 
trol of large space structures (LSS), such as satellites and spacecraft, 
where the structural parameters are poorly known or slowly time-varying. 

With the advent of the Space Shuttle Transportation System, It has 
become possible to conceive very large spacecraft and satellites which 
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could be carried Into space and deployed, assembled, or manufactured 
there. Such LSS would serve a variety of civilian and military needs 
Cl], [2], Including electrical energy generation from the solar power 
satellite - a structure nearly the size of Manhattan Island - to be 
constructed In space and operated In earth geosynchronous orbit [3]. The 
control technology needs for such LSS have been discussed In a variety 
of articles, e.g., [4], [5], and the developing LSS control theory and 
technology has been surveyed, for example, In [6]-[10]. 

The size of these structures, their low rigidity, and the small 
damping available In lightweight construction materials combine to make ■ 
LSS extremely mechanically flexible. In theory, LSS are distributed 
parameter systems whose dimension Is Infinite; however. In practice, 
their dynamics are usually modeled by large scale systems based on 
approximate elastic mode data. Active control schemes for LSS are often 
required to meet stringent requirements for their shape, orientation, 
alignment, and pointing accuracy. Such active control Is limited by the 
capacity of the on-board control computer, the modeling Inaccuracy In cur¬ 
rent finite element computer codes for analyzing structural dynamics, and 
available control devices (actuators and sensors); therefore, the con¬ 
troller must be based on some reduced-order model (ROM) of the LSS. 
Fundamental problems of LSS control Include: 

(1) selection of appropriate modes to control for desired system per¬ 
formance; 

(2) development of ROM for analysis and controller design; 

(3) computation of system model and control parameters; 

(4) controller design with multiple distributed actuators and sensors; 

(5) the number and location of sensors and actuators for efficient 
control; 

(6) the effect o', and compensation for, residual (unmodelled) modes 
and modeling error on the closed-loop system performance; 

(7) adaptive and self-tuning controllers for LSS with poorly known or 
changing parameters snd configurations. 

Item (7) is the basic topic of this paper but It must be considered in 
the context of the other Items with which It Is completely Intertwined. 

The need for adaptive control In LSS arises because of Ignorance of 
the system and changing control regimes . The former is due to 









(a) Ignorance of the system structure and order, and (b) Ignorance of the 
system parameters; the latter occurs because of changing configuration 
of the LSS. Changes In configuration may be due to construction In-space, 
thermal distortion, or reorientation of subsystems, e.g., rotating solar 
panels or sunshlelds; these changes usually produce slowly time-varying 
parameters. Ignorance of the LSS system structure and order Is due to the 
fundamental problem of modeling a distributed parameter system, e.g., 
faulty physics, reduced-order models, and Ignored nonlInearlties; this 
means that the order of the ROM Is lower than that of the actual LSS. 
Ignorance of the system parameters , while directly related to the system ' 
structure, Is due to the Inherent modeling error present In even the best 
structural analysis computer codes and to the limitation of testing huge, 
lightweight LSS on earth; this produces constant but poorly known system 
parameters. There Is a very clear need for an adaptive LSS control metho¬ 
dology that can begin with the best available computed parameters and 
self-tune Its way toward the correct parameters while stably controlling 
the LSS and, possibly, adapting to variable configurations. 

A great variety of adaptive control schemes exists for lumped param¬ 
eter, small scale systems [11]; In particular, model reference adaptive 
methods have achieved a great amount of success In producing stable, 
convergent adaptive controllers, and adaptive observers for systems whose 
structure Is known and whose parameters are constant but poorly known 
or slowly time-varying, e.g. : [12]-[25]. Adaptive schemes may be direct . 
I.e., the available control parameters are directly adjusted (adapted) to 
improve the overall system performance, e.g., [25]-[26], or Indirect , i.e., 
the system parameters are identified (based on the assumed system struc¬ 
ture) and the control commands are generated from these parameter esti¬ 
mates as though they were the actual values, e.g., [20], [24], [27], 

The abundance of adaptive control methods Is overwhelming and an 
understanding of the Interrelationships and structural commonality of these 
methods is desperately needed; see, e.g., [28], [29], for some beginnings 
in this direction. Furthermore, the use of such methods on distributed 
parameter or large scale systems, like LSS, is greatly limited by the ROM 
problem - the adaptive scheme must be based on a ROM of the actual system 
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and, hence, the order of the model Is, and must remain, substantially 
lower than the controlled system. In addition. It seems essential that the 
LSS adaptive controller provide a stabilizing control In such highly 
oscillatory systems as LSS . 

This paper develops a framework for LSS adaptive control problems and 
points out generic problems In the use of the most natural direct and 
Indirect adaptive approaches. In other forms, these problems will haunt 
every use of adaptive control on LSS and must be solved before the valu¬ 
able benefits of adaptive control can meet the needs of this new applica¬ 
tion area . A few preliminary attempts at adaptive control for specific 
distributed parameter systems or LSS have been made In [30]-[36]; also, 
for the corresponding parameter Identification problem for distributed 
parameter systems, see [37]. 

2.0 NONADAPTIVE LSS CONTROL 


Following [6], the LSS may be described by the partial differential 
equation: 


m(x) u tt (x,t) + 0 Q u t (x,t) + A q u(x,t) • F(x.t) 


( 2 . 1 ) 


where u{x,t) represents (possibly, a vector of generalized) displacements 
of the structure 0 off Its equilibrium position due to transient dis¬ 
turbances and the applied force distribution F(x,t). The mass distribu¬ 
tion m(x) Is positive and bounded on 3. The Internal restoring forces of 
the structure are represented by Ag u where Ag Is an appropriate differen¬ 
tial operator with domain D(A Q ) defined In a Hilbert space H Q with 
inner product (.,.)g. In most cases, Ag has discrete spectrum, 

l.e., 


A 0 \ “ “k 


( 2 . 2 ) 


where are the mode frequencies of vibration and 
shapes . The damping term Dg u t is generated by an 
differential operator and may represent gyroscopic 
very small (*\.%S critical) natural damping expected 
The applied force distribution Is given by 


^(x) are the mode 
appropriate Ag-bounded 
damping as well as the 
In the LSS. 


F(x.t) ■ Fg(x.t) + Fg(x,t) 


(2.3) 
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where Fq represents external disturbances and Fg represents the control 
forces due to M actuators: 

• M 

F c (x,t) ■ B q f • Y. fc 1 (x) f 1 (t) (2,4) 

1-1 


where b^ are the actuator Influence functions (usually point devices) and 
are the control commands. Observations are produced by P sensors: 

y ■ Cg u ♦ Cg u t (2.5) 

where y^(t) ■ (Cj,u) 0 ♦ (Cj,u t ) Q for 1 < j < P with c^ being the position 
sensor Influence functions and cj the velocity sensor ones (usually point 
devices). 

The state variable form of (2.1) and (2.3M2.5) Is obtained by 
taking 

v(x,t) • Cu(x.t). u t (x.t)]' 

In H s 0(Ag) x H q with energy norm: 

j] v H 2 ■ (mu t , u t ) ♦ (A^ u, A* u) (2.6) 


This produces 


v^ • Av ♦ Bf; v(0) ■ Vg 
y ■ Cv 


(2.7) 


where we consider the disturbance-free case (F Q = 0) and define 5 = [0 Bq] , 
r 0 I 

This distributed parameter system Is 

' A o ■°o > 

very oscillatory In the sense that the semigroup U(t) generated by A has 
very little damping: 


C = ]C 0 C 0 ] and A 


-€t 


U(t)i! < Mq e for t i 0 


( 2 . 8 ) 


where c >_ 0 and small and Mq >, 1. 

The desired performance of the actively controlled LSS greatly effects 
the design of the controller. Many desirable properties of the active 
structure can be obtained with constant feedback gains applied to the sys¬ 
tem state v(x,t); such solutions arise for regulator problems and sta¬ 
bilization (pole placement) problems for LSS. However, the full (infinite 
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dimensional) state v Is never available from a distributed parameter sys 
tern; only the P sensor outputs y are available. 

Imp!ententeble controllers for LSS (and most distributed parameter 
systems) must be based on finite dimensional on-board control computers 
which process the sensor outputs y and produce control commands f; thus, 
a reduced-order model (ROM) of the system (2.7) must be used for the 
controller design. A ROM can be obtained by projecting the system (2.7) 
in H onto an appropriate finite dimensional subspace H R ; the projections 
P (onto H, ; and Q (onto the residual subspace) are usually, but not 
always, orthogonal. Let v R « Pv and v R * Qv and, from (2.7), we obtain: 


*N " A N V N + A NR V R + B N f 
*R * A RN v N + A R V R + B R f 
* ’ C N V N + C R V R 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


where A N * PAP, A NR * PAQ, * PB, etc. The terms B R f and C R v R are 
called control and observation spillover ; the terms A RR v R are called 
model error . The ROM for this system is given by (2.9) and (2.11) with 
A nr . 0 and C R - 0: 


'N 


** V N * 8 N f 


y • C., v H i v M (0) • Pv„ 
resid 

with total energy }{ v[| given by: 


( 2 . 12 ) 


The ROM state v R and the residual state v R form the true system state v 

I „ !! 2 


j iMi 2 * 11 v N Sl 2 + Uv R lj 2 (if projection is orthogonal) 

All implementable controller designs based on any ROM must be evaluated 
in closed-loop with the actual LSS (2.7),and it is in this evaluation that 
the effects of model error and spillover due to the residuals become 
apparent. 

If the actual mode shapes 4^ are known, the modal ROM is a sensible 
choice: 


H N *^ sp *1 ••• V 

and the model error terms A^ R and A RN become zero. Of course, any collec¬ 
tion of modes could be used; usually, the most easily excited or critical 
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ones will be chosen. However, In many cases, the partial differential 
operator A Is too complex to provide closed-form mode shapes. Instead 
finite element approximations of the mode shapes are computed (e.g., 
via NASTRAN) i and these approximate mode shapes can be used to form the 
ROM for controller design; note that some model error is present when 
these approximations are used. Henceforth, we will assume the actual 
mode shapes are available to simplify the discussion but much of our 
analysis remains valid for approximate mode shapes and other types of 
ROM. 


Modern modal control (MMC) for LSS, as developed In [38], uses the 
modal (or approximate modal) ROM and develops a controller consisting 
of a state estimator based on the ROM and a constant gain control law: 


*N ' A N ? N + B H f + S'™* 
y * Vj^, v^(o) * o 


(2.13) 

(2.14) 


This controller design requires the ROM (A^, 8^, C^) to be controllable 
and observable for the calculation of control and estimator gains G^, 
K^. These conditions, in modal terms, provide insight into the number 
and location of actuates and sensors. From [38], (A^, B^, C N ) is con¬ 
trollable and observable, when position sensors are used (Cg * 0), if 
and only if 


(1) min(P,M) max mode frequency multiplicity in the ROM 

(2) each sub-block of and associated with a mode frequency ^ of 
multiplicity must have rank at least equal to c^,. 

Similar results hold for other types of sensors, e.g., velocity, acceler¬ 
ation, or mixtures of types [38], [391. These results are easy to in¬ 
terpret in terms of the mode shapes, e.g., if no repeated frequencies 
exist, then the above result says that a single actuator and sensor, 
not necessarily collocated, will do the job as long as neither is located 
at any of the ROM mode shape zeros. Since LSS have many symmetries and 
rigid body modes, it is not often that a LSS control problem will have a 
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controllable observable ROM with only one pair of devices; this has con¬ 
sequences for the adaptive control problem to be discussed later. 

Let the estimator error ^ 5 V N - v N be defined and, from eqs. (2.9) 
-(2.11) and (2.13)-(2.14), obtain 

*N " < A N + B N G fP V N + B N S N e N ^ 2 ‘ 15 ^ 

*N * ^ A N ' *N C N^ e M + h C R V R ^ 2 ‘ 16 ^ 

*R " B R G N v N + B R G N e N + A R V R (2.17) 

This shows the effect of spillover on the closed-loop system: even 
though A n + B n G n , A n - C N , and A R are stable, the closed-loop system 
need not be stable. When either control or observation spillover is 
absent (B R * 0 or C R * 0), then stability is assured; otherwise, spill¬ 
over causes pole-shifting and can induce instabilities [38], [40]. 

Bounds on the destabilizing effect of observation and control spillover 
were produced in [38] and can be extended to the case where some model 
error and small nonlinearities are present [41]. Such bounds give an 
indication of how much spillover the closed-loop system can tolerate. 

A variety of methods have been suggested to reduce spillover [6]. 

One obvious way would be to prefilter the sensor outputs, with a bandpass 
filter, to substantially reduce observation spillover. This alleviates 
the worst pola-shifting problem; bounds cn system performance with con¬ 
trol spillover alone can be found in [42]. Note that the post-filter 
of the controller outputs could do the same job by reducing control 
spillover; this interchange of filter and controller is possible due to 
linearity and time-invariance. The trade-off with this means of reducing 
spillover is that the prefilter introduces phase distortion which can 
have a destabilizing effect of its own. Therefore, a very high order 
filter may be required to keep the phase distortion acceptable; phase- 
locked-loop quadrature filters may be another solution [43]. E ven in 
nonadaptive LSS control, the spillover and model error problem is a funda 
mental one . 

Finally, we should note in this section that digitally implemented 
controls would be based on discrete-time versions of the distributed 
parameter system (2.7). One such version is obtained by using a uniform 
time step At: 
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v(k + 1) « $v(k) + E f fc 
y(k) - Cv(k) 


(2.18) 


where ♦ = U(At) and E = E Q B 


r At 

/ U(t) dr B and the control command is 
J 0 


a constant f k over the interval (k 1) At <, t < k At. Other versions of 
this could be obtained with nonuniform time steps. When the ROM procedure 
of projecting onto the subspace H N is used, we obtain the discrete-time 
versions of (2.9)-(2.11): 


v N (k‘ 4 1) 

■ «N * ♦«* V k > + E N f k 

(2.19) 

v R (k + 1) 

* *RN »N< k > * *R v R (k > + E » f k 

(2.20) 

y(k) 

* C N V k ' * C R V k - 

(2.21) 

H R is the modal subspace, the above become: 


v N (k + 1) 

" $ N + E 0 B N f k 

(2.22) 

v R (k + 1) 

• * R V R 00 + e 0 b r f k 

(2.23) 

yOO 

* C N V k ^' + C R V R^ 

(2.24) 


where (2.22) is the same as that obtained by directly discretizing the 
ROM in (2.12); if the exact mode shapes are not available, these two 
discretizations may yield different results. In addition, the sampling 
process can alias residual modes and increase observation spillover and 
the zero-order hold process can spread-out the control command signal 
spectrum and, hence, increase control spillover by increasing the energy 
in the residual mode spectrum; this has been observed and investigated 
in [44]. Therefore, the time discretization is a very important factor 
in the design of implementable LSS controllers. 


3.0 TOWARD ADAPTIVE CONTROL OF LSS 


In order to design MMC, or other controllers, for LSS as proposed in 
the previous section, we must have knowledge of the ROM parameters (A^, 
®N’ C N^* " rhese Parameters are obtained from modal data; they are the 
mode frequencies for A^ and the mode shapes ac actuator and sensor loca¬ 
tions for and C^, respectively. This data is required for three 


reasons: 






(1) to determine controllability and observability of the ROM and, 

hence, to help locate control devices effectively; 

(2) to design control and estimator gains; 

(3) to use In the state estimator's internal model. 

However, we have noted in Sec. 1.0 the sources of error for this data; 
consequently, a need arises for an adaptive version of the MMC of Sec. 

2 . 0 . 

The most logical and reasonable procedure to obtain adaptive con¬ 
trollers for a LSS seems to be the following: 

Procedure for Adaptive LSS Control 

(a) choose a “nice" reduced-order model (ROM); a modal ROM 
would be the obvious choice; 

(b) use your "favorite" lumped parameter adaptive control 
scheme; 

(c) design the adaptive controller as though the ROM were 
the actual LSS to be controlled, i.e., ignore the un¬ 
modeled residual part of the structure; 

(d) use this adaptive controller in closed-loop with the 
actual LSS and hope for the best. 

There is nothing wrong with following this procedure as a best first guess 
- in a way, there is little else that one can do to produce an implement- 
able adaptive LSS controller. 

In some cases, spillover is sufficiently small or enough other mathe¬ 
matical structure is present in the system, e.g., a high level of damping 
in the distributed parameter system [31] or low level of performance re¬ 
quired from the controller (increased damping via direct velocity feedback) 
[32], to allow the adaptive controller to operate. However, these situa¬ 
tions are rare with LSS and one would not like to count on the "generosity 
of nature" or the temporary suspension of Murphy's Law as part of the 
above design orocedure. Consequently, we would add the following items to 
that procedure: 

(e) analyze computer simulations of higher-order models of the 
LSS in closed-loop with the adaptive controller base-' on 
the lower-order ROM (e.g., simulate more modes than you 
plan to control); 















(f) investigate the specific mathematical mechanisms whereby 
the residual (unmodeled) part of the LSS couples Into a 
given adaptive control scheme (e.g., find out where and how 
spillover affects the adaptive controller); 

(g) obtain mathematical results on the amount of spillover 
and/or model error that can be tolerated In the closed- 

loop system and still provide adequate adaptation and control; 

(h) develop spillover and model error compensation schemes to 
augment the adaptive controller when the residuals cannot 
be tolerated (e.g., when the conditions of (g) are not sat- 
isfied); 

(i) recheck (g) with this compensation in the closed-loop system. 

We believe that, within the basic framework of LSS control as de¬ 
scribed in Sec. 2.0, this Augmented Procedure (a)-(i) will go a long way 
toward revealing the problems of adaptive LSS control (and indeed, most 
other adaptive control situations where control must be based on a ROM) 
and will help to focus needed attention on these crucial issues. For 
example, although Step (e) would be done most likely at some point in the 
system development phase of any project as complex as the construction and 
operation of a LSS, often it is done much too late and the design is 
"set in concrete (or in th-'s case, graphite-epoxy)"; the other steps 
(f)-(i) may not be done at all. Yet, ignoring the effects of the resid¬ 
ual unmodeled LSS can produce some disasterous behavior in the adaptive 
controller; this was pointed out quite clearly in the LSS example in [36]. 

When an adaptive scheme is applied to such a lightly damped, oscil¬ 
latory system as a LSS, the stability of the closed-loop system during 
adaptation is a necessity ; therefore, we do not view convergence and 
(global) stability results as luxuries for adaptive LSS control and shall 
only consider appropriate those lumped parameter adaptive schemes for 
which such results are available. However, even a globally stable adaptive 
scheme may prove to be unstable when it is used in closed-loop with the 
actual LSS instead of the ROM on which it was based . This is not a failure 
of the adaptive scheme; it is a failure to satisfy the mathematical hypoth¬ 
esis of the stability result associated with the scheme. 








daptive sch 
daptlve LSS 


use an adaptive observer and operate In continuous time, and the discrete 


time, direct or Indirect, schemes which are based on an Autoregressive 


Moving Average (ARMA) model of the controlled plant, e. 


These approaches represent a good cross-section of available lumped 
parameter adaptive control schemes which have been shown to possess the 
desired stability properties. We emphasize that the point of this sec¬ 
tion Is not to criticize or slander these schemes; rather, we mean to 
point out where the hypotheses of their stability results are violated, 
and must be modified, when we attempt to use them on LSS. We feel that 
consideration of these approaches within the context of the Augmented 
Procedure for Adaptive LSS Control (a)-(1) will illustrate the generic 
difficulties In the application of existing, well-behaved, lumped param¬ 
eter adaptive control schemes to LSS. 


3.1 Multivariable Systems Converted to Scalar Systems 

The results of many stable adaptive schemes, e.g., [13]-[14], [23]-[24], 
are limited to a single actuator and/or a single sensor; yet, we have seen 
In Sec. 2.0 that most LSS co.'.trcl problems will involve multiple actuators 
and sensors. One way to deal with this (although, admittedly it has its 
drawbacks) is to convert the controllable observable LSS problem via out¬ 
put feedback into one that is controllable and observable from a single 
actuator and/or single sensor; this can be done with almost any output 
feedback gains [45]-[47]. These gains would have to be based on the best 
available calculated ROM data and the designer must hope that they will 
continue to do their job during adaptation. 

The output feedback modifies the original system (2.9)-(2.11) to be- 


( A N + B N H N V V N + ( A NR + B N H N C R' V R + b N f 
( a dm + b R h N C n ) v n + (A 0 + S B H n C 0 ) v d + b D f 


R N R 


C N V N + C R V R 


where is the output feedback gain matrix, b^, b R , c^, c R are vectors, 
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and f,y have been renamed. Let H N C N be X N , etc. and we have that 

the new ROM (I N , b N , cj) Is a controllable, observable single-input, 
single-output system and (3.1) and (3.2) become: 

V N * *NR V R + b N f * 3,4 ) 

^R “ *RN V N + *R V R + b R f ^ 3 ‘ 5 ^ 


3.2 Indirect Adaptive Controller Design 


We apply the design of the adaptive controller In [24] directly to 
the ROM consisting of (3.4) and (3.3) with the assumption, for now, that 

A„ r - 0 and cj • 0. 

The control law Is given by 

f(t) ■ gj v N (t) + f c (t) (3.6) 

where g N is a constant gain vector, f is a "sufficiently rich" external 
signal (more about this later), and v N is derived from the following 
adaptive observer (or state estimator): 


v N * F v N + gy + hf;'v N (0) - V Q 


(3.7) 


where F is in arbitrary, stable matrix and g, h are unknown parameter 
vectors. The appropriate matching conditions are: 


( F + G * c 5 ■ 

I'D * V°> 


(3.8) 


where g*, h* are constant. 

Let p* = [g*^ h* T v N (0) T ] T and note that 

t>* • F p*; p*(0) * p Q (3.9) 

where F = diag [0 0 F] and, when g * g*, h * h*, we have from (3.8) and 


(3.7): 


v N (t) * [M(t) I N ] p*(t) (3.10) 

where M(t) 5 j 1 e F ^ [I N y(-r) 1^ f(t)] dr and we have used the fact 
that h, f are scalars. 






I 


Mow (3.7) can be rewritten as 
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v N (t) - (M(t) p(t) 

A(t) - F M(t) + [I N y(t) I N f(t)] 

M(0) ■ 0 


(3.11) 

(3.12) 

(3.13) 


where M(t) Is as defined In (3.10),and we have yet to produce the adaptive 
law to generate p(t). This adaptive law is given by the following: 


A(t) « F p(t) - a(t) [R(t) p(t) + r(t)] 


(3.14) 


where p(0) • Pg is arbitrary and 

R(t) - -q R(t)- F 7 R(t) - R(t) F + [M(t) I N ] 7 c N cj [M(t) I N ] 

' (3.15) 

R(0) ■ 0 


f(t) ■ -q r(t) - F 7 r(t) - [M(t) I N ] 7 c N y(t) 


r(0) » 0 


(3.16) 


where a(t) Is the adaptive gain and the constant q is chosen to exceed 
twice the absolute value of the real parts of the eigenvalues of F. The 
adaptive gain is chosen so that 


a(t) « Y + i ! U ( t )I| 


(3.17) 


u(t) ■ -a u(t) + (N/2) 1 * (;y(t); + |f(t)|) 

with X positive and F + F 7 <_ 2X I N . 


(3.18) 


3.3 Convergence Results: What Goes Wrong? 

All of the above is exactly as stated in [24] where It is also shown 
in Appendix I and II that 


u(t) > ' jM(t)i| 


(3.13) 


and, with f(t) sufficiently rich in frequencies, there is a t^ such that 
R(t) > o I N > 0 for all t > ^ (3.20) 
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In addition, It Is shown In Appendix III that 
R(t) p*(t) + r(t) s 0 (3.21) 

This Is very crucial to the stability results of [24] and It Is here that 
observation spillover (l.e., the fact that we are using a ROM) appears - 
(3.21) Is not valid when c R t 0; however, 

R(t) p*(t) ♦ r(t) • A R (t) (3.22) 

where A R (t) s J e’^ + ^ ( € " T ) [m( t ) I n ] T c n cj v r (t) dx 


Let e N (t) 5 v N (t) - v N (t) and Ap(t) 2 p(t) - p*(t); then 


e N (t) - [M(t) I N ] Ap(t) 

(3.23) 

and 


Ap(t) • [F - o(t) R(t)] Ap(t) - a(t) A R (t) 

(3.24) 

Consider V(t) = Ap(t) T Ap(t) and we obtain: 


V(t) < -y p V(t) - 2y A R (t) Ap(t) 

(3.25) 


where V(0) - ApJ Ap Q and a Pq ■ p Q - p*. 

This follows [24] except for the additional term In (3.25); also 
note that 

l!e N (t)!! 2 < V(t) [1 ♦ u(t)] (3.26) 

Let iv N (t) = v N (t) - vjj(t) and Av R (t) = v R (t) - vjj(t) where vj^ and v£ 
represent the ideal states of (3.4)-(3.5) when the parameters are 
exactly known: 

" <*N * b N ’!> V N * J m V R + b N f C (3 ' 27) 

*R " (r RN * b N 9|(> v N + *R V R + b R f c (3 ' 28) 

with vjj(0) » v R (0) and vjj(0) ■ v R (0). 

When the implementable adaptive control law (3.6) Is used, we obtain: 

* (Ajj + b^ o*j) -^v^ + + ^ 9 m (3.29) 

Av R “ ^RN + b N 9 N^ Av R + *R AV R + b R 9 N e N ^ 3,30) 

and we have the following result: 
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THEOREM 3.1 : Assume 

(D c R ■ 0 

(2} X N + t> N gj stable 

(3) X R stable 

( 4 ) J m ■ 0 

(5) f £ (t) sufficiently rich (l.e., It has at least 3N/2 
distinct frequencies) 

2 

Then there is a 5 > 0 such that for all IUpqH <. 6: 

(a) « N (t) Is bounded and (eventually) vanishes with an 
arbitrary exponential rate 

(b) lira Av N (t) ■ 0 

t— " 

(c) 11m Av 0 (t) ■ 0. 
t— K 


Therefore, even though the closed-loop system with the adaptive 
controller Is highly nonlinear, it Is stable while the adaptation 
Is taking place. In particular, (1), (3), and (4) are satisfied 
If there Is no observation spillover (C R ■ 0) and no model error 
(A^ r * 0) and some damping In the residuals (l.e., A R is stable). 


Proof : This result follows from the results of [24] because A R (t) = 0 In 
(3.22) when c R ■ 0. The stability of (3.29)-(3.30) Is determined by 
A N ♦ b N gj and A R when ^ NR ■ 0. 


Also, if t R • 0, then A^ R 


A A * A 

*NR’ *R V 


and c R ■ 0 since it is a row vector of C R . 


Note that the stability of the system with adaptive control Is 
determined by that of (3.27)-(3.28) - the Ideal case where the parameters 
are known and the external signal f is applled . This is natural, since 
adaptation cannot take place without f present; however, after adapta¬ 
tion, we would most likely want to turn off f . In addition, we could 
choose the f signal so as not to excite the residual frequencies whenever 
sufficient spectral separation is present. Still, we would need to turn 
on f now and then, in order to "tune-up" the controller. 


3.4 Spillover Compensation for the Indirect Adaptive Controller 

The above result is merely a slight extension of the results of [24] 
to a special case of the adaptive controller based on an ROM Instead of 
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the full system. However, It does suggest that some form of compensation 
should be used to eliminate the observation spillover. Such compensation 
must be essentially Independent of the parameters of the ROM; yet, most 
methods of spillover reduction require knowledge of the ROM (and some 
residual) parameters. 

One approach to spillover compensation already suggested In Sec. 2.0 
Is prefllterlng the sensor outputs to remove or greatly reduce the term 
C R v R (t). Such prefllterlng can be achieved with low-pass or band-pass 
filters when the ROM frequencies are known and separated from the residual 
frequencies. However, the modal frequency data Is part of the poorly 
known parameter data. 

In an attempt to resolve this predicament, we could try using phase- 
locked loop (PLL)-based filters with the center frequency of each loop 
tuned to the best approximation available of the corresponding ROM data. 

The Pll will adapt Itself until It tracks the actual mode frequency, and 
after "lock-on," It will behave as a narrow band-pass,11near filter which 
tunes out the observation spillover from the other frequencies [43]. Of 
course, sufficient spectral separation must be present, the calculated 
values of the ROM modal frequencies must be sufficiently good, the dis¬ 
tortion Introduced by the filter must be sufficiently small, and the 
adaptive controller must not shift the poles around too much. Thus, the 
PLL prefllter Is' not a panacea/ But, It might work to reduce spillover 
and. If It does. It has the added advantage that Its output could also 
reveal better estimates of the modal frequencies; this would take some 
of the load off the adaptive observer. If the adaptation mechanism causes 
too much pole shifting, the ROM frequencies could be excited via f and 
Identified In open-loop by the PLL filters before the adaptive controller 
Is turned on. Note that modal frequency data Is usually better known, 
via computer approximation, than modal shape data; hence, this approach 
might not be unreasonable. Another possible, but untried, approach to 
spillover compensation might be an adaptive version of the orthogonal 
filter In [48], 

Note that some prefllterlng (and postfllterlng) always takes place 
due to the bandwidth limitations of the sensors (and actuators). Whether 
this can be used to advantage In adaptive LSS control Is still a matter 
for speculation, e.g., [33]. 
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A problem that arises with the use of the results In [24] and their 
modification to LSS Is that the constant feedback gains gl must be cal- 
culated, In advance, to stabilize A^ ♦ b N g^. It would be better If 
these gains were adapted along with the parameters In the observer. Of 
course, after adaptation has taken place, they could be recalculated 
from the "tuned-up" parameters, but. In some cases, the adaptation phase 
Is never over, e.g., slowly-varying parameters. Other approaches could 
be used for adaptive pole-placement, e.g., [20], [21], but these also 
have their limitations. 


3.5 ARMA-Gettln' 

Many discrete-time adaptive control schemes depend on an Auto- 
Regressive Moving Average (ARMA) representation of the plant In discrete¬ 
time, e.g., [25]-[27]: 

N N 

y(k + N) - Y. a r y(k + r - 1) ♦ £ S p f(k ♦ r - 1) (3.31) 

r«l r«l 

for some H and appropriate matrices a . S„. What the ARMA says Is that, 

r r 

after N time steps, the present output Is only related to the past N out¬ 
puts and Inputs. Existence of an ARMA Is directly related to the finite 
dimensionality of the plant (N Is usually that dimension) and Is obtained 
using the Cayley-Hamllton theorem for matrices. For LSS, only a quasl- 
ARMA can exist; these were considered in detail In [4S]. From the 
Appendix of [49], we obtain the quasi-ARMA for the LSS (2.13) or (2.19)- 
( 2 . 21 ): 

N 

y(k + r - 1) + £ r r E n f(k + r - 1) ♦ R(k) (3.32) 
r*i 

H 

N) ♦ £ * R (k ♦ r - 1) 
r*l 


y(k + N) » £ a r 
r-1 


R(k) - C R v R (k + 


A r " r r ^NR ’ a r 


where 


and ? r Is defined In the Appendix of [49]. Since R(k) » 0 when C R ■ 0 
and & r ■ 0, we have the following result: 

THEOREM 3.2: When the observation spillover (C R ) end the 
model error term (4^) are both zero, the quasi-ARMA (3.32) 

Is a true ARMA for the LSS (2.18),and any stable adaptive 
scheme based on this ARMA will be globally stable when used 
In closed-loop with the actual LSS (2.18). 

When the rather stringent hypothesis of Theo. 3.2 Is not satisfied 
(as It may not be In practice), any adaptive LSS control scheme based on 
the quasl-ARMA (3.31) must Ignore R(k) In order to be Implamentable. How¬ 
ever, R(k) Is the term where the residual effects - spillover and model 
error - enter the scheme and can cause Instability. Again, as In Sec. 

3.4, prefllterlng or other compensation might be tried In an attempt to 
reduce or eliminate this term. 


4.0 CONCLUSIONS 

In an attempt to point out the crucial Issues and generic problems 
associated with adaptive control of large aerospace structures (LSS), we 
have reviewed the framework for nonadaptlve LSS control (Sec. 2.0) and, 
within this framework, have proposed a general procedure, based on 
reduced-order models (ROM) of the LSS, for obtaining and assessing the 
problems of adaptive LSS controllers (3.0). In addition, we have taken 
a closer look at the use of certain well-known, lumped parameter, stable 
adaptive control schemes In chis procedure. Taking these schemes as 
representative of the basic Ideas present In all lumped parameter adap¬ 
tive control approaches, we have obtained corresponding LSS adaptive 
controllers and found the following generic problems associated with 
adaptive LSS control : 

(1) LSS are distributed parameter or large scale systems; 
therefore, the plant dimension Is always larger than 
the dimension of the adaptive controller, which must bs 
based on a -.OM; 
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(2) LSS control must often be done with more than on: actua¬ 
tor and sensor; conversion of multivariable to scalar 
systems via output feedback introduces problems of its 
own (e.g., stability of the residuals); 

(3) LSS control problems are often non-minimum phase due to 
noncoilocated actuators and sensors; 

(4) Interaction of the residuals with the adaptive control¬ 
ler may negate the stabilizing properties of the con¬ 
troller due to observation spillover; this Interaction 
Is much worse due to the nonlinear nature of adaptive 
control; 

(5) Methods of spillover compensation for LSS often require 
knowledge of the ROM (and some residual) parameters - 
the very data that are poorly known; 

(6) The adaptation mechanism may shift the closed-loop fre¬ 
quencies around; this counteracts the benefits of any 
prefiltering unless sufficient spectral separation is 
maintained; 

(7) Indirect adaptive controllers need sufficient excitation 
from an external signal f ; however, this signal may 
substantially excite the residuals. 

(8) Discrete-time adaptive controllers can only be based 
on quasi-ARMA rather than strict ARMA representations 
of the LSS; this may negate the stability properties of 
such a controller. 

Stable adaptation is essential for such highly oscillatory systems 
as LSS, yet our preliminary stability results, Theos. 3.1 and 3.2, both 
require that observation spillover be somehow completely eliminated 
before it reaches the adaptive control logic; certainly, this is not 
an easy thing to do in general! Perhaps, global stability is too much 
to ask for LSS adaptive control because observation spillover will 
always be present to seme degree in LSS control. However, it seems rea¬ 
sonable to hope for the development of spillover bounds to give some 
idea of regions of stability for the successful operation of LSS adaptive 
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control. Some comparison should be made between stable adaptive con¬ 
trollers based on ROM and stable, robust control schemes, e.g., [50], 
[51]. 

In closing, m would like to say that it is not our Intent to present 
a gloomy picture for the application of adaptive control to LSS . In 
fact, the need for adaptive control in LSS is already becoming quite 
clear, and recognition of this need comes, for a change, at an appropri¬ 
ate time - before any LSS have been built and put into space. However, 
the development of adaptive control for LSS will not take place overnight 
and will not be done by one or two people. Consequently, what we have 
tried to stress here for the interested researcher are some of the funda¬ 
mental problems that arise and the basic steps which need to be taken 
toward the goal of successful adaptive control of LSS. In the long run, 
we have high hopes for the success of this endeavor and we expect that 
adaptive control theory will profit by its association with large aero¬ 
space structures, as well. 
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Tae quantitative effects of spillover on reduced 
crder adaptive coscrol are examined via a :tsc 
exams It. 

: H.TRODVCTION 

A major assumption in :.ia various rectnc pruofs 
(1) - !5) of luaped-paramettr system (IPS) adaptive 
control algorithms is the adequacy of the order of 
i.*.t ootorciitr. The mooel of the plant used ir. 
s ; a c i i yo. ng the orcer of the controller is required 
:: taro- or exceed shat o: rhe plant. For example 
, fcr the model reference adaptive control of the 


v(t) « ;■ a.vC'K-i) - b.u(k-j) (1) 

i-1 ‘ : 

ir. the absence nf accurate a priori specification of 
the plant parameters a^ arc an adjustable control 
lav of the form 

uU-1) • c, (k-l) r(k-l) * l c . (k-l)u(k-j) 

2 - 

- d 4 (k-l)y(k-i) (!) 


with a > p and n » q these similarities disappear 
along with the ability for the output error s - y to 
converge to zero. In this case it is expected that 
the results of indirect anc direct adaptive control 
will be quite eiscinguishacie, which counters she 
current trend of exclaiming their equivalence (7)-(10). 

The use of a reduced-order adaptive controller, 
such as (1), can arisa from a desire to reduce control 
complexity vlth a coosensurace reduction ir. controller 
sensitivity (_11) or from necessity in an attaept to 
use a finite dimensional controller or. ar. infinite 
dimensional system (12). A recant interest (_13) - (15) 
ir. develooing adaptive controllers for flexible struc¬ 
tures by aoplication of existing IPS adaptive control- 
la- strategies oo truncated expansion descriptions of 
the distributed parameter system (DPS) behavior of 
flexible structures has led to rvo qualitative descrip¬ 
tions of the misbehavior of reduced-order adaptive 
controllers. These interpretations of tht additional 
cifficulties facing reduced-order adaptive controllers, 
which arc bypassed by exact-order adaptive controllers, 
are suimnarized in the next section. The last section 
formulates a test problem, which Initializes attempts 
to quantify these qualitative msigr.es. The most 
current results of this investigation will be compiled 
in the conference presentation. 

II R£DUCED-ORDIR IFTiCTS IN .CAPTT.~ C0NTRC1 


is required cc insure the asymotocic convergence of y 
to the output s of the model 


Svk, ■ gr(k-i) ♦ _ f.s(k-i). (3) 

i«l * 

Not* cr.ar im (1) the susation limits m and n must 
correspond exactly to the moving-average and auto- 
ragression lengths in (1), respectively. As shown in 
(5), the problem of (1) - (3) is solvable via a 
variety of Indirect and dirset adaptive control tech¬ 
niques that prove quite similar in this exact order 
case. However, in the reduced order case where the 
centre! input is provided by 

P - 

u(k-l) - c, (k-l)r(k-l) * l c (k-l)u(k-j) 
j-2 J 
c 

- " 3. (k-l)y(k-i) (4) 

_ i »l 1 
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Though numerous strategies have been espoused 
( 16 ) for extracting a low order model from a too com¬ 
plex system description, none currently seem fully 
anclicable to the real-time, recursive requirements of 
on-line adaptive control algorithms. These reduced 
order modeling techniques seem to fall into rue broac 
categories: (i) extract cnose "modes" (or component- 
subsystems) from the full system description that are 
most influential la tha performance of the model in 
its subsequent use. This strategy is followed, for 
example, in (17) and (18). (il) Parameterize the 
reduced structure to provide the best prediction of the 
desired output. This latter strategy will produce a 
modal the "modes" of which need not correspond to any 
of chose of the full system as noted in (14). This 
latter strategy includes the model reference approach 
prominent in adaptive systems (19). 

The presentation in ( 14 ) emphasizes the difficulty 
of reduesd-order adaptive control by contrasting these 
two model reduction strategies. Consider the approxi¬ 
mation of ^ 

d(9,k) ■ l H (k)9 (6) , (5) 

J-l J J 


ORIGINAL PAGE IS 
OP POOR QUALITY 












vh«r« d can serve as the deflection of a flexible 
structure, w« the nodal amplitudes, and tj th« nod* 
shapaa for dH or with t< ! 1 u cha output of a 
dacouplad (partial fraction axpanaion) daacripcion of 
a IPS with Uj aa the individual "nodal" outputa (20). 
Assuming that the nodal asplltudaa obey a diffarence 
aquation aueh aa 

n 

V,(k) - 1 [a 4 ,W (k-i) ♦ b,.F,(k-i)], (6) 

J i.j ji J J* J 

vhara cha Fj arc the nodal inputa and the ij< and bj* 
are cha unknown paranatera of incaraat, the first 
strategy for adaptive identification of (6) require* 
extraction of the appropriate Wj fros (5) foraing 

Mof'[!,»] 

d R (e,k) - l Mjfkja^e), (7) 

which aa noted in (12) , (1£) and (_2£) remains an open 
problen without knowledge o: N and the ty the 
second strategy would estimate the Uj such that 

Mof[l,S]. 

d(e,k) • l W_(k)o.(j) (8) 

j 

closely approximated d in (5) in, «.*., a least- 
squares sense. The estimated nodal aaplitudea 
need not provide useful values for the estimation of 
the a-iJ and b<< in (6) with W. repiacad by Wj as noted 
in (l £). In IS) c.na chosen at each tine instant 
to closely natch a to 5 say not even obey an under¬ 
lying difference ecuation such as (6). The problen 
with closing the acaptive control loop is tnen inten¬ 
sified by the difficulty in idencifying (6) and than 
in feeding back the w< instead of the for modal 
control. 

The difference in various identification tech¬ 
niques in the reduced-order IPS case can be emphasized 
by a simple exaaplt froc (21). Consider identifying 
the stable, second-oruer plant 

y(k) • 1.131-y(k-i; - 0.25y(k-2: * Q.05u(k-1) - C.-uik-2) 

(9) 

using a first order modal and three conson adaptive 
identification schemes, tr.e equation error formulatec 
IMS algorithm (_22), the gradient based output error 
formulated Steam's algorithm (SA) (21). and the 
stability based output error formulated SHARP (2<0. 

The algorltnms art 


UIS: 

y(k) • 

*(k)y(k-l) * b(k)u(k-l) 

(10) 


a(k+l) 

■ a(k) + u(y(k) - y(k) ]y(k-l) 

(11) 


b(k+l) 

“ b(k) - o[y(k) - v(k)]u(k-l) 

(12) 

SA: 

y(k) - 

a(k)y(k-l) * b(k)u(k-l) 

(13) 


a(k+l) 

- a(k) + u(y(k) - y(k)]X(k) 

(1 i) 


l(k) • 

y(k-l) + a(k)X(k-l) 

(13) 


b(k*l) 

* b(k) - b[v(k) - y(k))y(k) 

(16) 


y(k) - 

u(k-l) + a(k)y(k-l) 

(17) 


SHARP: y(k) • 

a(k)y(k-l) * b(k)u(k-l) 

(18) 

a(k+l) 

■ *(k) ♦ uv(k)y(k-l) 

(19) 

bOc+1) 

• b(k) * ov(k)u(k-l) 

(20) 

v(k) ■ 

y(k) - y(k) * X[y(k-1) - y(k-l)] 

. (21) 


If the y in (9) were generated by a first order modal 
as are (10), (12) and (18), than each of the three 
schemes (if u, s, and 1 ware ehoeen to asaure stabi¬ 
lity) converge to the unique correct velues given 
sufficiently rich u. However for the reduced-order 
identification tasks of (9) - (21) each scheme yields 
different mean convergence points as shown Id Figure 
1. The Input u in this case is white and xaro-maan. 
The contours in Figure 1 represent constant normalized 
mean square error 

t(e 2 (k)} - S([y(k) - ay(k-l) - bu(k-l)] 2 ) (22) 

for various choices of a and b. Soca that tha 
graditne-based SA is suscaptibla to local cenvergenca 
and for this problem the SA mean eonvergance point is 
dependant on the initial guesaaa for a ar.c b. Though 
if any of these schemes is used for simultaneous 
identification and control the input would r.o longer 
be white, clearly tach would have a quit* difftrtnc 
control oerformanct. 1: becomes unclear as to which 
would adaotivaly provide the best reauced-orcer modal. 

An alternate interpretation of tha "dual" ill- 
effects o: reduced-order adaptive control has emerged, 
also from the flexible structure toncrol problem (12? 
(15). “or this viewpoint consider separation of a 
single-input, single-output system state-space de¬ 
scription of the reduced-oroer model and tha remainder 
(or residual) as 


"V^i" 

A; 

a "s' 

' i + i 'u(k) 

(23) 

i x^(k+l) j 


sjLv k) j ^rJ 


y ( k > — 

V 

"s U) l 

(2-) 


_V k) _ ■ 

Prom (23) - (2d) the "spillover" of the residual modes 
into the observation of y(k) via eqxq(k) and the spill¬ 
over of the control designed for the reducee model 
into r.hc residual modes xq(k+l) via bqu(k) arc clearly 
displayed. Also the possible coupling of the reduced 
and modal modes via A;;q and Aqx is immediately 
eeparenc. Shown m (12) *nc ( 15 ) i* tne predictable 
fact that if As* - 0, and cq ■ 0, then, assuming the 
residual modes are noc unstable, any full-order-stabl* 
adaptive controller identifying A»|, bs, end c«g 
explicitly or implicitly from y and u alone would 
retain its stability. If Ajfq ■ 0, but eq and bq ar* 
nonzero then the degradation of reducad-order adaptive 
control ha* two sources. Unmodeled eomponancs in / 
via cqxq will generate an error that is indistinguish¬ 
able from parameter error causing adaptation. The 
application of control u to xq via bqu will cause 
further excitation of the residual modes further per¬ 
turbing the reduced modeling of y. Not only will the 
parameter estimates become incorrect but also the 
state estimates from an adaptive observer will be 
incorrect leading to an unpradlctabla "controlled” 
response. 











(30) 


# These two probleas of inappropriate pirutuc 
and state estimation are the.sice ones noted in the 
earlier interpretation with Vj rather chan Vj 
txtraction. The ability to obtain dp corraaponda to 
an effective itroing of cr. Sinca auch ability doaa 
not appaar forthcoming in ganaral, tha affacca of 
apillovar on reduced-order adaptiva control via tha 
f taesne stratagy notad at tha beginning of thia 
” taction Boat ba avaluatad. a taat axaapla follows 
wnich it daaignad to provida intight into chit 
problaa. 

III TEST EXAMPLE 


Conaldar controlling tha second-order plant with 
f tranafar function 



and tharafora diffaranea aquation daacrlption 
^ y(k) • (a.*a,)y(k-i) - (a,a,)y(k-2) 

- (b 1 *b : t)u(k- ; ) - (b 1 a,-cb,a,)u(k-2) (26) 


V k) ‘ bOt+1) 


•iilU 


d(k) - f 1 (k)r(k) + f 2 (k)y(k). (31) 

Tha connon adaptiva concrollar froa (£) 

f (k) - f.(k-l) * or(k-l)[dr(k-l) ♦ cy(k-l) - y(k)] 

. - (32) 

f,(k) • f-(k-l) yy(k-l)[dr(k-l) ♦ ey(k-l) - y(k)) 

i (33) 

with j and p appropriataly coaacralnad will para- 
aatarlxa (31) for tha diract caaa. 

Pralialnary raaulta ara aablguoua. Mona of tha 
alaulacad taatt daarly aupport tha unlvaraal useful- 
naaa of naglacting high frequency aodat, tha tupari- 
orlcy of aodal rajaction or actaaptad "full" adaptive 
tracking, or tha clear benefit of indirect or diract 
adaptiva tachniquat. In other wordt, one approach it 
batter in one axaapla than another, which require! 
further invattigaticn of raducad-ordar adaptive 
control. The mixed evidence will ba interpreted in 
tha conference pratantation. 


under tha attuaption of a firtt order aodal. Using 
tha tenaration technique of (23) - (21), (26) can oa 
rewritten at 


» 


“.(k-l)“ 
x.(k-l) : 


a. 0 ;x,(ki; ;b, : 

i ,; I ” ! *iu(k) 
3 a,' : x, (k M b , , 


(27) 


y(k) • J. c.;x. (k) 


(28) 
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